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Abstract 

Grain  boundaries  exhibit  phase-like  behavior  in  which  their  structure,  chemistry  and  properties  may  change  discontinuously  at  critical 
values  of  thermodynamic  parameters  such  as  temperature,  pressure  and  chemical  potential.  Therefore,  grain  boundaries  (and  other  inter¬ 
faces  such  as  surfaces  and  heterophase  boundaries)  can  be  treated  as  thermodynamically  stable  interfacial  states.  To  differentiate  these 
interfacial  states  from  bulk  phases,  the  term  “complexion”  has  been  introduced.  A  variety  of  terminology  has  been  used  to  describe  com¬ 
plexions  and  complexion  transitions.  In  many  cases,  several  terms  exist  that  describe  essentially  the  same  phenomenon.  We  give  an  over¬ 
view  of  complexion-related  terminology,  suggest  a  preferred  nomenclature  and  discuss  a  classification  framework  that  can  be  used  to 
categorize  complexions  and  complexion  transitions.  The  field  of  grain  boundary  complexions  has  evolved  rapidly  in  the  past  decade 
due  to  advances  in  experimental  equipment  -  in  particular,  aberration-corrected  transmission  electron  microscopy  -  and  progress  in  com¬ 
putational  simulation  methods.  Grain  boundary  complexion  transitions  are  the  root  cause  of  a  wide  variety  of  materials  phenomena  - 
such  as  abnormal  grain  growth,  grain  boundary  embrittlement  and  activated  sintering  -  that  have  defied  mechanistic  explanation  for 
years.  In  this  overview,  we  review  the  history  and  theory  of  grain  boundary  complexion  transitions,  their  role  in  materials  processing 
and  their  effect  on  materials  properties. 

©2013  Acta  Materialia  Inc.  Published  by  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 

Grain  boundaries  strongly  influence  the  properties  and 
behavior  of  polycrystalline  materials  during  processing 
and  in  service  [1],  Grain  boundary  energy,  mobility,  diffu- 
sivity,  cohesive  strength  and  sliding  resistance  -  all  of 
which  depend  upon  local  structure  and  chemistry  -  deter¬ 
mine  bulk  materials  behavior  and  properties  such  as 
superplasticity,  creep,  fatigue,  corrosion,  strength  and 
conductivity  [2],  Grain  boundary  specialists  have  long 
recognized  that  grain  boundaries  can  be  considered  as 
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quasi-two-dimensional  “phases”  that  may  undergo  phase¬ 
like  transitions  in  which  their  structure  and  chemistry 
changes  abruptly  at  critical  values  of  thermodynamic 
parameters  [1,3-11],  As  predicted  by  Hart  [5],  the  grain 
boundary  properties  most  strongly  affected  by  these  transi¬ 
tions  will  be  the  non-equilibrium  properties,  such  as  mobil¬ 
ity,  diffusivity,  intergranular  cohesive  strength  and  grain 
boundary  sliding  resistance.  Equilibrium  grain  boundary 
properties,  such  as  energy,  specific  volume  and  adsorption, 
will  also  be  affected,  but  in  general  are  less  accessible  to 
experiment  than  non-equilibrium  grain  boundary  proper¬ 
ties.  If  a  significant  fraction  of  the  grain  boundaries  in  a 
polycrystalline  material  undergo  a  transition,  the  cumula¬ 
tive  effect  can  be  a  dramatic  and  rather  sudden  change  in 
macroscopic  properties  [5], 
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Grain  boundary  complexion  transitions  have  tradition¬ 
ally  been  called  grain  boundary  “phase  transitions”,  in 
analogy  to  bulk  phase  transitions.  Interfaces  such  as  grain 
boundaries  can  be  analyzed  using  equilibrium  thermody¬ 
namics  [12]  and  may  transform  in  a  manner  analogous  to 
bulk  phase  transitions,  but  equilibrium  interfacial  states 
do  not  satisfy  the  Gibbs  definition  of  a  phase  because  they 
are  inhomogeneous  and  may  have  gradients  of  structure, 
composition,  properties  and  other  order  parameters. 
Therefore,  Tang  et  al.  [13]  introduced  the  term  “complex¬ 
ion”2  to  denote  an  equilibrium  interfacial  state.  A  complex¬ 
ion,  concisely  defined,  is  interfacial  material  or  strata  that 
is  in  thermodynamic  equilibrium  with  the  abutting  phase(s) 
and  has  a  stable,  finite  thickness  that  is  typically  on  the 
order  of  0.2-2  nm  [14],  A  complexion  cannot  exist  indepen¬ 
dently  of  the  abutting  phases  and  its  average  composition 
and  structure  need  not  be  the  same  as  the  abutting  phases. 

com-plex-ion  noun  Xlom-'plek-shanX 

Interfacial  material  or  strata  that  is  in  thermodynamic 
equilibrium  with  its  abutting  phase(s) 


The  term  “complexion”  has  been  adopted  by  the 
authors  of  this  overview  [11,15-36],  and  has  been  recog¬ 
nized  [37-41]  and  adopted  [42-45]  by  others.  Therefore, 
in  this  overview  the  term  “grain  boundary  complexion” 
will  be  used  rather  than  “grain  boundary  phase”  to 
describe  an  equilibrium  grain  boundary  state,  even  though 
the  original  “phase”  terminology  may  have  been  used  by 
the  original  author(s).  Similarly,  although  the  “surface 
phase”  terminology  is  widely  used  in  the  surface  science 
community,  equilibrium  interfacial  states  at  free  surfaces 
will  be  referred  to  here  as  “surface  complexions”. 

Grain  boundary  complexion  transitions  are  difficult  to 
predict  a  priori  and  typically  occur  independently  of  bulk 
phase  transformations,  i.e.  at  different  values  of  tempera¬ 
ture,  pressure  or  composition.  Therefore,  they  often  result 
in  unexpected  and  seemingly  inexplicable  phenomena.  In 
1968,  Hart  postulated  that  grain  boundary  complexion 
transitions  could  be  responsible  for  temper  embrittlement 
in  steel  [3-5],  although  he  recognized  that  insufficient 
experimental  evidence  existed  at  the  time  to  support  the 
hypothesis  [4],  It  was  not  until  20  years  after  Hart’s  first 
paper  that  the  existence  of  grain  boundary  complexion 
transitions  was  deemed  to  be  conclusive  [8],  Since  then, 
advances  in  transmission  electron  microscopy  and  other 
experimental  techniques  have  provided  strong  evidence 
that  grain  boundary  complexion  transitions  are  responsible 
for  materials  phenomena  as  diverse  as  abnormal  grain 
growth  [17,46],  solid-state  activated  sintering  [47]  and 
liquid-metal  embrittlement  [48], 

Interfaces  can  be  separated  into  two  groups:  external 
interfaces  (e.g.  surfaces)  and  internal  interfaces  (e.g.  grain 
boundaries,  heterophase  boundaries,  stacking  faults  and 


2  The  term  “complexion”  was  tentatively  suggested  by  W.C.  Carter  and 
embraced  by  R.M.  Cannon. 


antiphase  boundaries)  [9],  All  of  these  interfaces  play  an 
important  role  in  materials  science,  and  understanding 
their  phase-like  behavior  has  been  identified  as  one  of  eight 
grand  challenges  in  ceramic  science  [49],  In  this  overview, 
we  will  focus  on  grain  boundary  complexions  and  their 
impact  on  materials  properties  and  processing.  We  will  also 
occasionally  discuss  surface  complexions  and  heterophase 
boundary  complexions  when  appropriate  to  illustrate 
important  concepts. 

2.  The  fundamentals  of  grain  boundary  complexions 

Seminal  articles  on  the  fundamentals  of  grain  boundary 
transitions  include  those  by  Hart  [3,5],  Cahn  [6]  and  Rott- 
man  [8,9].  Grain  boundary  transitions  are  also  discussed  in 
the  textbook  on  crystalline  interfaces  by  Sutton  and  Balluffi 
[1],  In  this  section,  we  will  summarize  the  main  points  and 
conclusions  of  these  thermodynamic  treatments. 

While  the  vast  majority  of  research  on  grain  boundary 
complexions  has  focused  on  non-pure  systems  in  which  grain 
boundary  segregation  plays  a  critical  role,  grain  boundary 
complexion  transitions  may  occur  in  even  pure  materials. 
We  will  first  discuss  grain  boundary  complexions  in  pure 
materials  to  illustrate  their  fundamental  characteristics. 
Then  we  will  discuss  the  more  complicated  complexion  tran¬ 
sitions  in  non-pure  systems  that  can  involve  changes  in  grain 
boundary  composition  as  well  as  structure. 

2.1.  Thermodynamic  parameters  of  grain  boundaries 

The  fundamental  thermodynamic  quantity  that  charac¬ 
terizes  grain  boundaries  (and  interfaces  in  general)  is  the 
energy  per  unit  area,  y,  which  represents  the  amount  of 
work  required  to  create  one  unit  of  grain  boundary  area 
via  a  reversible  process.  In  other  words,  y  represents  the 
additional  (excess)  free  energy  per  unit  area  that  exists  in 
the  system  due  to  the  presence  of  the  grain  boundary. 
The  value  of  y  is  a  function  of  both  bulk  and  interfacial 
thermodynamic  parameters.  The  usual  bulk  thermody¬ 
namic  parameters  are  temperature  (T),  pressure  (P)  and/ 
or  chemical  potential  (pi)  (or,  alternatively,  composition 
[6]).  There  are  five  interfacial  thermodynamic  parameters 
that  describe  the  five  macroscopic  degrees  of  freedom  of 
a  grain  boundary.  Several  different  geometric  conventions 
exist  for  describing  these  five  parameters  [50],  For  ease  of 
conceptual  understanding,  we  use  the  simple  convention 
followed  by  Rottman  [9],  in  which  three  variables  specify 
the  misorientation  vector  R  between  the  two  crystals  and 
two  independent  variables  in  the  form  of  a  unit  vector,  it, 
specify  the  average  orientation  of  the  grain  boundary  plane 
normal  (i.e.  the  grain  boundary  plane  inclination).  The 
misorientation  vector  R  is  defined  by  a  rotation  of  angle 
m  around  a  direction  common  to  both  grains  specified  by 
a  unit  vector  «,  i.e.  R  =  urn.  Given  a  misorientation  vector 
R,  the  rotation  angle  and  axis  can  be  recovered  using  the 
relations  co  =  \R\  and  it  =  R/co. 

The  energy  of  a  grain  boundary  defined  using  these 
parameters  depends  upon  6  +  C  independent  variables, 
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which  includes  the  five  geometric  parameters  of  the  grain 
boundary,  where  C  is  the  number  of  components  [6].  The 
grain  boundary  energy  y  may  also  depend  upon  applied 
electric  and  magnetic  fields,  for  example,  but  their  influence 
on  complexion  stability  and  transitions  has  not  been  trea¬ 
ted  in  detail  and  will  not  be  discussed  in  this  overview. 

In  addition  to  these  6  +  C  variables,  there  are  three 
more  geometric  parameters  that  specify  the  relative  transla¬ 
tional  displacement  of  the  two  crystals  and  an  additional 
geometric  parameter  that  specifies  the  translational  dis¬ 
placement  of  the  grain  boundary  plane  along  its  normal 
direction  with  respect  to  the  two  crystal  lattices  [9].  How¬ 
ever,  in  most  analyses,  the  grain  boundary  is  assumed  to 
relax  to  equilibrium  with  respect  to  these  four  microscopic 
parameters,  and  thus  they  are  typically  ignored  [6,9], 

The  grain  boundary  energy,  y,  is  always  a  continuous 
function  of  one  of  the  thermodynamic  parameters  (e.g.  T, 
P,  R  or  h)  for  a  complexion  transition  that  occurs  under 
equilibrium  conditions.  This  implies  that  two  (or  more) 
complexions  may  coexist  at  equilibrium  in  a  manner  anal¬ 
ogous  to  the  coexistence  of  bulk  phases  at  solidus  and  liq- 
uidus  lines  on  phase  diagrams. 

A  first-order  complexion  transition  occurs  when  there  is 
a  discontinuity  in  the  first  derivative  of  y  as  a  function  of 
one  of  the  thermodynamic  parameters  (T,  P,  q,-,  R  or  «). 
A  second-order  (or  higher-order)  complexion  transition 
occurs  when  there  is  a  discontinuity  in  a  second  (or 
higher-order)  derivative  of  y  [8,9].  This  thermodynamic 
definition  of  a  complexion  transition  is  fundamental  and 
universal,  and  is  directly  analogous  to  the  definition  of  a 
bulk  phase  transition  [9],  It  is  the  key  to  recognizing  when 
one  complexion  is  “thermodynamically  distinct”  from 
another  complexion,  in  the  same  way  that  two  bulk  phases 
can  be  recognized  as  thermodynamically  distinct  from  each 
other. 

We  emphasize  that  complexion  transitions  are  not  lim¬ 
ited  to  changes  in  the  atomic  structure  or  thickness  of  the 
grain  boundary  core  (which  are  sometimes  called  “struc¬ 
tural  transitions”)  but  may  also  include  transitions  in  com¬ 
position,  chemical  bonding,  roughening,  atomic 
reconstructions  of  the  grain  boundary  core  even  if  they 
do  not  result  in  a  thickness  change,  and,  indeed,  any  phys¬ 
ical  change  that  results  in  a  discontinuity  of  the  first-  or 
higher-order  derivative  of  y  as  a  function  of  one  of  the  ther¬ 
modynamic  parameters  (  T,  P.  ph  R  or  /*). 

Unfortunately,  in  real-world  experiments,  it  is  extraordi¬ 
narily  difficult  to  detect  complexion  transitions  based  on 
this  thermodynamic  definition,  so  approximations  are 
made  and  observable  features  such  as  grain  boundary 
thickness  and  solute  adsorption  are  used  as  proxies  to  iden¬ 
tify  distinct  complexions.  The  use  of  such  proxies  is  not  an 
infallible  method  and  may  sometimes  lead  to  two  complex¬ 
ions  being  identified  as  distinct  from  one  another  when 
they  are  in  fact  not  thermodynamically  distinct.  Theoreti¬ 
cal  models  and  computer  simulations  are  more  amenable 
to  applying  the  rigid  thermodynamic  definition  of  a  com¬ 
plexion  transition  to  identify  distinct  complexions. 


Fig.  1.  An  early  theoretical  grain  boundary  complexion  diagram  for  a 
pure  material  showing  the  pressure-temperature  ( P-T)  locus  of  a 
complexion  transition  between  two  hypothetical  complexions,  complexion 
a  (blue  region)  and  complexion  (S  (red  region).  Hart  [5]  referred  to  ft  as  the 
“high  temperature  grain  boundary  phase  (H-phase)”  and  a  was  referred  to 
as  the  “low  temperature  grain  boundary  phase  (L-phase)”.  Adapted  from 
Ref.  [5]  (For  interpretation  of  the  references  to  color  in  this  figure  legend, 
the  reader  is  referred  to  the  web  version  of  this  article.). 

2.2.  Grain  boundary  complexions  in  pure  materials 

Hart  [5]  conducted  a  theoretical  thermodynamic  analy¬ 
sis  of  a  grain  boundary  with  fixed  geometry  in  a  pure  mate¬ 
rial  and  developed  a  Clausius-Clapeyron  equation  that 
defined  the  pressure-temperature  locus  of  a  grain  bound¬ 
ary  complexion  transition.  This  analysis  resulted  in  a  grain 
boundary  complexion  diagram  (Fig.  1)  showing  regions  of 
complexion  stability  for  two  complexions  in  P—T  space  and 
the  complexion  transition  line  dividing  these  two  regions. 
Hart  further  showed  that,  during  a  first-order  complexion 
transition,  the  magnitude  of  the  slope  discontinuity  in  y 
as  a  function  of  temperature  or  pressure  can  be  directly 
related  to  a  discontinuity  in  the  specific  excess  volume 
Avy  and  entropy  A.v.,  of  the  grain  boundary  as  follows 
(where  Avy  and  A sy  represent  the  difference  in  specific  vol¬ 
ume  and  entropy  between  the  two  complexions  coexisting 
at  equilibrium)  [5]: 

A  (dy/dT)p  =  -As,,  (1) 

A(dy/dP)T  =  Avy  (2) 

In  other  words,  the  slope  of  y-T  and  y-P  curves  (at  con¬ 
stant  pressure  and  temperature,  respectively)  will  be  differ¬ 
ent  above  and  below  the  complexion  transition  point,  as 
shown  schematically  in  Fig.  2.J  There  is  no  discontinuity 
in  the  value  of  y  itself  at  a  complexion  transition  point, 
which  follows  from  equilibrium  thermodynamics  [6].  As 
Hart  points  out,  the  type  of  first-order  complexion  transi¬ 
tion  shown  in  Fig.  2  is  associated  with  a  latent  heat  of 
transformation,  T  As..,  which  in  principle  can  be  measured 


3  An  apparent  drafting  error  led  to  a  reversal  of  slopes  in  the  original 
plot  of  y  vs.  P  in  Ref.  [5],  which  has  been  corrected  in  Fig.  2(a). 
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Fig.  2.  A  schematic  illustration  of  grain  boundary  excess  free  energy  y  in  a 
pure  material  as  a  function  of  (a]  pressure  and  (b)  temperature  for  two 
“competing”  complexions:  a  (blue  line)  and  /?  (red  line).  Solid  lines 
represent  the  stable  complexion  and  dashed  lines  represent  the  metastable 
complexion.  A  first-order  grain  boundary  complexion  transition  occurs  at 
the  vertical  dotted  black  lines  and  is  recognizable  by  the  discontinuity  in 
the  slope  of  y.  (a)  Complexion  /I  has  a  larger  excess  specific  volume  than  a. 
(and  hence  a  larger  slope,  (dy/SP)T),  and  is  therefore  the  stable  complexion 
at  low  pressure,  (b)  Complexion  fl  has  a  larger  excess  specific  entropy  than 
a  (and  hence  a  larger  negative  slope,  (6y/ST)P),  and  is  therefore  the  stable 
complexion  at  high  temperature.  Hart  referred  to  (1 ,  which  has  larger 
excess  specific  volume  and  larger  excess  specific  entropy,  as  the  “high 
temperature  phase  (H-phase)”,  and  a  was  referred  to  as  the  “low 
temperature  phase  (L-phase)”  [5],  Adapted  from  Hart  [5],  (For  interpre¬ 
tation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred 
to  the  web  version  of  this  article.) 

experimentally,  and  in  fact  has  been  measured  for  complex¬ 
ion  transitions  at  liquid-liquid  interfaces  [51,52].  It  is  also 
associated  with  a  change  in  the  specific  excess  volume  of 
the  grain  boundary,  Avy,  which  in  principle  can  be  detected 
by  sensitive  dilatometry  measurements  [5], 

Conceptually,  we  can  understand  the  mathematical 
behavior  of  y  during  a  first-order  complexion  transition 
by  envisioning  that  several  different  metastable  grain 
boundary  complexions  could  potentially  exist  at  a  given 
grain  boundary.  Each  metastable  complexion  may  have  a 
different  atomic  structure  (and,  for  multicomponent  sys¬ 
tems,  may  also  have  a  different  composition)  and  therefore 
a  different  value  of  y.  At  thermodynamic  equilibrium,  the 
complexion  with  the  smallest  value  of  y  is  the  stable  com¬ 
plexion  and  will  be  present  in  lieu  of  the  other  possible 
(metastable)  complexions.  As  any  given  thermodynamic 
parameter  is  varied,  the  y  curve  of  a  previously  metastable 
complexion  may  cross  that  of  the  stable  complexion, 
thereby  forcing  a  complexion  transition.  This  behavior  is 
illustrated  in  Fig.  2  for  complexions  a  and  ft.  The  a  com¬ 
plexion  is  stable  at  low  temperature  and  high  pressure, 
while  the  [>  complexion  is  stable  at  high  temperature  and 
low  pressure.  There  is  a  discontinuity  in  the  first  derivative 
of  y  at  the  complexion  transition  point  denoted  by  the  ver¬ 
tical  dotted  lines,  but  y  is  continuous  at  the  complexion 
transition  point. 


The  discussion  thus  far  has  focused  on  grain  boundary 
complexion  transitions  in  pure  materials  with  fixed  grain 
boundary  shape  and  geometrical  parameters  R  and  h.  Such 
complexion  transitions  are  called  “congruent”  transitions 
[6]  and  are  characterized  by  changes  in  atomic  structure 
in  the  grain  boundary  core.  At  equilibrium,  two  complex¬ 
ions,  a  and  /?,  may  coexist,  as  shown  schematically  in 
Fig.  3(a).  Grain  boundary  complexion  transitions  that 
result  in  a  change  in  the  geometrical  parameters  of  the 
grain  boundary  (the  “grain  boundary  character”)  are  called 
“non-congruent”  transitions  [6]  and  can  be  categorized  as 
either  faceting  transitions  or  dissociation  transitions.  In  a 
faceting  transition,  the  grain  boundary  plane  normal  it 
decomposes  into  h\  and  h2,  the  area-weighted  average  of 
which  is  equivalent  to  h.  The  coexistence  of  a  and  [l  during 
a  faceting  transition  is  shown  schematically  in  Fig.  3(b).  In 
a  dissociation  transition,  a  single  grain  boundary  dissoci¬ 
ates  into  two  new  interfaces,  separated  by  a  bulk  phase, 
with  misorientation  R\  and  R2,  the  combined  misorienta- 
tion  of  which  is  equal  to  the  original  misorientation  R 
(Fig.  3(c)).  The  dissociation  transition  is  more  commonly 
known  as  a  “wetting  transition”  and  occurs  when  the 
sum  of  the  energy  of  the  two  new  grain  boundaries  is  less 
than  that  of  the  original  grain  boundary.  If  the  wetting 
phase  is  not  crystalline,  then  R\  and  R>  are  undefined, 
but  the  original  misorientation  R  is  still  maintained  across 
the  phase. 

Cahn  [6]  discussed  the  geometry  of  grain  boundary  com¬ 
plexion  transitions  in  great  detail  and  reached  several 
important  conclusions,  which  we  summarize  here:  (i)  a 
smoothly  curving  grain  boundary  can  be  considered  as  a 
single  grain  boundary  complexion;  (ii)  grain  boundary 
complexion  transitions  usually  occur  at  “singularities”  in 
grain  boundary  geometry,  e.g.  at  edges  and  facets;  (iii) 
complexion  coexistence  and  transitions  are  almost  always 
non-congruent  and  thus  involve  discontinuities  in  grain 
boundary  geometry  (faceting  or  dissociation);  (iv)  congru¬ 
ent  complexion  transitions  are  possible  but  rare;  (v)  the 
observation  of  a  congruent  complexion  transition  might 
be  an  artifact  of  the  transition  occurring  far  from  equilib¬ 
rium  conditions.  For  further  details,  the  reader  is  referred 
to  the  original  text  [6].  However,  recent  experiments  and 
simulations  have  shown  that  congruent  complexion  transi¬ 
tions  are  not  as  uncommon  as  Cahn  suggested;  some  exam¬ 
ples  are  discussed  in  the  following  paragraph.  It  should 
also  be  pointed  out  that  Cahn  made  the  above  statements 
with  pure  systems  in  mind.  On  the  other  hand,  essentially 
all  materials  contain  impurities,  and  examples  of  congruent 
transitions  are  much  more  abundant  in  non-pure  systems, 
such  as  the  complexion  transitions  from  a  clean  complex¬ 
ion  to  a  monolayer,  bilayer  or  trilayer  complexion  (see  Sec¬ 
tion  2.3.3).  Therefore,  while  faceting  and  dissociation  are 
two  important  markers  for  grain  boundary  complexion 
transitions  in  pure  materials,  complexion  transitions  in 
general  entail  a  rich  variety  of  behaviors,  including  those 
that  maintain  the  grain  boundary  geometry. 
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Fig.  3.  Schematic  illustration  of  the  coexistence  of  two  grain  boundary  complexions,  a  and  /i,  at  three  different  grain  boundary  geometries:  (a)  coexistence 
of  two  complexions  (on  the  same  macroscopic  plane)  at  grain  boundary  with  the  same  misorientation  R  and  inclination  n;  (b)  coexistence  of  two 
complexions  with  the  same  grain  boundary  misorientation  R  but  different  grain  boundary  inclinations,  hl  and  hi.  resulting  in  a  faceted  structure;  (c) 
coexistence  of  two  complexions  (separated  by  a  bulk  phase)  with  the  same  grain  boundary  inclination  it  but  different  grain  boundary  misorientations  R\ 
and  R2 ,  as  might  exist  after  a  dissociation  transition  (i.e.  a  wetting  transition).  Complexion  coarsening  may  occur  in  (a)  and  (b)  to  reduce  the  total  length  of 
the  one-dimensional  complexion  boundaries.  Adapted  front  Ref.  [9], 


Despite  decades  of  research,  efforts  to  identify  grain 
boundary  complexion  transitions  in  pure  metals  via  high- 
resolution  transmission  electron  microscopy  (HRTEM) 
and  computer  simulations  have  been  largely  unsuccessful. 
Nevertheless,  experimental  evidence  suggesting  grain 
boundary  complexion  transitions  in  pure  metals  has  existed 
for  decades.  For  example,  researchers  have  reported 
anomalies  and  discontinuities  in  the  activation  energy  as 
a  function  of  temperature  for  grain  boundary  mobility  in 
zone-refined  lead  [53],  grain  boundary  diffusivity  in  pure 
copper  [54]  and  grain  boundary  sliding  in  zinc  bicrystals 
[55].  A  fair  criticism  of  these  types  of  experimental  studies 
is  that  they  could  be  influenced  by  trace  impurities. 
However,  a  recent  molecular  dynamics  (MD)  computer 
simulation  employing  an  innovative  methodology  has 
successfully  demonstrated  a  grain  boundary  complexion 
transition  in  pure  Cu  grain  boundaries  at  temperatures  well 
below  the  melting  point,  providing  strong  support  for  the 
existence  of  grain  boundary  transitions  in  pure  materials 
[39].  The  transition  demonstrated  by  the  MD  simulation 
was  a  congruent  transition,  in  which  the  atomic  structure 
of  the  grain  boundary  core  changed  but  the  grain  boundary 
geometry  (character)  remained  invariant.  The  coexistence 
of  these  two  intrinsic  grain  boundary  complexions  is  shown 
in  Fig.  4(a)  and  the  atomic  structure  of  each  complexion  is 
shown  in  Fig.  4(b  and  c).  Furthermore,  the  complexion 
transition  exhibited  a  discontinuity  in  grain  boundary 
density,  consistent  with  Hart’s  analysis  [5].  The  authors 
suggest  that  this  type  of  grain  boundary  complexion  tran¬ 
sition  may  play  a  role  in  the  healing  of  radiation  damage, 
and  that  grain  boundary  complexion  transitions  in  pure 
metals  may  be  more  common  than  previously  thought  [39]. 

2.3.  Grain  boundary  complexions  in  non-pur e  materials 

Grain  boundary  adsorption  (segregation  )  may  occur  in 
non-pure  materials,  causing  the  composition  of  grain 
boundaries  to  differ  from  the  overall  composition  at  equi¬ 
librium  [2],  Differences  in  grain  boundary  composition 
may  also  arise  from  kinetic  processes  such  as  grain  growth. 


4  The  terminologies  “grain  boundary  adsorption”  and  “grain  boundary 
segregation”  are  used  interchangeably  to  refer  to  the  same  phenomenon. 


Such  non-equilibrium  segregation  is  referred  to  as  enrich¬ 
ment  and  will  not  be  discussed  in  detail  here. 

Grain  boundary  adsorption  plays  a  dominant  role  in 
complexion  transitions  in  non-pure  systems.  Adsorption 
phenomena  must  also  be  considered  when  conducting 
experiments  on  nominally  pure  systems,  which  invariably 
contain  impurities.  When  Hart  first  suggested  the  concept 
of  grain  boundary  complexion  transitions  in  1968,  he  pro¬ 
posed  that  they  might  be  responsible  for  temper  embrittle¬ 
ment  of  steel  via  a  first-order  discontinuity  in  the  grain 
boundary  adsorption  of  solute  elements  [3].  Hart  made  this 
prediction  based  on  an  analogy  to  two-dimensional  phase 
transitions  on  free  surfaces,  a  concept  already  well  known 
at  the  time  from  the  work  of  Fowler  and  Guggenheim 
[56]  and  others  [57,58], 

Early  theories  of  grain  boundary  adsorption  were  based 
on  the  McLean  model,  which  assumes  that  segregation  is 
constrained  to  a  monolayer  of  atomic  sites  at  the  grain 
boundary  [59].  Later  models  acknowledged  that  multilayer 
grain  boundary  segregation  can  occur  and  may  be  accom¬ 
panied  by  structural  changes  in  the  grain  boundary  core 
[13,60,61].  When  multilayer  segregation  occurs,  entirely 
new  structures  may  manifest  at  the  grain  boundary  and 
induce  a  grain  boundary  complexion  transition.  Recent 
experimental  work  has  demonstrated  that  such  adsorp¬ 
tion-based  grain  boundary  complexion  transitions  are 
linked  to  phenomena  such  as  abnormal  grain  growth  in 
ceramics  [17],  activated  sintering  [47]  and  liquid  metal 
embrittlement  [48],  Various  terms  have  been  used  to 
describe  multilayer  adsorption  and  structural  transitions 
at  grain  boundaries.  We  will  direct  most  of  our  focus  on 
premelting  and  prewetting  transitions,  which  are  two  of 
the  most  widely  studied  types  of  grain  boundary  complex¬ 
ion  transitions. 

In  this  section,  we  discuss  the  fundamentals  of  grain 
boundary  complexion  transitions  in  non-pure  systems,  with 
an  emphasis  on  theoretical  models  for  grain  boundary 
adsorption.  We  first  review  the  conditions  necessary  for  a 
first-order  complexion  transition  based  on  classical  models 
of  grain  boundary  adsorption,  then  move  onto  more  mod¬ 
ern  models  of  premelting  and  prewetting  transitions. 
Finally,  we  discuss  the  fundamentals  of  discrete 
complexions. 
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Fig.  4.  MD  and  Monte  Carlo  simulation  results  that  demonstrate  intrinsic  (a-c)  and  extrinsic  (d,  e)  grain  boundary  complexion  transitions  at  a  £5(310) 
boundary  in  pure  Cu  (a-c)  and  Ag-doped  Cu  (d,  e).  A  congruent  complexion  transition  is  shown  in  (a),  in  which  the  two  intrinsic  grain  boundary 
complexions  (“kites”  and  “split  kites”)  coexist  and  meet  at  a  one-dimensional  complexion  boundary.  The  atomic  structure  of  the  intrinsic  “kites” 
complexion  is  shown  in  (b)  and  the  atomic  structure  of  the  “split  kites”  complexion  is  shown  in  (c).  When  the  Cu  is  doped  with  Ag,  a  similar  complexion 
transition  occurs  from  (d)  “kites”,  with  a  monolayer  complexion  with  Ag  adsorption,  to  (d)  “split-kites”,  with  bilayer  complexions  with  Ag  adsorption. 
These  two  complexions  can  be  considered  as  two  different  discrete  Dillon-FIarmer  complexions  (i.e.  monolayer  and  bilayer  complexions).  The 
complexions  in  (b)  and  (c)  are  two  manifestations  of  the  clean  complexion  in  the  Dillon-Harmer  scheme.  Parts  (a-c)  reprinted  from  Ref.  [39]  with 
permission  and  parts  (d,  e)  reprinted  from  Ref.  149]  with  permission. 


We  cannot  hope  to  give  a  fully  rigorous  and  complete 
account  of  all  theories  and  models  used  to  describe  these 
phenomena,  which  are  numerous  and  highly  detailed. 
More  complete  accounts  can  be  found  in  various  textbooks 
[1,41,56,59,62]  and  review  articles  [2,15,22,23,32,63],  and 
the  original  research  articles  referenced  therein.  Our  aim 
here  is  to  give  an  overview  of  grain  boundary  complexions 
and  transitions  in  non-pure  systems,  with  an  emphasis  on 
providing  physical  insight  and  understanding. 

2.3.1.  Classical  grain  boundary  adsorption! segregation 
models 

In  1957,  McLean  [59]  proposed  a  classical  grain  bound¬ 
ary  adsorption  model  via  an  analogy  of  the  Langmuir  sur¬ 
face  adsorption  model  [64],  This  Langmuir-McLean  model 
can  be  derived  with  the  following  assumptions: 

(1)  There  are  a  fixed  number  of  (homogeneous)  excess 
adsorption  sites  per  unit  area  at  a  grain  boundary. 


(2)  The  grain  boundary  can  be  treated  as  a  separate  ther¬ 
modynamic  entity  that  is  in  equilibrium  with  the  bulk 
phase  (with  a  given  fraction  of  solute  of  Xbldk). 

(3)  Both  the  bulk  phase  and  the  grain  boundary  can  be 
treated  as  ideal  solutions. 

The  grain  boundary  atomic  fraction  of  solute  XGB  is 
then  given  by  the  following  Langmuir-McLean  isotherm: 

XqB  Xhutk  J^qds 

i  x  e  (3) 

where  R  is  the  universal  gas  constant  and  A Gads  is  the  stan¬ 
dard  molar  Gibbs  free  energy  of  adsorption,  which  is  de¬ 
fined  as  a  negative  value  for  positive  adsorption  of  solute 
at  the  grain  boundary. 

Great  efforts  have  been  made  to  refine  and  extend  the 
Langmuir-McLean  grain  boundary  adsorption  model  in 
various  aspects;  readers  are  referred  to  elegant  articles  by 
Hondros  and  Seah  [65],  Wynblatt  and  Chatain  [66]  and 
Lejcek  and  Hofmann  [67]  for  elaboration  of  these  grain 
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boundary  adsorption  models.  Perhaps  the  simplest  way  to 
illustrate  the  possible  occurrence  of  a  first-order  grain 
boundary  complexion  transition  is  by  analyzing  the  gener¬ 
alized  Langmuir-McLean  isotherm  that  was  given  in  Ref. 

[67]: 


nM 

UGB 


1 bulk 
~M  ‘ 
lbulk 


(4) 


where  and  a^ulk  are  the  activities  of  the  impu¬ 

rity  (i.e.  solute,  I)  and  matrix  (M)  elements  in  the  grain 
boundary  (GB)  and  bulk  phase,  respectively.  It  is  worth 
noting  that  this  formulism  [67]  was  derived  by  treating 
the  grain  boundary  as  a  stand-alone  “phase”  instead  of  a 
complexion  by  ignoring  the  gradient  energy  term  associ¬ 
ated  with  the  “cross”  bonds  between  the  grain  boundary 
atoms  and  matrix  atoms.  A  more  rigorous  treatment  that 
considers  the  multilayer  nature  and  compositional  gradi¬ 
ents  is  given  by  Wynblatt  and  Chatain  [66],  Nonetheless, 
a  special  case  of  this  less-rigorous  generalized  isotherm 
can  illustrate  the  possible  occurrence  of  a  first-order  grain 
boundary  complexion  transition  as  well  as  its  origin.  Here, 
we  assume  that  the  grain  boundary  can  be  represented  as  a 
regular  solution  (with  an  interaction  parameter,  QGB)  and 
the  bulk  phase  is  still  an  ideal  solution  (if  the  bulk  phase 
is  not  an  ideal  solution,  a  similar  result  can  still  be  obtained 
in  the  dilute-solution  limit  of  the  bulk  phase,  assuming 
Henry’s  law  applies).  The  generalized  Langmuir-McLean 
isotherm  in  Eq.  (4)  can  then  be  simplified  to: 


Xgb 
1  —XGb 


X bulk 
1  X  bulk 


&G^l~2aGBxGB 


(5) 


where  Ag££(=  AG*ads  +  QGB).  Eq.  (5)  is  essentially  equiva¬ 
lent  to  the  Fowler-Guggenheim  isotherm  [56],  which  is 
commonly  expressed  as: 


X  GB  Xbulk  _AGli+agyg 8 

- = - e  RT 

\—xGB  1  —Xbuik 


(6) 


where  z  the  coordination  number  and  a  (=—2QGB/z)  is  a 
parameter  characterizing  adsorbate-adsorbate  interac¬ 
tions.  The  Fowler-Guggenheim  isotherm  was  also  first 
proposed  for  surface  adsorption  and  then  adapted  for 
grain  boundary  adsorption  [3,65],  In  a  classic  article  in 
1977,  Hondros  and  Seah  [65]  pointed  out  that  this  Fow¬ 
ler-Guggenheim  isotherm  can  produce  a  first-order  grain 
boundary  complexion  transition  if  the  adsorbate-adsor¬ 
bate  interaction  is  strongly  attractive  (under  the  condition 
that  2 Qgb  =  —a z  >  ~4 RT)\  this  is  illustrated  in  Fig.  5(a).  It 
is  interesting  to  note  that  this  condition  for  the  occurrence 
of  a  first-order  complexion  transition  is  identical  to  that  for 
a  phase  separation  of  a  regular  solution  (Q  >  2 RT),  which 
provides  a  clear  hint  of  the  physical  origin  of  this  first-order 
grain  boundary  complexion  transition. 

We  should  further  note  that  this  first-order  complexion 
transition  as  a  result  of  an  attractive  adsorbate-adsorbate 
interaction  occurs  because  of  the  same  physical  principle 
that  produces  the  first-order  gas-to-liquid  bulk  phase 


transition  in  the  van  der  Waals  isotherm  for  non-ideal 
gas  with  an  attractive  intermolecular  interaction 
(Fig.  5(b)).  In  fact,  in  two  original  articles,  published  in 
1968  and  1972  [3,5],  Hart  proposed  that  grain  boundaries 
can  undergo  first-order  complexion  transitions,  in  part 
based  on  analogy  to  the  three-dimensional  gas-liquid  bulk 
phase  transformation,  in  addition  to  the  analogy  to  surface 
phase  transitions  of  the  type  described  theoretically  by 
Fowler  and  Guggenheim  [56], 

We  point  out  that  the  Fowler-Guggenheim  model, 
when  applied  to  grain  boundary  adsorption,  in  essence 
demonstrates  a  first-order  complexion  transition  within 
a  grain  boundary  of  fixed  structural  width  (or,  strictly 
speaking,  of  fixed  number  of  adsorbate  sites)  in  which 
there  is  a  jump  from  low  adsorbate  content  to  high 
adsorbate  content.  With  more  advanced  models,  it  is  also 
possible  to  have  similar  complexion  transitions  between 
grain  boundaries  of  different  structural  widths  (and  hence 
of  different  numbers  of  adsorbate  sites),  in  which  the 
thickness  of  the  grain  boundary  core  changes  together 
with  the  amount  of  solute  adsorption.  The  ability  of  a 
boundary  to  undergo  structural  and  compositional  tran¬ 
sitions  within  a  boundary  core  of  fixed  width  or  between 
boundaries  of  different  fixed  widths  indicates  that  there 
could  be  many  more  or  many  fewer  complexions  in  a 
given  system  than  suggested  by  the  discrete  Dillon- 
Harmer  complexion  categories  that  are  introduced  later 
in  Section  2.3.3. 

The  above  discussion  clearly  demonstrates  an  analogy 
between  the  bulk  phase  transitions  and  interfacial  complex¬ 
ion  transitions.  Although  the  Fowler-Guggenheim  iso¬ 
therm  (i.e.  a  modified  Langmuir-McLean  model  where 
the  grain  boundary  behaves  like  a  regular  solution)  can 
produce  a  first-order  complexion  transition,  this  model 
has  several  limitations,  because  it  does  not  consider  that: 

•  the  grain  boundary  adsorption  can  occur  in  a  multilayer 
fashion  where  the  adsorption  sites  are  not  all  equivalent 

[66]; 

•  in  addition  to  chemical  adsorption,  grain  boundaries 
can  undergo  transitions  in  atomic  structure  that  will 
change  the  adsorption  free  energies  (i.e.  an  interfacial 
disordering  transition  can  be  coupled  with  an  adsorp¬ 
tion  transition  [16,60,61]);  and 

•  rigorously,  a  grain  boundary  cannot  be  treated  as  a  sep¬ 
arate,  homogeneous  thermodynamic  entity  (e.g.  even  if 
we  assume  that  adsorption  only  occurs  within  one 
monolayer  at  the  grain  boundary,  there  will  be  an  extra 
gradient  energy  term  that  will  depend  on  both  bulk  and 
grain  boundary  compositions)  [66], 

Thus,  more  realistic  models  should  consider  the  spatial 
gradients  in  both  composition  and  structure.  These  can 
be  better  treated  in  lattice-gas  models  [60,66,68]  or  dif¬ 
fuse-interface  (phase-field)  models  [13,16,61,69],  which  will 
be  discussed  in  the  following  sections. 
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Fig.  5.  Analogy  between  (a)  a  first-order  complexion  transition  in  which  there  is  a  discontinuous  jump  in  adsorbed  solute  content  and  (b)  a  first-order  bulk 
phase  transition  from  gas  to  liquid.  Both  of  these  first-order  transitions  arise  when  the  attractive  interactions  between  particles  (adsorbate  atoms  or  gas 
molecules,  respectively)  reach  a  critical  strength,  (a)  The  Fowler-Guggenheim  isotherm  is  plotted  for  increasing  values  of  2 QaB/RT  (=—azl(RT)).  When 
this  value  is  equal  to  zero,  Langmuir-McLean  behavior  occurs,  (b)  The  van  der  Waals  equation  of  a  non-ideal  gas  (P  =  RT/(Vm  —  h)  +  a/V is  plotted 
for  increasing  values  of  “a”,  where  “a”  is  the  strength  of  intermolecular  attraction  and  “b”  is  the  volume  excluded  by  a  mole  of  gas  molecules.  When  “a”  is 
equal  to  zero,  the  gas  is  ideal.  First-order  transitions  in  (a)  and  (b)  are  represented  by  vertical  lines  with  arrows. 


2.3.2.  Prewetting  and  premelting  complexion  transitions 

Prewetting  and  premelting  transitions  are  two  types  of 
complexion  transitions  that  have  been  studied  extensively 
during  the  past  few  decades.  Premelting  refers  to  the  for¬ 
mation  of  a  disordered,  liquid-like  film  on  a  crystalline  sur¬ 
face  (or  at  a  grain  boundary,  phase  boundary  or  free 
surface)  at  a  temperature  below  the  bulk  melting  tempera¬ 
ture  or  solidus  of  the  underlying  crystal.  Extensive  discus¬ 
sions  of  premelting  can  be  found  in  review  articles 
[70,71],  Strictly  speaking,  premelting  is  defined  only  for 
unary  systems  in  which  there  is  a  single  melting  tempera¬ 
ture.  However,  premelting  may  also  refer  to  high-tempera¬ 
ture  disordering  of  grain  boundaries  in  multicomponent 
systems. 

Prewetting  is  a  type  of  complexion  transition  that  occurs 
when  a  layer  of  material  of  fixed  equilibrium  thickness 
forms  at  the  interface  in  the  thermodynamic  vicinity  of  a 
wetting  transition,  i.e.  near  the  temperature  or  composition 
at  which  a  wetting  transition  would  occur.  The  discovery  of 
the  prewetting  transition  is  attributed  to  Calm  [69]  and 
Ebner  and  Saam  [72],  As  the  term  is  commonly  used,  it 
refers  to  a  first-order  adsorption  transition  in  which  the 
“film”  of  material  separating  two  phases  undergoes  a  dis¬ 
continuity  in  both  adsorption  and  thickness.  Premelting 
can  be  considered  as  a  special  case  of  prewetting  in  which 
there  is  a  discontinuity  in  crystalline  order  and  thickness 
but  not  necessarily  in  adsorption.  These  transitions  are 
strongly  related  to  wetting  transitions  and  hence  their  anal¬ 
ysis  is  rooted  in  a  wetting  analysis  that  compares  relative 
interfacial  energies  between  bulk  phases.  Therefore  we  will 
review  the  fundamentals  of  prewetting  and  premelting 
transitions  -  both  of  which  are  complexion  transitions  - 
in  the  context  of  wetting  transitions. 

When  a  liquid  droplet  is  placed  on  top  of  a  surface,  it 
will  partially  wet  the  surface  with  a  finite  contact  angle 


when  yh  +  ysi  >  ysv,  where  y/v  (yOT)  is  the  surface  energy  of 
the  liquid  (solid)  phase  and  ys/  is  the  liquid/solid  interface 
energy.  On  the  other  hand,  when  the  liquid  will  completely 
wet  the  surface  by  spreading  on  the  surface  to  form  a  con¬ 
tinuous  film,  which  has  an  arbitrary  thickness  that  is  only 
determined  by  the  volume  of  the  liquid  droplet,  we  have 
ysv  =  yi„  +  ys/  thermodynamically  (by  definition).  Prewett¬ 
ing  refers  to  a  transition  that  occurs  between  partial  and 
complete  wetting,  in  which  a  thin  layer  of  material  forms 
to  cover  the  surface  but  has  an  equilibrium  thickness  that 
is  controlled  by  thermodynamic  state  variables  (e.g.  tem¬ 
perature,  pressure)  and  does  not  depend  on  the  amount 
of  liquid  available.  The  prewetting  transition  was  first  pre¬ 
dicted  by  Cahn  in  his  seminal  paper  entitled  “Critical  Point 
Wetting”  [69],  In  this  paper,  Cahn  considers  a  binary  fluid 
with  a  miscibility  gap  that  is  in  contact  with  a  flat  container 
surface.  He  applies  a  diffuse-interface  model  to  describe  the 
free  energy  of  this  system,  which  in  one  dimension  is 
expressed  as 

/*oo 

F[c,  T]  =  <P(cs)  +  /  [f(c,T)  +  K(dc/dx)2]dx  (7) 

Jo 

where  <P  represents  a  short-range  surface  interaction  energy 
that  is  a  function  of  the  liquid  composition  at  the  surface, 
cs  =  c(x  =  0),  /(c  T)  is  the  homogeneous  free  energy  density 
of  the  liquid  and  K(dc/dx)2  is  the  gradient  energy  term  that 
penalizes  the  development  of  a  composition  gradient  within 
the  fluid.  For  the  purpose  of  discussion,  suppose  that  the 
surface  energy  with  the  container  has  a  larger  value  for  li¬ 
quid  phase  a  (rich  in  component  A),  yac,  than  for  phase  [3 
(rich  in  component  B),  ypc.  Mathematically,  this  is  reflected 
by  (l>  being  a  decreasing  function  of  c  and  having  a  mini¬ 
mum  at  c  =  1  (pure  B).  For  compositions  within  the  misci¬ 
bility  gap  (the  two-phase  coexistence  region),  Cahn  uses  a 
simple  scaling  argument  to  show  that  a  first-order  complete 
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Fig.  6.  Illustration  of  Cahn's  surface  wetting  diagram  in  (a)  the  composition-temperature  (c-T)  space  and  (b)  the  excess  chemical  potential-temperature 
(A  fi-T)  space,  where  A  n  =  n  —  )ic(T)  is  the  chemical  potential  of  B  (/i)  relative  to  its  value  at  two-phase  coexistence  which  is  a  function  of  temperature. 
The  prewetting  line  extends  into  the  single-phase  region  from  the  two-phase  region,  beginning  at  Twet  (the  bulk  wetting  temperature)  and  terminating  at  the 
surface  critical  point  (marked  by  a  red  dot  in  (a)  and  (b)).  The  bulk  critical  temperature  is  Tc  (green  dot  in  (a)  and  (b))  and  the  two-phase  coexistence 
condition  is  represented  in  (b)  by  a  single  vertical  line,  A/x  =  0,  in  A/i -T  space.  Adapted  from  Ref.  [15].  (For  interpretation  of  the  references  to  color  in  this 
figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


wetting  transition  will  occur  upon  raising  temperature  to¬ 
wards  the  critical  point  of  the  bulk  liquid.  Above  a  transi¬ 
tion  temperature  Twet,  ypc  +  yap  =  yac  and  phase  [3  will 
perfectly  wet  the  container  surface  and  separate  phase  a 
from  the  container,  as  illustrated  in  Fig.  6. 

For  compositions  in  the  single  a-phase  region,  the  com¬ 
petition  between  the  short-range  surface  interaction  fP. 
which  favors  /?,  and  the  homogeneous  free  energy  f{c,  T), 
which  prefers  a,  results  in  the  formation  of  a  thin,  B-rich 
surface  adsorption  layer.  When  the  bulk  fluid  composition 
reaches  the  solubility  limit  at  T  >  Twet,  the  adsorption  layer 
will  transform  into  bulk  j3  phase  and  completely  wet  the 
surface.  When  gradually  increasing  the  bulk  liquid  compo¬ 
sition  at  F>  Twet,  Cahn  uses  an  elegant  and  general  gra¬ 
phic  analysis  of  Eq.  (7)  to  show  that  a  first-order 
prewetting  transition  can  develop  in  the  single-phase  region 
before  the  adsorption  layer  becomes  bulk  [I  phase  and  per¬ 
fectly  wets  the  surface  at  the  solubility  limit.  The  conditions 
for  the  transition  are  represented  by  a  prewetting  line  cPn(- 

on  the  phase  diagram,  which  is  determined  by  solving 
the  equilibrium  equation  of  the  system.  As  shown  in 
Fig.  6(a),  the  prewetting  line  meets  the  solubility  line  at 
Twet  and  terminates  at  a  surface  critical  point  on  the  other 
end.  On  cP^{T),  two  surface  complexions  coexist  in  equi¬ 
librium  -  a  thicker  complexion  with  more  B  segregation 
level  and  a  thinner  complexion  with  less  B  segregation. 
Both  the  B  segregation  amount  and  the  layer  thickness  will 
undergo  a  discontinuous  jump  across  this  line  (Fig.  6), 
which  is  characteristic  of  a  first-order  transition.  The  prew¬ 
etting  transition  predicted  by  Cahn  has  been  confirmed  in 
liquid  helium  [73],  organic  [74]  and  liquid  metal  [75] 
systems. 

The  superposition  of  the  prewetting  line  onto  the  bulk 
phase  diagram  represents  perhaps  the  first  example  of  a 
complexion  diagram.  Such  diagrams  supply  information 


on  the  dependence  of  the  equilibrium  quantities  of  a  micro¬ 
structure  defect  (e.g.  interface),  such  as  width  and  compo¬ 
sition  on  thermodynamic  state  variables,  but  these 
microstructural  aspects  do  not  produce  additional  degrees 
of  freedom  that  would  modify  the  rules  governing  the 
topology  of  bulk  phase  diagrams.  We  discuss  the  creation 
of  grain  boundary  complexion  diagrams  later  in 
Section  2.5. 

Cahn’s  critical  point  wetting  theory  (CPWT)  not  only 
inspired  substantial  efforts  in  the  experimental  search  of 
wetting  transitions,  it  also  laid  the  groundwork  for  many 
subsequent  theoretical  analyses  of  prewetting  and  related 
phenomena.  As  we  shall  show  below,  many  of  the  models 
later  proposed  to  explain  analogous  phenomena  are  in  fact 
mathematically  isomorphic  to  CPWT.  CPWT  can  also  be 
applied  to  model  vapor  adsorption  on  an  inert  substrate 
near  the  liquid-vapor  critical  point  by  simply  replacing 
composition  c  with  the  density  in  Eq.  (7)  [63,76],  Here 
the  prewetting  transition  gives  rise  to  a  dense,  liquid-like’ 
adsorption  layer  at  the  interface  between  the  substrate 
and  the  bulk  vapor  phase,  and  the  model  predictions  on 
the  conditions  for  wetting  transitions  are  identical  to  those 
given  by  Fig.  6. 

Wetting  transitions  can  also  occur  at  interfaces  in  solid 
materials  such  as  grain  boundaries  and  phase  boundaries. 
Furthermore,  a  grain  boundary  may  be  wetted  by  either 
a  liquid  or  solid  phase  with  a  different  composition  from 
the  bulk,  with  the  former  associated  with  the  liquid  metal 
embrittlement  phenomenon  [77,78]  and  evidence  of  the 
latter  (“solid-state  wetting”)  found  in  Al-Zn  alloys  [79], 
Wynblatt  and  Chatain  [80]  modeled  solid-state  wetting 


5  In  Cahn’s  model  for  a  demixed  liquid,  the  prewetting  layer  is  liquid¬ 
like.  However,  in  the  general  case,  the  prewetting  layer  does  not 
necessarily  have  to  be  liquid-like  (e.g.  for  prewetting  in  solids). 
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transitions  in  a  binary  alloy  with  a  miscibility  gap. 
Although  their  investigation  is  based  on  a  regular  solution 
lattice-gas  model,  it  is  essentially  equivalent  to  CPWT,  with 
all  the  energy  terms  in  Eq.  (7)  represented  by  correspond¬ 
ing  lattice  interactions.  The  calculated  wetting  diagram 
(see  Fig.  7)  is  qualitatively  identical  to  Fig.  6.  Fike  grain 
boundary  energy,  the  wetting  temperature  at  a  grain 
boundary  is  found  to  be  anisotropic  and  inversely  related 
to  grain  boundary  energy.  Flnlike  wetting  by  a  liquid 
phase,  solid-state  wetting  transitions  may  incur  additional 
strain  energy  due  to  the  lattice  mismatch  between  solid 
phases  with  different  compositions,  which  is  not  included 
in  the  calculations  of  Wynblatt  and  Chatain.  Cahn  dis¬ 
cussed  the  implications  of  such  a  difference  on  wetting  tran¬ 
sitions  in  a  recent  paper  [81],  He  pointed  out  that,  when  the 
two  solid  phases  near  the  critical  point  of  a  miscibility  gap 
form  an  incoherent  interface  or  semi-coherent  interface  (i.e. 
lattice  mismatch  accommodated  by  interface  dislocations), 
wetting  transitions  will  not  occur.  A  wetting  transition  is 
possible  in  the  case  of  a  coherent  interface.  However, 
because  the  misfit  strain  energy  scales  with  the  thickness 
of  the  wetting  layer,  the  wetting  film  can  only  have  a  finite 
thickness  and  complete  wetting  is  thus  forbidden. 

Surface  or  interface  premelting  offers  another  example 
of  complexion  transitions.  Melting  of  a  solid  usually  starts 
from  its  surface  or  interface.  Free  surfaces  have  been 
observed  to  start  melting  at  a  temperature  Tsm  below  the 
bulk  melting  point  Tm  in  various  materials,  including  ice 
[70,71]  and  metals  [82,83],  At  Tsm  <  T  <  Tm,  the  equilib¬ 
rium  surface  configuration  switches  from  a  “dry”  surface 
state  to  a  liquid-like  surficial  film.  Such  a  film  represents 
a  specific  complexion,  and  has  an  equilibrium  thickness 
and  a  structure  less  ordered  and  less  dense  than  the  crystal 
phase.  When  raising  temperature  towards  Tm,  the  structure 
of  the  surficial  film  approaches  that  of  bulk  liquid  phase 
and  its  thickness  diverges  upon  reaching  Tm.  The  thermo¬ 
dynamic  driving  force  underlying  premelting  is  also  related 
to  minimization  of  surface  energy.  The  liquid  phase  of  a 
material  usually  has  a  lower  surface  energy  (yfc)  than  its 
crystalline  counterpart  (yW),  where  the  superscript  (0) 


Fig.  7.  (a)  Prewetting  line  (solid  blue)  superimposed  on  the  bulk  phase 
diagram  with  a  miscibility  gap.  Arrows  indicate  three  paths  of  increasing 
composition,  which  correspond  to  (1)  no  prewetting  transition,  (2)  first- 
order  prewetting  transition  and  (3)  continuous  prewetting  transition  that 
occurs  above  the  prewetting  critical  temperature,  TPWc .  (b)  Adsorption 
vs.  bulk  composition  diagram  for  three  paths  shown  in  (a);  note  the  first- 
order  discontinuous  jump  in  adsorption  for  path  2.  Adapted  with  from 
Ref.  [80],  (For  interpretation  of  the  references  to  color  in  this  figure 
legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


denotes  that  it  always  represents  a  pristine  “dry”  surface 
without  the  premelting  layer  while  the  true  ysv  represents 
a  premelted  surface  in  cases  where  premelting  does  occur, 
and  the  inequality  y^  >  ylB  +  yst  often  holds  for  high 
energy  surface  orientations  of  the  crystal.  Therefore,  it 
may  become  energetically  favorable  to  form  a  liquid-like 
surficial  film  on  top  of  the  crystal  phase  provided  that 
the  molar  free  energy  of  the  liquid  phase  is  not  much 
greater  than  that  of  the  stable  crystalline  phase,  i.e.  near 
Tm  (although  we  should  recognize  that  the  premelting  layer 
should  exhibit  partial  order  and  contain  gradient  energies 
that  make  its  energy  substantially  different  from  that  of  a 
bulk  liquid).  Similarly,  grain  boundary  premelting  may 
also  be  expected  when  the  grain  boundary  energy  is  larger 
than  twice  the  liquid/solid  interface  energy. 

Surface  (or  interface)  melting  can  be  viewed  as  a  special 
case  of  a  wetting  transition  in  which  the  interface  between  a 
crystal  and  a  third  phase  is  wetted  by  its  own  liquid  phase 
of  arbitrary  thickness  at  the  bulk  melting  point.  Similar  to 
the  way  that  wetting  and  prewetting  transitions  are  related, 
premelting  transitions  are  the  precursor  to  complete  melt¬ 
ing,  and  can  be  analyzed  and  predicted  by  CPWT-like 
mean-field  models.  Here  we  focus  on  models  for  grain 
boundary  prenrelting  transitions;  analogous  models  for 
surface  premelting  have  been  reviewed  in  Refs.  [78]  and 
[23],  Fobkovsky  and  Warren  [84]  and  Tang  et  al.  [13,85] 
studied  grain  boundary  prenrelting  using  a  diffuse-interface 
model  for  polycrystalline  materials  developed  by  Kobay- 
ashi,  Warren  and  Carter  (KWC)  [86,87],  In  the  KWC 
model  for  a  single-component  system,  a  two-dimensional 
polycrystalline  structure  is  described  by  two  field  variables, 
i.e.  a  crystallinity  field  rj(x)  and  a  crystallographic  orienta¬ 
tion  field  9(x).  i]  is  a  coarse-grained  measure  of  local  struc¬ 
tural  order  [88,89];  rj  =  1  and  0  are  defined  as  the  perfect 
crystalline  or  disordered  states,  respectively.  6  is  a  local 
measure  of  the  “most  likely”  crystallographic  orientation 
within  fixed  coordinates.  The  excess  free  energy  of  a  planar 
grain  boundary  located  at  x  =  0  is  given  by 


FGB[n,e-,T\  = 


v"  /  dtj 


4/("'r)+2U 


g{n) 


dx 


(8) 


which  comprises  the  homogeneous  free  energy  density  A f 
and  energy  terms  related  to  the  gradients  of  r]  and  0.  The 
orientation  gradient  energy,  g(tj)  |^|,  characterizes  the  en¬ 
ergy  penalty  associated  with  creating  a  misorientation 
across  the  grain  boundary.  Its  prefactor,  g,  is  required  to 
be  a  monotonically  increasing  function  of  ij  and  vanishes 
at  r]  =  0.  Such  a  requirement  originates  from  a  simple  phys¬ 
ical  intuition:  there  should  be  no  penalty  for  generating  an 
orientation  gradient  in  a  completely  disordered  phase  that 
has  a  randomly  orientated  structure  and  the  penalty  should 
increase  with  the  increasing  structural  order  of  the  mate¬ 
rial.  It  is  shown  that  an  additional  quadratic  term  of  the 
orientation  gradient  should  be  included  in  the  free  energy 
functional  to  model  curvature-driven  grain  boundary 
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motion  [90],  but  may  be  omitted  for  planar,  stationary 
boundaries. 

Kobayashi  and  Giga  [91]  proved  that,  for  an  equilib¬ 
rium  bicrystal  which  minimizes  Eq.  (8),  the  corresponding 
0  field  must  be  a  step  function  that  localizes  all  of  its 
change  at  the  grain  boundary  core  x  =  0  where  ;/  has  a  min¬ 
imum  t]GB.  Fig.  8(a)  shows  the  schematic  profiles  of  tj  and  0 
for  a  bicrystal  that  are  admitted  by  the  KWC  model.  With 
this  knowledge,  Tang,  Carter  and  Cannon  (TCC)  [13,85] 
wrote  Eq.  (8)  in  another  form: 


FGB[t];  AO ,  T]  =  g(r]GB)A0  +  2 


X 


A /(f?,  T) 


(9) 


where  AO  is  the  difference  between  the  orientations  of  the 
two  neighboring  grains,  i.e.  the  misorientation,  and  a  sym¬ 
metrical  i]  profile  about  the  grain  boundary  is  assumed. 
Notably,  Eq.  (9)  has  the  same  form  as  the  free  energy 
expression  in  CPWT.  The  surface  energy  term  <P{  cs)  in 
Eq.  (7)  favors  a  surface  composition  that  is  different  from 
the  bulk  liquid  composition  preferred  by  the  homogeneous 
free  energy  /.  Similarly,  the  misorientation  energy  penalty 
g(y]GB)A0  in  Eq.  (9)  favors  a  grain  boundary  more  disor¬ 
dered  =  0)  than  the  bulk  crystalline  phase  stabilized 
by  A/.  Inspired  by  such  an  analogy,  a  graphic  construction 
method  similar  to  Cahn’s  approach  is  employed  in  Ref.  [13] 
to  analyze  the  equilibrium  grain  boundary  structure  near 
Tm.  As  summarized  in  Fig.  8(b),  the  analysis  predicts  three 
types  of  grain  boundary  behavior  near  Tm  that  depend  on 
the  grain  boundary  misorientation  value. 


(i)  For  small  misorientations  (AO  <  A0wet),  the  grain 
boundary  energy  is  smaller  than  twice  the  liquid/solid 
interface  energy,  yGB  <  2ysh  The  boundary  retains  a 
relatively  ordered  structure  up  to  the  bulk  melting 
point  and  can  be  superheated  above  Tm. 


(ii)  A  first-order  premelting  transition  can  occur 
below  Tm  at  intermediate  misorientation  values 
(A 0wet  <  A 0  <  A0crit),  which  produces  a  discontinuous 
increase  in  the  width  and  structural  disorder  of  the 
grain  boundary.  The  premelting  line  on  the  T  AO  dia¬ 
gram  starts  at  ( Tm ,  A 0wet)  and  terminates  at  a  grain 
boundary  critical  point  (Tcrit  A 0crit). 

(iii)  At  large  misorientations  (AO  >  A 0cril),  a  grain  bound¬ 
ary  premelts  continuously,  i.e.  the  grain  boundary 
crystallinity  rjGB  decreases  to  0  and  the  grain  bound¬ 
ary  width  diverges  gradually  without  a  discontinuity 
when  approaching  Tm. 

Because  the  KWC  model  is  isomorphic  to  CPWT,  it  is 
not  surprising  to  see  that  the  different  types  of  predicted 
grain  boundary  premelting  behavior  have  a  parallel  rela¬ 
tion  with  Cahn’s  predictions  of  the  prewetting  phenomena. 
The  grain  boundary  premelting  diagram,  Fig.  8(b),  looks 
much  like  the  prewetting  diagram  in  the  /i  T  space 
(Fig.  6(b)). 

Although  grain  boundary  premelting  in  pure  materials 
has  been  predicted  by  theory  [13]  and  atomistic  simulations 
(see  Section  2.4)  and  suggested  by  some  indirect  experimen¬ 
tal  evidence,  such  as  abrupt  changes  in  grain  boundary  dif- 
fusivity  and  sliding  [55,92,93],  its  direct  observation  (e.g.  by 
HRTEM)  in  real  material  systems  has  been  elusive  [94],  In 
contrast,  evidence  of  grain  boundary  premelting  in  non- 
pure  systems  is  much  more  abundant.  For  example,  nano¬ 
meter-thick  disordered  interfacial  films,  known  as  “inter¬ 
granular  films”  (IGFs),  have  been  widely  observed  at 
grain  boundaries  in  various  ceramics  that  contain  inten¬ 
tionally  (e.g.  as  sintering  aids)  or  unintentionally  added 
impurities  such  as  Si3N4-Si02  (impurity),  ZnO  Bi203, 
SrTi03-Ti02,  etc.  Similar  nanoscale  films  have  also  been 
confirmed  or  deduced  in  metallic  alloys  (e.g.  W-Ni, 
Mo-Ni,  Cu-Bi,  etc.).  Comprehensive  reviews  of  the  IGF 
phenomenon  can  be  found  in  Refs.  [63,95].  Numerous 


Fig.  8.  (a)  Schematic  diagram  of  equilibrium  profiles  of  )/  (degree  of  crystallinity)  and  6  (orientation)  fields  for  a  grain  boundary  at  x  =  0.  AO  is  the 
misorientation  between  the  two  grains;  its  magnitude  represents  the  grain  boundary  energy.  The  degree  of  crystallinity,  r/,  reaches  a  minimum  at  the  center 
of  the  grain  boundary  core  (x  =  0).  (b)  A  schematic  grain  boundary  premelting  diagram  for  a  pure  system  in  temperature-misorientation  ( T-AO)  space. 
The  premelting  line  terminates  at  a  critical  point  marked  by  the  red  dot.  Grain  boundaries  with  large  misorientation  (i.e.  high  energy)  above  A 0cril 
continuously  premelt  as  the  bulk  melting  temperature  T„,  is  approached,  while  grain  boundaries  with  intermediate  AO  (energy)  exhibit  first-order 
premelting  as  the  temperature  is  raised.  Small  AO  (low  energy)  grain  boundaries  with  misorientations  smaller  than  the  A 0we,  do  not  undergo  premelting. 
Adapted  front  Ref.  [13], 
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experiments  have  confirmed  that  IGFs  are  in  thermody¬ 
namic  equilibrium  with  the  bulk  phase.  They  exhibit  a 
liquid-like  structure  and  equilibrium  thickness,  and  exist 
under  conditions  off  solid-liquid  coexistence,  which  agree 
with  the  characteristics  of  premelting  transitions.  Com¬ 
pared  to  the  abutting  grains,  IGFs  are  enriched  in  impuri¬ 
ties;  their  formation  causes  the  original  grain  boundary  to 
develop  multilayer  solute  adsorption,  which  is  reminiscent 
of  prewetting.  IGFs  can  hence  be  viewed  as  products  of 
coupled  premelting/prewetting  transitions. 

To  provide  a  thermodynamic  description  of  premelting 
in  non-pure  materials,  the  TCC  analysis  for  pure  systems 
[13]  was  extended  to  binary  systems  [16],  where  the  excess 
free  energy  of  a  grain  boundary  is  given  as 


FGB[,hc-A9,T] 
=  g(r 1gb)ao 


+  2 


A/(c,  1,  T) 


v2  / drj\ 
2  \dx ) 


k2  ( dc\ 
2  \dx) 


dx 


(10) 


Fig.  9.  A  schematic  grain  boundary  premelting  diagram  for  an  A-B 
binary  system  which  has  a  eutectic  transition  at  Te.  The  prewetting  lines 
(solid  red)  of  three  grain  boundaries  with  different  misorientation  AO 
(AOS<A0,„  <  AO,)  are  shown  in  the  oc  single-phase  region.  Also  shown  are 
the  complete  wetting  ( Twer )  or  premelting/prewetting  (Tpw)  temperatures 
(solid  blue  lines)  for  these  boundaries  in  the  a  +  liquid  or  a  +  ft  two-phase 
regions.  (For  interpretation  of  the  references  to  color  in  this  figure  legend, 
the  reader  is  referred  to  the  web  version  of  this  article.) 


Compared  to  the  free  energy  of  pure  systems  (Eq.  (9)), 
A f  is  now  also  composition-dependent  and  an  additional 
composition  gradient  energy  term  appears  in  Eq.  (10).  An 
analogous,  but  more  involved,  graphic  construction  was 
developed  to  determine  the  equilibrium  grain  boundary 
state  that  minimizes  Eq.  (10)  [16],  While  the  binary  model 
also  predicts  three  categories  of  melting  behavior  (super¬ 
heating,  first-order  and  continuous  premelting)  in  alloys, 
conditions  for  premelting  display  more  intricate  depen¬ 
dence  on  the  thermodynamic  state  variables. 

The  TCC  analysis  of  grain  boundary  premelting  is  based 
on  simple  model  systems  in  which  the  grain  boundary 
energy  increases  monotonically  with  misorientation.  In  real 
materials,  the  misorientation  dependence  of  grain  bound¬ 
ary  energy  is  usually  more  complex,  e.g.  there  exist  low- 
energy  cusps  at  large  misorientations  that  correspond  to 
special  grain  boundary  geometries  and  low  index  interface 
planes.  It  is  thus  proper  to  view  Ad  as  a  measure  of  grain 
boundary  energy  when  interpreting  the  predicted  grain 
boundary  premelting  behavior  as  a  function  of  grain 
boundary  misorientation.  Mishin  et  al.  [61]  recently  exam¬ 
ined  premelting  transitions  in  Cu-Ag  in  a  multi-phase-field 
model,  using  realistic  bulk  thermodynamic  description  of 
Cu-Ag  alloys  from  the  CALPHAD  approach.  Their  calcu¬ 
lations  reveal  several  classes  of  melting  behavior  at  bound¬ 
aries  of  different  energies  that  agree  with  the  generic 
predictions  of  Ref.  [16], 

Fig.  9  illustrates  a  schematic  grain  boundary  premelting 
diagram  predicted  by  the  TCC  model  for  a  eutectic  alloy 
[13],  to  which  many  IGF-forming  systems  belong.  It  shows 
that,  for  states  in  the  solid-liquid  two-phase  coexistence 
region,  a  grain  boundary  can  be  perfectly  wetted  by  the 
equilibrium  liquid  phase  with  a  composition  on  the  liqui- 
dus  above  a  complete  wetting  temperature  Twet.  Similar 
to  prewetting  in  binary  liquid,  a  premelting/prewetting  line 


extends  into  the  single-phase  region  from  Twet.  Upon  cross¬ 
ing  this  line  from  low  to  high  solute  (B)  saturation  (or  tem¬ 
perature),  the  grain  boundary  transforms  to  a  more 
disordered  structure  with  a  higher  solute  segregation  via 
a  first-order  transition.  Importantly,  there  exists  a  distinct 
premelting/prewetting  line  for  each  grain  boundary  with 
a  different  misorientation  (i.e.  different  grain  boundary 
energy),  and  Twet  increases  with  decreasing  misorientation 
AO  (i.e.  with  decreasing  grain  boundary  energy).  When  a 
boundary’s  Twe,  is  below  the  eutectic  temperature  Te,  the 
premelting/prewetting  line  intersects  the  solubility  line  of 
the  solid-solid  two-phase  region  at  Tpw.  In  this  case,  com¬ 
plete  wetting  does  not  occur  to  the  grain  boundary  in  the 
two-phase  region  below  Te  as  the  liquid  phase  is  metastable 
there.  Instead,  the  grain  boundary  will  undergo  a  first- 
order  premelting/prewetting  transition  in  the  two-phase 
region  at  Tpw,  as  shown  in  Fig.  9.  The  premelted  intergran¬ 
ular  film  will  thicken  and  completely  wet  the  grain  bound¬ 
ary  when  T  reaches  Te,  at  which  a  stable  liquid  phase  exists. 

The  predictions  of  the  diffuse-interface  model  are  in 
qualitative  agreement  with  the  experimental  observations 
of  IGFs.  The  model  prediction  that  thicker,  disordered 
grain  boundaries  can  form  in  subeutectic  and  single-phase 
regions  is  consistent  with  the  experimental  conditions 
under  which  IGFs  were  found.  The  coupled  premelting/ 
prewetting  transition  leads  not  only  to  structural  disorder¬ 
ing  but  also  to  a  higher  segregation  level  of  the  minority 
species  at  the  grain  boundary.  Nevertheless,  the  diffuse- 
interface  model  suggests  that  the  equilibrium  compositions 
of  the  premelted  boundary  differ  from  either  solid  or  liquid 
bulk  phase,  which  also  agrees  with  experiments  [96-100], 
As  shown  in  Fig.  9,  a  grain  boundary’s  premelting/prewett¬ 
ing  line  may  terminate  at  a  grain  boundary  critical  point 
and  not  extend  to  the  end  member  A,  which  means  that 
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the  grain  boundary  can  undergo  a  first-order  premelting 
transition  in  alloys  but  not  in  the  pure  system.  This  demon¬ 
strates  that  the  conditions  for  a  first-order  premelting  tran¬ 
sition  could  be  less  stringent  in  non-pure  systems  than  in 
pure  materials. 

One  limitation  of  Cahn’s  CPWT  and  similar  diffuse- 
interface  models  for  prewetting  and  premelting  transitions 
is  that  only  short-range  interactions  are  considered  in  the 
free  energy  of  the  systems.  On  the  other  hand,  long-range 
interactions,  such  as  London  dispersion  and  electrostatic 
forces,  are  known  to  influence  the  development  and  charac¬ 
teristics  of  interface  wetting  or  melting  transitions  in  signif¬ 
icant  ways  [70,71,76],  especially  for  inorganic  materials. 
For  example,  the  presence  of  an  attractive  dispersion  force 
at  a  grain  boundary  could  prevent  it  from  being  completely 
wetted  by  the  liquid  phase  where  a  liquid  and  a  solid  coex¬ 
ist  [101],  Bishop,  Carter  and  Cannon  developed  a  diffuse- 
interface  model  formulation  that  incorporates  the  electro¬ 
static  interaction  between  charged  ion  species  in  the  system 
[15,102],  The  model  was  used  to  numerically  study  the  for¬ 
mation  and  stability  of  IGFs  in  Si3N4-Si02  [102],  Ebner 
and  Saam  [103]  considered  the  effects  of  long-range  disper¬ 
sion  force  on  wetting  phenomena  in  a  mean-field  theory  of 
an  Ising  lattice-gas  system.  They  found  that  the  inclusion  of 
the  long-range  interaction  significantly  complicates  the 
wetting  behavior  of  an  adsorbate  on  a  substrate  near  the 
bulk  critical  point.  Alternatively,  the  effects  of  long-range 
forces  on  IGFs  can  conveniently  be  taken  into  account  in 
so-called  force-balance  models,  which  will  be  discussed  in 
Section  2.4. 

2.3.3.  Discrete  complexions 

Grain  boundary  complexions  often  occur  in  discrete 
thicknesses  and  can  be  classified  by  their  width  in  terms 
of  the  number  of  atomic  layers.  This  complexion  categori¬ 
zation  scheme  has  its  origin  in  a  grain  growth  kinetics  study 
on  doped  alumina,  in  which  high-angle  annular  dark  field 
scanning  transmission  electron  microscopy  (HAADF- 
STEM)  imaging  was  used  to  observe  the  grain  boundary 
complexions  of  specific  grain  boundaries  with  different 
mobility.  As  shown  in  Fig.  10  (a-f),  the  discrete  complex¬ 
ions  fall  into  six  groups  based  on  their  thickness,  as 
observed  in  either  HAADF-STEM  or  HRTEM  imaging. 
The  physical  origin  of  these  different  complexion  thick¬ 
nesses  will  vary  from  one  materials  system  to  the  next; 
the  categorization  scheme  is  simply  based  on  the  approxi¬ 
mate  thickness.  These  six  discrete  complexions  are  referred 
to  as  Dillon-Harmer  complexions  [25],  and  include  clean 
complexions  (with  no  observable  segregation  of  solute 
atoms),  monolayer  complexions,  bilayer  complexions,  tri¬ 
layer  complexions,  nanolayer  complexions  and  wetting  lay¬ 
ers.  These  six  complexion  types  were  first  discovered  in  Ref. 
[46]  and  formally  categorized  as  complexions  in  Ref.  [17], 
Within  the  Dillon-Harmer  complexion  categorization 
framework,  there  is  a  general  trend  of  increasing  adsorp¬ 
tion  and  structural  and  chemical  disorder  with  increasing 
thickness.  In  practice,  the  Dillon-Harmer  complexions 


are  defined  and  categorized  based  on  STEM  observations 
of  thickness  and  whether  or  not  adsorbed  solute  is  visible 
at  the  grain  boundary. 

We  emphasize  that  the  phrase  “monolayer  complexion” 
should  not  be  confused  with  the  surface  science  term 
“monolayer”  (ML);  the  former  phrase  is  a  type  of  complex¬ 
ion  within  the  Dillon-Harmer  scheme,  while  the  latter 
word  is  a  quantity  of  adsorbate  that  can  be  defined  in  sev¬ 
eral  different  ways  [104].  To  avoid  confusion,  when  discuss¬ 
ing  Dillon-Harmer  complexions  it  is  best  to  include  the 
word  complexion.  In  this  work,  we  will  use  the  abbrevia¬ 
tion  ML  when  referring  to  the  concept  of  a  specific  quan¬ 
tity  of  adsorbate. 

Nanolayer  complexions  are  greater  than  three  atomic 
layers  in  thickness  but  are  not  a  bulk  phase;  these  complex¬ 
ions  are  typically  on  the  order  of  one  or  two  nanometers 
thick.  Historically,  nanolayer  complexions  have  been 
referred  to  as  IGFs.  The  last  category,  wetting,  is  actually 
a  bulk  phase  rather  than  a  complexion.  It  has  an  arbitrary 
thickness  while  the  other  categories  all  have  a  fixed,  finite 
thickness  set  by  equilibrium  thermodynamics.  While  the 
six  discrete  complexions  are  a  useful  method  for  categoriz¬ 
ing  complexions,  it  should  not  be  inferred  that  all  nanolay¬ 
ers,  for  example,  have  the  same  thickness,  or  that  there  can 
only  be  a  single  type  of  nanolayer  complexion  within  a 
given  system:  there  may  be  multiple  types  of  nanolayers 
with  different  thicknesses  within  the  same  system. 

The  six  discrete  Dillon-Harmer  complexions  exhibit 
similarities  to  a  simplified  model  where  hard  spheres  are 
sandwiched  between  two  hard  walls  (Fig.  11(b)).  In  such 
a  case,  well-known  colloidal  theories  suggest  that  an  oscil¬ 
latory  structural  interaction  will  produce  free  energy  min¬ 
ima  that  correspond  to  discrete  thicknesses  (at 
D  =  na  +  1/4 a  for  the  case  of  hard  spheres  sandwiched 
between  two  hard  walls,  where  D  is  the  distance  between 
two  walls,  a  is  the  diameter  of  the  hard  sphere  and  n  is  a 
positive  integer  (n=  1,2,3,...);  Fig.  11(c))  [105],  It  was 
proposed,  based  on  a  phenomenological  model  [22,106], 
that  a  similar  mechanism  may  produce  a  series  of  discrete 
Dillon-Harmer  complexions  (Fig.  11(a)).  Since  the  statisti¬ 
cal  model  suggests  that  the  (n  +  l)th  energy  barrier  should 
decay  exponentially  by  a  factor  of  l/<?  as  compared  with  the 
nth  energy  barrier,  the  thickness  of  a  complexion  should  be 
nearly  continuous  beyond  the  trilayer  complexion  (how¬ 
ever,  rigorously,  it  can  still  be  discrete  with  smaller  barri¬ 
ers).  This  will  result  in  approximately  six  complexions,  as 
shown  in  Fig.  1 1(a).  It  is  also  possible  to  have  complexions 
with  discrete  thicknesses  equivalent  to  four,  five  or  more 
atomic  layers;  however,  because  the  energy  barriers 
between  these  thicker  discrete  complexions  are  small,  they 
are  considered  to  all  be  nanolayer  complexions.  Beyond  the 
“roughening  point”,6  the  effective  thickness  of  the  complex¬ 
ion  will  be  continuous,  where  nanolayers  with  continuous 


6  In  this  context,  “roughening”  refers  to  the  changeover  in  behavior 
from  discrete  changes  in  thickness  to  continuous  changes  in  thickness. 
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Fig.  10.  The  six  discrete  Dillon-Harmer  complexions  as  originally  discovered  in  undoped  and  doped  (CaO,  MgO,  Si02,  Nd203)  A1203  (a-f)  and 
analogous  examples  of  discrete  Dillon-Harmer  complexions  in  metallic  systems  (g  1).  The  relationship  between  dopant  type  and  grain  boundary 
complexion  type  for  the  A1203  complexions  can  be  seen  in  Fig.  21.  Parts  (a-f)  reprinted  from  Ref.  [17]  with  permission;  parts  (g,  h)  reprinted  from  Ref. 
[163]  with  permission;  part  (i)  reprinted  from  Ref.  [48]  with  permission;  part  (j)  reprinted  from  Ref.  [344]  with  permission;  parts  (k,  1)  reprinted  from  Ref. 
[205]  with  permission. 


but  finite  “equilibrium”  thickness  are  thermodynamically 
stable  (Fig.  12).  On  the  other  hand,  discrete  complexions 
with  a  thickness  of  approximately  zero,  one,  two  or  three 
atomic  layers  appear  below  the  roughening  points  where 
the  so-called  stepwise  “layering  transitions”  occur. 

A  seminal  article  by  Pandit,  Schick  and  Wortis  provided 
further  insight  regarding  the  origin  of  discrete  complexions 
[107],  Via  a  lattice-gas  model  for  multilayer  adsorption  of 
inert  gas  molecules  on  an  attractive  substrate,  Pandit 
et  al.  showed  that  layering  transitions  can  occur  below 
roughening  transition  points,  producing  discrete  interfacial 
states  with  thicknesses  of  approximately  a  single  atomic 
layer,  two  atomic  layers,  three  atomic  layers,  etc. 
(Fig.  12(a)).  A  recently  proposed  model  [25]  further  sug¬ 
gested  that  similar  layering  transitions  can  occur  at  grain 
boundaries  to  produce  discrete  Dillon-Harmer  complex¬ 
ions  (Fig.  12(b)).  Semi-grand  canonical  Monte  Carlo  simu¬ 
lations  have  also  recently  demonstrated  a  transition 
between  two  different  discrete  Dillon-FIarmer  complexions 
(monolayer  and  bilayer  complexions)  at  a  £5(310)  grain 
boundary  in  Ag-doped  Cu,  as  shown  in  Fig.  4(d  and  e). 

In  reality,  the  adsorbates  at  grain  boundaries  are  not 
“hard  spheres”.  Thus,  the  ideal  situation  predicted  by  the 
simplified  models  shown  in  Figs.  11  and  12  is  not  typically 
observed  in  experiments.  For  example,  a  recent  study  [34] 
revealed  a  first-order  grain  boundary  transition  from  a 
bilayer  complexion  to  a  clean  complexion  in  Si-Au  in  the 
absence  of  the  intermediate  state  of  a  monolayer  complex¬ 
ion.  Specifically,  Fig.  13(a)  shows  the  coexistence  of  a 
bilayer  and  a  clean  complexion  at  a  low-symmetry  (£43) 
twist  grain  boundary;  the  atomically  abrupt  transitions 
between  the  two  complexions  suggests  that  this  complexion 
transition  is  likely  to  be  a  first-order  transition,  because  it  is 
associated  with  a  discontinuity  in  the  interfacial  excess  of 
the  solute  (Au)  adsorption. 


Fig.  13(b)  schematically  illustrates  possible  relative  sta¬ 
bilities  of  a  bilayer,  a  monolayer  and  a  clean  complexion. 
It  was  proposed  that  a  bilayer  complexion  in  Si-Au  is  sta¬ 
bilized  because  Au  atoms  bond  strongly  to  the  Si  atoms  on 
the  adjacent  Si  grain  surface  (since  the  regular  solution 
parameter  for  Au-Si  amorphous  alloy  was  estimated  to 
be  Q  ~  —  8  kJ  mol-1,  Au-Si  bonds  are  energetically 
favored  to  form)  [34],  similar  to  the  case  of  Ni-Bi  [48], 
As  illustrated  in  Fig.  13(b),  it  may  be  energetically  more 
expensive  to  form  an  Au  monolayer  complexion  at  this 
low-symmetry  boundary  because  the  Au  monolayer  com¬ 
plexion  cannot  grow  coherently  with  respect  to  both  grain 
surfaces,  i.e.  some  strong  Au-Si  bonds  must  be  broken  at 
(at  least)  one  Si-Au  interface.  In  other  words,  the  mono- 
layer  complexion  may  represent  a  high  energy  state  because 
both  Si  grain  surfaces  would  tend  to  impose  structural 
order  onto  the  monolayer  complexion  and  this  structural 
ordering  is  incompatible,  leading  to  a  structural  frustration 
that  destabilizes  the  monolayer  complexion  [34], 

In  summary,  the  grain  boundary  complexion  forma¬ 
tion  and  transitions  can  be  complex  and  material  spe¬ 
cific.  For  ceramic  materials,  directional  bonds,  the 
presence  of  both  cations  and  anions,  electrostatic  space 
charges  and  van  der  Waals  forces  can  make  the  interfa¬ 
cial  interactions  much  more  complicated.  For  example, 
an  analogous  surface  complexion  transition  from  the 
clean  complexion  to  the  nanolayer  complexion  (i.e.  the 
free-surface  counterpart  to  nanolayer  complexions  at 
grain  boundaries  [23])  was  observed  in  V205  on  Ti02 
[24,108],  In  many  ceramic  systems,  complete  wetting 
may  be  inhibited  by  the  presence  of  long-range  attrac¬ 
tive  van  der  Waals  interactions  [63,109],  More  realistic 
models  that  can  predict  the  specific  complexion  stability 
and  transitions  in  any  (or  a  broad  range  of)  real  mate¬ 
rials  are  yet  to  be  developed. 
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(a) 


Clean  Monolayer 


general  trend  of  increasing  adsorption  and 
structural  &  chemical  disorder 


Fig.  11.  The  six  discrete  Dillon-Harmer  complexions  can  be  viewed  as  somewhat  analogous  to  a  simplified  model  where  hard  spheres  are  sandwiched 
between  two  hard  walls,  (a)  Schematic  depiction  of  the  six  discrete  Dillon-Harmer  complexions;  (b)  hard  spheres  of  diameter  a  sandwiched  between  hard 
walls  of  spacing  D\  (c)  pressure  (red  line)  and  energy  (blue  line)  as  a  function  of  wall  spacing  D.  As  D  increases,  there  are  periodic  pressure  and  energy 
fluctuations.  The  equilibrium  hard  wall  spacing  D  occurs  at  values  of  minimum  energy,  where  D  =  na  +  'Aa  (open  circles  on  blue  line).  Parts  (b,  c)  adapted 
from  Ref.  [105],  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


2.4.  Other  theories,  models  and  important  concepts 

Historically,  two  types  of  models  have  been  employed  to 
describe  wetting  transitions:  diffuse-interface  and  force-bal¬ 
ance  models.  In  Section  2.3,  we  discussed  diffuse-interface 
models  for  interface  prewetting  and  premelting  complexion 
transitions,  in  which  interfaces  are  considered  to  be  intrin¬ 
sically  diffuse  with  a  characteristic  finite  width.  In  contrast, 
a  wetting  (or  melting)  configuration  in  force-balance  mod¬ 
els  is  phenomenologically  treated  as  a  wetting  (or  melting) 
layer,  i.e.  a  bulk  phase,  sandwiched  by  two  sharp  interfaces 
of  zero  thickness.  The  equilibrium  thickness  h  of  the  wet¬ 
ting  (or  melting)  film,  and  hence  the  wetting  state  (h  =  0, 
h  =  oo,  0  <  h  <  oo),  is  determined  by  the  balance  of  various 


forces  between  the  two  interfaces  that  are  thickness- 
dependent. 

A  force-balance  model  was  introduced  by  Clarke  et  al. 
[110,111]  to  explain  the  presence  of  stable  IGFs  at  grain 
boundaries.  A  fundamental  assumption  of  this  model  is  that 
IGFs  are  quenched  liquid  films  of  the  secondary  phase  con¬ 
fined  by  two  grains  of  the  primary  phase.  In  the  first  paper 
[110],  Clarke  assumes  that  IGFs  are  mainly  stabilized  by  an 
attractive  long-range  dispersion  force  ndisp  vs.  a  repulsive 
short-range  steric  force  nsteric  with  the  following  expressions: 


nsteric  =  -4 cj)Q  exp (~h/£)  (12) 
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Fig.  12.  (a)  A  schematic  surface  complexion  diagram  (top)  and  the  corresponding  surface  excess  adsorption  vs.  chemical  potential  curves  (bottom)  for 
multilayer  surface  adsorption  of  noble  gas  molecules  on  an  attractive  inert  substrate,  (b)  A  computed  grain  boundary  complexion  diagram  (top)  and  the 
corresponding  grain  boundary  excess  adsorption  vs.  normalized  bulk  composition  curves  (bottom).  Panel  (a)  is  replotted  from  Ref.  [107]  with  permission 
from  APS  and  panel  (b)  is  replotted  from  Ref.  [25]  with  permission  from  AIP. 


where  h  is  the  film  thickness,  A12 1  is  the  Hamaker  constant 
of  the  grain/IGF/grain  composite  structure,  £  is  the  corre¬ 
lation  length  of  the  liquid  phase  and  4>0  is  a  strength  con¬ 
stant.  The  repulsive  steric  force  is  also  known  as  the 
structural  disjoining  pressure.  The  origin  of  such  a  force 
is  related  to  the  phenomenon  that  a  liquid  layer  of  micro¬ 
scopic  thickness  confined  between  two  solid  interfaces 
tends  to  have  a  more  ordered  structure  and  thus  a  higher 
free  energy  density  than  the  bulk  liquid;  its  structural  order 
decreases  with  increasing  interface  separation,  which  is 
equivalent  to  an  effective  repulsive  interaction  between 
the  two  opposing  interfaces.  Interestingly,  the  expression 
of  IIsteric  was  derived  from  a  diffuse-interface  model  [110], 
Under  Clarke’s  simplifying  assumption  that  steric  and 
London  dispersion  forces  dominate  the  equilibrium  IGF 
thickness,  the  equilibrium  thickness  (which  occurs  at  zero 
net  force,  F)  can  be  found  by  setting  Eqs.  (11)  and  (12) 
equal  to  each  other  and  solving  for  h.  Using  parameters 
estimated  for  the  Al203-Si02  system,  Clarke  shows  that 
the  net  force,  F  =  ndisp  +  nsteric,  vanishes  at  an  IGF  thick¬ 
ness  of  ~1  nm,  which  is  in  qualitative  agreement  with 
experimental  observations.  The  possible  contribution  of 
the  electric  double-layer  interaction  to  force  balance  was 
also  considered  by  Clarke  et  al.  later  [111].  In  contrast  to 
IGFs  at  ceramic-ceramic  interfaces  in  which  dispersion 
forces  are  the  primary  attractive  force,  recent  research  on 
IGFs  at  metal-ceramic  interfaces  has  shown  that  electro¬ 
static  image  forces  are  the  dominant  attractive  force  at 
gold-sapphire  interfaces  [42], 


Since  the  Clarke  model  was  proposed  more  than  two 
decades  ago,  there  has  been  an  ever  increasing  amount  of 
evidence  suggesting  that  IGFs  are  in  thermodynamic  equi¬ 
librium  with  the  surrounding  bulk  phases  rather  than  being 
a  quenched  liquid.  There  are  also  cases  in  which  IGFs  rep¬ 
resent  a  metastable  equilibration,  e.g.  as  shown  for  Si02- 
based  IGFs  at  Ti02  grain  boundaries  by  Adder  and 
Chiang  [112].  This  suggests  that  IGFs  can  also  form  under 
sub-eutectic  or  sub-solidus  conditions,  where  bulk  liquid  is 
metastable;  they  are  associated  with  an  excess  volumetric 
free  energy  A fvo/h  relative  to  the  stable  bulk  solid  phase, 
which  translates  into  a  significant  attractive  force  77„o/  =  - 
A fvoi  on  the  films  and  should  be  accounted  for  in  the  calcu¬ 
lation.  Such  a  modification  has  been  applied  to  study  the 
stability  of  IGFs  and  surflcial  amorphous  films  [63,113], 
IGFs  of  an  equilibrium  thickness  can  also  form  above 
the  solidus  line,  in  equilibrium  with  a  partial-wetting 
(non-wetting)  bulk  liquid  phase.  In  such  a  case,  the  average 
composition  of  IGFs  can  differ  markedly  from  the  bulk 
liquid  phase,  despite  a  thermodynamic  equilibrium 
between  the  liquid-like  IGF  complexion  and  the  bulk  liquid 
(see  Refs.  [23,63]  and  references  therein). 

In  addition  to  continuum-level  modeling,  atomistic 
modeling  has  also  been  extensively  used  to  study  complex¬ 
ion  transitions,  which  provides  detailed  insights  into  the 
atomic  mechanisms  of  such  transitions.  There  exists  a  large 
body  of  atomistic  simulation  studies  of  surface  melting 
transitions  in  the  literature  (see  e.g.  the  review  in  Ref. 
[114]).  In  an  early  modeling  study  of  grain  boundary 


P.R.  Cantwell  et  al.  I  Acta  Materialia  62  (2014)  1-48 


17 


•  Si 

•  Au 


structural 

frustration 


\ 


coherent 

(Au-Si) 


\  1 


t 

incoherent 

(Si-Si ) 


ft 

Au-Si  coherency  cannot  be 
maintained  on  both  sides 


t 

incoherent 

(Au-Au) 


Fig.  13.  (a)  HAADF-STEM  image  of  bicrystal  boundary  in  Au-doped  Si  showing  the  abrupt  transition  region  between  the  bilayer  complexion  and  the 
clean  complexion.  The  abrupt  transition  indicates  the  occurrence  of  a  first-order  complexion  transition.  The  bicrystal  is  a  (1  1  1)  15°  twist  boundary,  i.e. 
E43.  ( b)  Schematic  illustration  of  the  possible  free  energy  states  that  may  lead  to  this  first-order  transition.  It  was  proposed  lhat  the  monolayer  complexion 
is  a  very  shallow  metastable  state  (and  hence  is  essentially  unstable)  because  it  is  energetically  favorable  to  form  strong  Au-Si  bonds,  (c)  Schematic 
diagrams  of  clean,  monolayer  and  bilayer  complexions  showing  that  the  Au-rich  monolayer  complexion  is  structurally  metastable.  Reprinted  from  Ref. 
[34]  with  permission. 


premelting,  Kikuchi  and  Cahn  [68]  analytically  solved  a 
lattice  gas  model  with  a  simple  interatomic  potential  and 
found  that  continuous  premelting  occurs  at  a  tilt  grain 
boundary  that  exhibits  a  logarithmic  divergence  in  thick¬ 
ness  when  approaching  the  bulk  melting  point.  A  similar 
conclusion  was  obtained  from  Monte  Carlo  (MC)  simula¬ 
tions  of  the  same  model  [115], 

Grain  boundary  melting  behavior  in  elemental  systems 
has  been  extensively  investigated  by  molecular  dynamics 
simulations  using  (semi-)empirical  interatomic  potentials 
[116-128],  where  grain  boundary  structural  transitions 
are  detected  by  measurement  of  changes  in  grain  boundary 
structural  order  parameter(s)  or  other  derived  properties, 
such  as  grain  boundary  diffusivity  [117,125]  and  sliding 
resistance  [126],  The  existence  and  nature  of  premelting 
transitions  are  not  unequivocally  resolved  by  these  studies, 
which  sometimes  yield  inconsistent  results.  A  number  of 
MD  simulations  show  that  high-angle  grain  boundaries 
tend  to  increase  continuously  in  disorder  with  increasing 


temperature,  indicating  continuous  premelting 
[116,118,128],  while  some  simulations  find  that  the  grain 
boundaries  still  retain  significant  crystallinity  at  tempera¬ 
tures  close  to  Tm  [121-123],  Other  simulations  suggest  that 
a  first-order  premelting  transition  occurs  below  Tm 
[119,120,124],  As  pointed  out  by  Tang  et  al.  [13],  such 
inconsistencies  might  originate  from  numerical  artifacts, 
failure  in  accurate  determination  of  bulk  melting  points, 
differences  in  atomic  potentials  used  and  challenges  in 
obtaining  ground  state  grain  boundary  structures 
[129,130].  Recently,  von  Alfthan  et  al.  [128]  examined  the 
structures  of  several  Si  twist  grain  boundaries  near  Tm  with 
improved  simulation  techniques.  Among  the  three  (001) 
twist  grain  boundaries  (£25,  £5  and  £29,  corresponding 
to  twist  angles  of  16°,  37°  and  44°)  they  studied,  evidence 
of  continuous  premelting  was  found  for  high-angle/high- 
energy  grain  boundaries  (£5  and  £29),  starting  at  0.7- 
0.8 Tm,  but  not  at  the  low-energy  £25  grain  boundary.  This 
finding  reaffirms  the  dependence  of  the  grain  boundary 
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melting  behavior  on  the  “dry”  grain  boundary  energy,  as 
suggested  by  phenomenological  models  [13,61],  Neverthe¬ 
less,  no  first-order  grain  boundary  transition  was  observed 
in  their  simulations. 

It  is  much  more  challenging  to  rigorously  determine  the 
existence  and  characteristics  of  premelting  transitions  in 
multicomponent  systems  from  MD  simulations  because 
of  the  often  impractically  long  simulation  time  required 
to  establish  equilibrated  composition  profiles  in  the  sys¬ 
tems.  Existing  MD  studies  of  IGFs  typically  prepare  inter¬ 
granular  films  through  simulated  high-temperature 
equilibration  and  quenching  processes.  While  the  struc¬ 
tures  obtained  may  not  necessarily  correspond  to  the  low¬ 
est  energy  states,  they  provide  valuable  insights  on  the 
detailed  bonding  structures  and  solute  distributions  within 
the  films  as  well  as  their  implications  for  various  properties 
[131-138],  For  example,  Garofalini  and  co-workers  [131  — 
136]  have  systematically  studied  IGFs  between  silicon 
nitride  and  alumina  grains  using  classical  MD.  Their  work 
reveals  the  presence  of  structural  order  and  composition 
gradients  across  the  films  and  the  segregation  of  impurity 
atoms  at  IGF/grain  interfaces,  which  agree  with  experi¬ 
mental  observations.  Compared  to  MD,  it  is  more  feasible 
to  obtain  thermodynamically  equilibrated  structures 
through  MC  simulations.  Williams  and  Mishin  [139] 
recently  applied  the  semi-grand  canonical  MC  method  with 
an  embedded-atom  potential  to  study  grain  boundary  pre- 
melting  in  binary  Cu-Ag  alloys.  They  observed  the  forma¬ 
tion  of  a  thin  disordered  layer  with  a  composition  close  to 
the  liquid  composition  at  a  £5  grain  boundary  when  the 
temperature  nears  the  solidus  line  from  below.  However, 
the  thickness  of  the  liquid-like  layer  remains  finite  at  the 
solidus  and  the  grain  boundary  can  be  overheated  or  over¬ 
saturated  above  the  solidus,  which  indicates  nonmelting 
behavior. 

Advances  in  computing  power  in  the  last  few  decades 
have  enabled  the  application  of  first-principles  calculations 
to  almost  every  area  of  materials  science  including  grain 
boundary  complexions.  Density  functional  theory  has  been 
employed  to  investigate  various  aspects  of  IGFs,  including 
the  bonding  and  electronic  structures  [140,141],  the  specific 
adsorption  sites  of  rare-earth  elements  at  IGF/grain  inter¬ 
faces  [142-144],  and  the  viscosity  [145]  and  mechanical 
strength  [146-148]  of  intergranular  films.  These  calcula¬ 
tions  provide  quantitative  information  on  the  atomistic 
structures  and  properties  of  IGFs,  and  facilitate  the  inter¬ 
pretation  of  experimental  observations  [144], 

In  addition  to  grain  boundary  premelting  transitions, 
other  types  of  grain  boundary  complexion  transitions  have 
also  been  studied  and  revealed  by  MD  or  MC  simulations. 
A  recent  study  predicted  an  extrinsic  complexion  transition 
in  £5(310)  boundaries  in  Ag-doped  Cu  as  a  function  of 
temperature,  from  a  monolayer  to  a  bilayer  complexion, 
as  shown  in  Fig.  4(d  and  e)  [149].  In  this  case,  the  bilayer 
complexion  contained  less  total  adsorbate  (Ag)  than  the 
monolayer  complexion,  in  contrast  to  the  assumed  general 
trend  of  increasing  adsorbate  content  with  increasing 


complexion  thickness.  Other  examples  include  the  rough¬ 
ening  transition  [150],  the  defaceting  transition  [151],  a  dis¬ 
location  pairing  transition  at  low-angle  grain  boundaries  in 
body-centered  cubic  (bcc)  Fe  [152]  and  transitions  between 
two  ordered  grain  boundary  structures  [39,153],  A  recent 
experimental  study  has  revealed  solid-solid  interface  recon¬ 
struction  at  an  Ni-Al203  boundary  [45],  Solid-solid  inter¬ 
facial  reconstruction  transitions  at  grain  boundaries  and 
heterophase  boundaries  could  potentially  be  a  very  impor¬ 
tant  type  of  complexion  transition.  As  the  study  of  these 
types  of  transitions  is  particularly  difficult  in  a  transmission 
electron  microscope  (because  the  images  are  a  two-dimen¬ 
sional  projection  parallel  to  the  interface),  atomistic  model¬ 
ing  will  play  a  critical  role  in  revealing  the  detailed 
structures  and  stability  of  such  complexions. 

A  continuum  field  theory,  referred  to  as  the  phase-field 
crystal  model,  was  recently  developed  [154,155]  to  describe 
phenomena  on  atomic  length  and  diffusive  time  scales. 
Unlike  the  diffuse-interface  models  introduced  in  Sec¬ 
tion  2.3.2,  a  phase-field  crystal  model  constructs  the  free 
energy  functional  using  the  local-time-averaged  atomic 
density  as  the  order  parameter  and  thus  resolves  the  atomic 
structure  of  a  crystalline  material.  At  the  same  time,  the 
mean-field  approach  allows  it  to  cover  time  scales  much 
larger  than  MD,  which  makes  it  computationally  more  effi¬ 
cient  in  obtaining  equilibrium  configurations  of  grain 
boundary  structure.  With  such  advantages,  new  insights 
on  grain  boundary  complexion  transitions  have  been 
obtained  using  this  modeling  approach.  For  example,  Mel- 
lenthin  et  al.  [156]  performed  a  phase-field  crystal  study  of 
grain  boundary  premelting  in  a  two-dimensional  single¬ 
component  crystal  with  hexagonal  ordering.  A  liquid  film 
was  found  to  form  at  high-angle  grain  boundaries  and  its 
thickness  diverges  logarithmically  when  the  melting  point 
is  reached,  which  is  consistent  with  the  continuous  premelt¬ 
ing  behavior  depicted  by  diffuse-interface  models.  How¬ 
ever,  local  melting  only  starts  around  dislocations  at  low- 
angle  grain  boundaries;  the  grain  boundary  structure  con¬ 
sists  of  isolated  liquid  pools  separated  by  solid  bridges, 
which  can  be  overheated.  At  intermediate  misorientations, 
a  first-order  dislocation  pairing  transition  occurs  below  the 
melting  point,  which  is  also  observed  in  MD  simulation 
[152].  Similar  misorientation-dependent  melting  behavior 
was  also  found  in  a  phase-field  crystal  study  of  a  three- 
dimensional  bcc  crystal  by  Berry  et  al.  [157], 

2.5.  Developing  grain  boundary  complexion  diagrams 

The  development  of  grain  boundary  complexion  dia¬ 
grams  is  a  critical  aspect  of  complexions  research.  Com¬ 
plexion  diagrams  have  the  potential  to  be  as  important 
and  useful  as  bulk  phase  diagrams,  which  are  indispensable 
for  many  aspects  of  materials  science  and  engineering. 
Complexion  diagrams  can  be  developed  based  on 
experimental  data  and/or  thermodynamic  models  and 
simulations,  and  we  discuss  both  approaches  in  this  sec¬ 
tion.  In  general,  constructing  complexion  diagrams  is  very 
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challenging  from  a  scientific  perspective,  but  it  can  also  be 
very  technologically  rewarding,  given  the  important  roles 
of  grain  boundaries  in  controlling  both  materials  fabrica¬ 
tion  and  properties.  Thus,  this  section  focuses  on  discuss¬ 
ing  the  construction  and  potential  applications  of  grain 
boundary  complexion  diagrams. 

First,  we  comment  on  some  general  characteristics  of 
complexion  diagrams.  Because  grain  boundary  complex¬ 
ions  have  6  +  C  degrees  of  thermodynamic  freedom,  a  full 
complexion  diagram,  even  for  a  single-component  material 
(Cm  1),  cannot  be  plotted  in  two  dimensions  as  is  typically 
done  for  bulk  phase  diagrams,  which  have  C  degrees  of 
freedom  (hence,  binary  bulk  phase  diagrams  can  be  plotted 
in  two  dimensions).  However,  two-dimensional  complex¬ 
ion  diagrams  can  be  plotted  in  which  the  axes  represent 
two  of  the  five  grain  boundary  geometrical  parameters, 
as  done,  for  example,  in  Ref.  [9],  In  such  diagrams,  the 
other  three  geometrical  parameters  and  the  bulk  composi¬ 
tion  and  temperature  are  held  constant.  Alternatively,  if 
the  five  macroscopic  geometric  grain  boundary  parameters 
(and  pressure)  are  held  constant,  then  complexion  dia¬ 
grams  will  have  C  degrees  of  freedom  and  can  therefore 
be  overlaid  onto  bulk  phase  diagrams.  For  example,  the 
premelting,  prewetting  and  wetting  lines  shown  in  Fig.  9 
are  lines  of  complexion  transition  that  represent  the  com¬ 
plexion  transitions  of  three  different  grain  boundaries,  each 
with  a  fixed  misorientation  ( AS/,  A0m,  and  A 6S).  Additional 
lines  of  complexion  transition  could  be  drawn  for  other 
grain  boundary  misorientations,  but  such  a  complexion 
diagram  would  be  overly  complicated  and  hence  not  very 
useful. 

When  complexion  diagrams  are  overlaid  onto  bulk 
phase  diagrams,  it  is  helpful  to  emphasize  regions  of  com¬ 
plexion  existence  in  addition  to  lines  of  complexion  transi¬ 
tion,  the  same  way  that  regions  of  bulk  phase  coexistence 
are  marked  on  phase  diagrams  in  addition  to  the  solvus, 
solidus  and  liquidus  lines.  Such  complexion  diagrams  are 
overlaid  onto  a  standard  bulk  phase  diagram  (Fig.  14(a)) 
by  drawing  lines  of  complexion  transition  and  indicating 
regions  of  complexion  stability  (Fig.  14(b)).  In  the  sche¬ 
matic  example  in  Fig.  14(b),  it  is  understood  that,  for  a  real 
sample  with  no  restriction  on  allowed  grain  boundary 
geometry,  not  all  grain  boundaries  in  the  bilayer  region 
actually  have  a  bilayer  complexion;  rather,  some  of  the 
boundaries  have  the  bilayer  complexion  (e.g.  the  grain 
boundaries  with  an  initially  higher  energy),  while  the  other 
grain  boundaries  have  a  monolayer  complexion.  Con¬ 
versely,  in  the  monolayer  region,  some  grain  boundaries 
might  have  the  bilayer  complexion.  This  somewhat  impre¬ 
cise  situation  necessarily  results  from  the  reduction  of  the 
6  +  C  dimensional  space  of  complexions  into  a  subspace 
that  can  be  plotted  in  two  dimensions.  We  additionally 
point  out  that  Fig.  14(b)  is  highly  schematic  and  for  pur¬ 
poses  of  illustration  only;  whether  or  not  all  features 
depicted  in  the  diagram  actually  exist  in  a  real  materials 
system  is  an  open  question  that  has  not  been  sufficiently 
explored  experimentally.  Furthermore,  for  simplicity,  some 


features  have  been  omitted;  for  example,  the  transition 
from  a  monolayer  to  a  clean  complexion  as  XB  — >  0  is 
not  shown,  although  it  is  understood  that  for  a  material 
in  which  XB  =  0  and  X A  =  1  (the  far  left  side  of 
Fig.  14(b))  there  cannot  be  a  monolayer  complexion 
because  there  are  no  solute  atoms. 

2.5.1.  Experimental  approach 

In  practice,  complexion  diagrams  can  be  constructed 
based  on  experimental  data  from  bicrystals  or  polycrystals. 
As  suggested  by  Hart  [5],  complexion  transitions  can  be 
identified  by  discontinuous  changes  in  properties  such  as 
adsorption,  diffusivity,  grain  boundary  mobility  and 
mechanical  properties.  Additionally,  electron  microscopy 
can  be  used  to  identify  the  structure  and  chemistry  of  dif¬ 
ferent  complexions,  as  done,  for  example,  in  Ref.  [17], 
Thus,  the  construction  of  complexion  diagrams  can  pro¬ 
ceed  in  a  manner  analogous  to  the  construction  of  bulk 
phase  diagrams,  in  which  properties  and  structure  are  stud¬ 
ied  as  a  function  of  thermodynamic  parameters  (usually 
temperature  and  bulk  composition)  and  discontinuities 
represent  transitions  between  regions  of  thermodynamic 
stability. 

The  solidus  line  on  bulk  phase  diagrams  (representing  the 
line  of  solid-liquid  bulk  phase  transition)  has  often  been 
determined  by  the  dramatic  and  discontinuous  changes  in 
properties  between  the  solid  and  liquid  phases  [5];  mechani¬ 
cal  properties  testing  is  generally  a  reliable  indicator  of  a 
bulk  phase  transition.  However,  since  grain  boundary  com¬ 
plexion  transitions  can  also  cause  dramatic  changes  in  mac¬ 
roscopic  properties,  changes  in  mechanical  properties  do  not 
always  correspond  to  a  bulk  phase  transition;  sometimes 
they  can  correspond  to  a  complexion  transition.  Where  they 
do,  the  changes  in  mechanical  properties  (or  other  proper¬ 
ties)  can  be  used  to  create  complexion  diagrams. 

To  elucidate  these  ideas  and  give  a  concrete  example  of 
an  experimentally  derived  complexion  diagram,  we  will  dis¬ 
cuss  the  Cu-Bi  system,  a  classic  grain  boundary  embrittle¬ 
ment  system.  The  bulk  solubility  of  Bi  in  Cu  is  extremely 
low  and  the  determination  of  the  solidus  line  has  been 
the  subject  of  many  experimental  studies  (see  e.g.  Ref. 
[158]  and  references  therein).  An  early  study  using 
notched-bar  mechanical  testing  concluded  that  the  solid 
solubility  of  Bi  in  Cu  is  30  ppm  at  800  °C  and  3  ppm  at 
600  °C;  at  lower  Bi  concentrations,  embrittlement  did  not 
occur,  presumably  due  to  the  absence  of  a  bulk  liquid 
phase  [159],  It  is  now  known  that  embrittlement  in  the 
Cu-Bi  system,  just  as  in  the  Ni-Bi  system  [48],  is  caused 
by  a  grain  boundary  complexion  transition  that  results  in 
a  Bi-rich  bilayer  complexion  at  the  grain  boundary  [35], 

A  more  recent  study  in  which  Cu  single  crystals  were 
annealed  in  the  presence  of  a  Bi  vapor  source,  followed 
by  mass  spectrometry  to  measure  the  solid  solubility  of 
Bi,  provides  a  more  accurate  and  complete  solidus  line 
[160].  Studies  have  also  measured  Bi  adsorption  at  Cu  grain 
boundaries  with  Auger  electron  spectroscopy  (AES)  [161] 
and  diffusivity  of  Cu  and  Bi  in  Bi-doped  Cu  [162]  as  a 
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(a)  Phase  Diagram  (b)  Complexion  Diagram 


Fig.  14.  (a)  A  schematic  bulk  phase  diagram  for  a  binary  eutectic  system;  (b)  a  schematic  complexion  diagram  overlaid  onto  the  bulk  phase  diagram  in  (a), 
showing  regions  of  monolayer  and  bilayer  complexion  existence.  Such  complexion  diagrams  represent  an  “average”  grain  boundary  behavior  (because 
they  can  only  be  strictly  true  for  a  single  set  of  grain  boundary  geometrical  parameters.  R  and  h).  It  is  understood  that  some  grain  boundaries  in  the  bilayer 
region  might  actually  have  a  monolayer  complexion  and  vice  versa.  The  complexion  diagram  here  is  for  general  illustration  purposes  only  and  not  all 
features  necessarily  exist  in  a  real  materials  system,  and  additional  features  that  are  not  shown  will  be  present  in  a  real  system.  For  example,  as  the 
composition  tends  toward  XB  =  0  (XA  =  1),  i.e.  toward  the  left  side  of  the  diagram,  the  monolayer  complexion  will  transition  to  a  clean  complexion.  This 
transition  may  be  of  first  or  higher  order,  but  for  simplicity  is  not  shown  on  the  diagram. 


function  of  temperature  and  bulk  concentration.  These 
experimental  data  can  be  used  in  conjunction  with  electron 
microscopy  studies  [35,163]  to  create  a  complexion  diagram 
showing  the  regions  of  monolayer  and  bilayer  complexion 
stability  in  the  Cu-Bi  system.  Selected  diffusion  data  are 
shown  in  Fig.  15(a  and  b),  along  with  microscopy  results 
(Fig.  15(c  and  d))  and  a  map  showing  the  regions  of 
~1  ML  and  ~2  ML  Bi  adsorption  based  on  AES  data 
(Fig.  15(e)).  In  fact,  Fig  15(e)  itself  is  a  complexion  dia¬ 
gram,  but  the  adsorbate  arrangement  and  atomic  structure 
of  the  complexion  was  not  known  at  the  time.  Taken 
together,  these  data  can  be  used  to  plot  the  Cu-Bi  complex¬ 
ion  diagram  shown  in  Fig.  15(f).  Interestingly,  the  onset  of 
embrittlement  in  the  notched-bar  mechanical  test  results 
discussed  above  [159]  (indicated  by  two  yellow  dots  on 
Fig.  15(e))  correlates  better  with  the  complexion  transition 
line  (monolayer-to-bilayer  complexion  transition)  than  it 
does  with  the  bulk  solidus  line,  consistent  with  the  idea  that 
the  formation  of  a  bilayer  complexion  in  Bi-doped  Cu  is 
the  root  cause  of  embrittlement  in  this  system. 

There  is  an  important  caveat  to  note  here:  the  STEM 
image  in  Fig.  15(c)  shows  a  monolayer  complexion  at  a 
high-symmetry  grain  boundary  in  a  bicrystal.  The  grain 
boundary  complexion  at  this  bicrystal  boundary  (and  in 
most  bicrystals,  unless  special  care  is  taken  to  form  a  gen¬ 
eral  boundary)  is  a  special  case  that  may  not  be  representa¬ 
tive  of  general  polycrystalline  grain  boundaries.  In  fact, 
general  grain  boundaries  in  Bi-doped  Cu  exhibit  the  clean 
complexion  at  low  Bi  content  and  the  bilayer  complexion 
at  higher  Bi  content,  but  there  is  no  evidence  for  a  mono- 
layer  complexion  at  a  general  boundary  in  Bi-doped  Cu 
[35],  This  situation  is  similar  to  Bi-doped  Ni,  in  which 
the  clean  and  bilayer  complexions  were  observed  at  general 
grain  boundaries  but  a  monolayer  complexion  was  not 
observed  [48],  The  clean-to-bilayer  complexion  transition 


(with  no  observed  monolayer  complexion)  was  also 
observed  in  a  low-symmetry  bicrystal  experiment  in  the 
Au-doped  Si  system,  as  shown  in  Fig.  13  [34],  Although 
there  is  not  enough  evidence  to  make  a  confident  general¬ 
ization,  it  is  possible  that  the  monolayer  complexion  rarely 
occurs  at  general  grain  boundaries  in  polycrystalline  mate¬ 
rials  even  though  it  is  frequently  observed  in  bicrystals  and 
high-symmetry  grain  boundaries.  More  experimental  work 
on  this  topic  is  needed,  but  the  results  at  hand  suggest  that 
multilayer  adsorption  models  (e.g.  [164])  are  substantially 
more  physically  relevant  than  monolayer  adsorption  mod¬ 
els.  Indeed,  it  is  physically  intuitive  that  each  of  the  two 
atomic  planes  that  terminate  the  adjacent  crystals  at  a  gen¬ 
eral  grain  boundary  will  contain  adsorbed  solute  (as  in  a 
bilayer  complexion)  and  that  the  solute  would  only  be 
localized  to  a  single  atomic  plane  (as  in  a  monolayer  com¬ 
plexion)  if  the  grain  boundary  is  highly  symmetrical  and/or 
if  the  terminating  planes  are  crystallographically  similar  or 
identical,  as  would  occur  in  bicrystals  or  at  special  grain 
boundaries. 

In  Fig.  16,  we  present  a  complexion  diagram  for  the 
Y203-doped  A1203  system  based  on  experimental  results 
from  grain  growth  studies.  Although  this  diagram  is  in  the 
initial  stages  of  development  and  needs  further  work  and 
more  data  for  refinement,  it  illustrates  the  manner  in  which 
a  complexion  diagram  for  a  ceramic  material  can  be  con¬ 
structed.  The  foundation  of  this  diagram  is  the  idea  that 
abnormal  grain  growth  is  caused  by  the  coexistence  of  two 
(or  more)  grain  boundary  complexions  with  dramatically 
different  grain  boundary  mobilities,  as  has  been  shown  to 
exist  in  several  different  doped  alumina  systems  [17]. 

There  are  three  distinct  regions  of  the  complexion  dia¬ 
gram  in  Fig.  16:  a  low-T  complexion  region  with  low  grain 
boundary  mobility  (blue),  a  high-T  complexion  region 
with  high  mobility  (red)  and  a  region  in  which  these  two 
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Fig.  15.  By  consideration  of  a  collection  of  experimental  data  (a-e),  a  grain  boundary  complexion  diagram  can  be  constructed  (e).  (a,  b)  Grain  boundary 
diffusivity  for  Cu  and  Bi  radioisotopes  in  Bi-doped  Cu  as  a  function  of  Bi  content  at  1 1 16  K  (843  °C);  there  is  a  sharp  jump  in  diffusivity  at  the  complexion 
boundary  shown  by  the  dotted  vertical  line,  (c,  d)  STEM-HAADF  images  of  the  Bi-rich  monolayer  complexion  and  bilayer  complexion,  which  correspond 
to  small  and  large  overall  Bi  concentration,  respectively,  (e)  A  plot  of  AES  results  superimposed  on  the  Cu-Bi  phase  diagram;  blue  squares  represent 
approximately  1  ML  of  adsorbed  Bi  and  red  squares  represent  approximately  2  ML  of  Bi  (color  added),  (f)  A  grain  boundary  complexion  diagram 
showing  regions  of  complexion  stability  (green  dots  correspond  to  diffusivity  measurements  shown  in  (a,  b)  and  yellow  dots  correspond  to  embrittlement 
of  notched-bar  test  specimens  from  Ref.  [159]).  The  phrase  "monolayer  or  clean”  in  the  blue  region  means  that  there  is  conflicting  information  about  which 
complexion  is  dominant  in  this  region.  A  study  employing  STEM-HAADF  images  suggests  that  grain  boundaries  at  very  low  Bi  content  typically  have  the 
clean  complexion  [35],  However,  AES  results  show  significant  segregation  in  this  region  of  approximately  1  ML  161],  which  would  be  consistent  with  the 
monolayer  complexion  shown  in  (c),  although  we  note  that  the  grain  boundary  in  (c)  is  a  high-symmetry  bicrystal  boundary  and  might  not  be 
representative  of  most  general  boundaries  in  a  polycrystal.  Parts  (a,  b)  adapted  from  Ref.  [  162];  part  (c)  adapted  from  Ref.  [163];  part  (d)  adapted  from 
Ref.  [35];  part  (e)  adapted  from  Ref.  [161],  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this 
article.) 

complexions  coexist  (purple).  Experiment  has  shown  that  where  a  single  complexion  type  dominates  the  grain  growth 
normal  grain  growth  with  a  unimodal  grain  size  distribu-  behavior.  In  the  region  of  complexion  coexistence,  abnor- 

tion  occurs  in  the  low-T  and  high-T  complexion  regions,  mal  grain  growth  occurs  and  the  grain  size  distribution 
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Al2°3  mol  ppm  Y203 

Fig.  16.  Experimental  complexion  diagram  for  Y2O3-AI2O3  showing 
regions  in  which  a  low-temperature  complexion  with  low  grain  boundary 
mobility  (blue)  and  a  high-temperature  complexion  with  high  mobility 
(red)  are  stable.  Normal  grain  growth  occurs  in  both  regions.  There  is  a 
region  of  overlap  where  it  is  believed  that  both  complexions  may  coexist 
(purple).  It  is  believed  that  abnormal  grain  growth  occurs  in  this  region 
because  the  high-T  complexion  has  significantly  higher  grain  boundary 
mobility  than  the  low-T  complexion.  The  diagram  is  still  in  the  early  stages 
of  development  and  further  experimental  investigation  is  needed  to  refine 
it.  It  was  constructed  from  a  limited  data  set  (10  points)  compiled  from 
various  sources:  data  points  1-3  (yellow)  are  from  Refs.  [165-171],  points 
4—7  (green)  are  from  Ref.  [172]  and  points  8-10  (orange)  are  from  Ref. 
[173],  (For  interpretation  of  the  references  to  color  in  this  figure  legend, 
the  reader  is  referred  to  the  web  version  of  this  article.) 

becomes  bimodal.  The  experimental  data  points  in  the  fig¬ 
ure  (labeled  1-10)  are  gathered  from  several  sources:  data 
points  1-3  (yellow)  are  from  Refs.  [165-171],  points  4-7 
(green)  are  from  Ref.  [172]  and  points  8-10  (orange)  are 
from  Ref.  [173], 

The  observations  corresponding  to  these  10  data  points 
are  summarized  in  Table  1  and  can  be  briefly  described  as 
follows.  Data  points  1,  4,  5  and  6  all  demonstrate  a 
unimodal  grain  size  distribution  and  the  absence  of  a 


second  phase  (YAG),  indicating  that  they  have  a  low-T 
complexion  in  the  single-phase  region  of  the  bulk  phase 
diagram.  Data  points  2  and  8  also  have  a  unimodal  grain 
size  distribution,  and  data  point  2  exhibits  the  presence 
of  a  second  phase  (YAG).  No  second  phase  was  detected 
in  the  experiment  represented  by  data  point  8,  though  a 
second  phase  appeared  at  the  same  doping  level  at  higher 
annealing  temperatures  (i.e.  with  increasing  grain  size), 
indicating  that  the  solvus  line  may  be  grain  size  dependent. 
Abnormal  grain  growth  occurred  in  the  samples  repre¬ 
sented  by  data  points  3,  7  and  9,  suggesting  that  two  com¬ 
plexions  with  dramatically  different  mobilities  exist  in  this 
region.  Of  particular  importance  is  data  point  7,  which 
contains  no  second  phase  and  shows  that  a  complexion 
transition  may  occur  in  the  single-phase  region  of  the  phase 
diagram,  leading  to  abnormal  grain  growth. 

The  two  experimental  complexion  diagrams  discussed  in 
this  section  demonstrate  a  general  approach  to  construct¬ 
ing  experimental  grain  boundary  complexion  diagrams  by 
combining  results  from  several  different  types  of  studies, 
such  as  mechanical  testing,  grain  growth  studies,  chemical 
analysis  and  electron  microscopy.  This  approach  is  analo¬ 
gous  to  the  method  by  which  bulk  phase  diagrams  are  con¬ 
structed,  and  can  be  applied  to  many  other  systems  in 
which  complexion  transitions  play  a  key  role  in  materials 
processing  and  performance. 

2.5.2.  Thermodynamic  calculation  approach 

In  this  section,  we  discuss  efforts  to  calculate  complexion 
diagrams.  Since  few  prior  studies  have  been  conducted  to 
compute  rigorous  grain  boundary  complexion  diagrams 
for  real  systems  that  can  be  verified  by  experiment,  we  will 
discuss  recent  efforts  to  compute  a  type  of  diagram  that  is 
less  rigorous  but  still  practically  useful,  known  as  a  “/.-dia¬ 
gram".  The  calculation  of  /.-diagrams  can  be  viewed  as  a 
first  step  toward  computing  more  rigorous  grain  boundary 
complexion  diagrams  for  real  systems. 

We  first  note  the  existence  of  a  vast  literature  [43,174— 
197]  on  constructing  surface  complexion  diagrams,  which 
are  often  called  “surface  phase  diagrams"  by  the  surface 
physics  community.  These  include  surface  reconstructions 
for  both  pure,  adsorbate-free  surfaces  and  surfaces  with 
adsorbed  impurity  atoms.  However,  the  majority  of  this 
work  has  been  done  for  either  adsorbate-free  surfaces  or 
adsorbate  amounts  in  the  ML  and  sub-ML  regime 
(yet  with  the  observations  of  a  rich  variety  of  surface 


Table  1 

Experimental  data  used  to  construct  the  grain  boundary  complexion  diagram  for  Y-doped  AI2O3  (Fig.  16). 


1350 °C 

1450  °C 

1475  °C 

1500 °C 

1550 °C 

1600 °C 

1700 °C 

100  ppm 

NGG/1 

NGG/1 

NGG/1 

200  ppm 

450  ppm 

1000  ppm 

NGG/1 

NGG/1 

AGG/2 

AGG/2 

NGG/1 

NGG/2 

AGG/1 

The  type  of  grain  growth  (normal  -  NGG  or  abnormal  -  AGG)  and  number  of  phases  observed  (1  or  2)  are  indicated  for  each  level  of  Y  dopant  and 
temperature.  These  data  points  have  been  compiled  front  Refs.  [165-173], 
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complexions).  However,  recent  studies  suggest  the  possibil¬ 
ity  of  thicker  complexions.  A  surface  complexion  diagram 
for  nanoscale  equilibrium-thickness  “surface  amorphous 
films”  (i.e.  the  free-surface  counterpart  to  the  nanolayer 
complexion  at  grain  boundaries)  has  been  constructed 
experimentally  for  Bi203  on  ZnO  (1120)  surfaces 
[23,198-200],  and  a  corresponding  thermodynamic  model 
(analogous  to  the  grain  boundary  complexion  model  [13]) 
has  been  developed  to  construct  surface  complexion  dia¬ 
grams  theoretically  [113],  Complexion  diagrams  for  hetero¬ 
phase  interfaces  have  also  been  constructed  by  first- 
principles  calculations,  e.g.  for  V-doped  WC-Co  in  a  recent 
study  [201].  Lines  of  grain  boundary  complexion  transition 
have  also  been  overlaid  onto  bulk  phase  diagrams  based  on 
experimental  results  [77,202], 

A  few  prior  studies  have  been  conducted  to  compute  rig¬ 
orous  grain  boundary  complexion  diagrams  (with  transition 
lines  and  critical  points)  for  simplified  ideal  systems,  includ¬ 
ing  regular  solution  alloys.  These  studies  have  computed  a 
prewetting  complexion  diagram  for  a  binary  alloy  that 
exhibits  a  solid-state  miscibility  gap  (Fig.  7(a))  [80],  a  pre¬ 
melting  complexion  diagram  for  unary  systems  (Fig.  8(b)) 
[13],  a  coupled  premelting  and  prewetting  complexion  dia¬ 
gram  for  a  binary  alloy  with  a  eutectic  reaction  (Fig.  9) 
[16],  and  a  complexion  diagram  for  a  system  with  a  hard 
sphere  type  potential  that  leads  to  the  formation  of 
Dillon-Harmer  complexions  (Fig.  12(b))  [25],  A  prior 
phase-field  modeling  study  has  computed  a  premelting-like 
complexion  diagram  for  Cu-Ag  [61],  which  has  not  yet  been 
verified  experimentally  (see  more  discussion  subsequently). 

Since  prior  studies  to  compute  rigorous  grain  boundary 
complexion  diagrams  that  can  be  experimentally  verified 
are  lacking,  some  recent  efforts  have  been  made  to  develop 
2-diagrams.  These  diagrams,  which  do  not  have  well- 
defined  transition  lines,  essentially  treat  the  grain  boundary 
as  a  nanoscale,  liquid-like  state  of  interfacial  matter,  simi¬ 
lar  to  the  nanolayer  complexion  shown  in  Fig.  10.  Conse¬ 
quently,  it  has  been  proven  that  such  computed  2- 
diagrams  can  effectively  predict  the  thermodynamic  ten¬ 
dency  for  general  grain  boundaries  to  disorder 
[22,31,32,105],  Although  they  are  not  rigorous  complexion 
diagrams,  the  ability  for  2-diagrams  to  predict  some  useful 
trends  is  demonstrated  by  the  prior  studies  of  solid-state 
activated  sintering  of  oxide  ceramics  [47]  and  refractory 
metals  [22,31,32,106,203-205], 

Solid-state  activated  sintering  refers  to  the  enhancement 
of  sintering  rates  due  to  solid-state  additives  (prior  to  the 
formation  of  a  bulk  liquid  phase).  After  its  discovery  in 
the  1950s,  the  exact  underlying  mechanism  for  subeutectic 
activated  sintering  puzzled  the  materials  community  for 
decades.  Recent  studies  showed  that  bulk  phase  diagrams 
are  not  adequate  for  designing  optimal  activated  sintering 
protocols.  This  is  because  liquid-like  complexions  (nano¬ 
layer  complexions)  can  form  at  thermodynamic  equilibra¬ 
tion  well  below  the  bulk  solidus  line  and  enhanced 
diffusion  in  these  complexions  can  result  in  enhanced 
sintering  rates,  similar  to  the  case  of  liquid-phase  sintering 


but  occurring  under  conditions  where  the  bulk  liquid  is  not 
stable.  Specifically,  formation  of  liquid-like  grain  boundary 
complexions  and  enhanced  sintering  can  initiate  below  60% 
of  bulk  solidus  temperatures.  Recent  HRTEM  studies  have 
directly  revealed  the  stabilization  of  liquid-like  nanolayer 
complexions  well  below  the  bulk  solidus  line  in  both  cera¬ 
mic  [23,47,63,206]  and  metallic  [203,207]  activated  sintering 
systems. 

To  explain  why  liquid-like  nanolayer  complexions  can 
form  below  the  bulk  solidus  line,  we  should  emphasize  that 
these  films  should  not  follow  the  same  thermodynamic 
relations  as  bulk  phases  in  general.  It  is  now  well  estab¬ 
lished  that  nanoparticles  can  often  melt  at  hundreds  of 
degrees  below  the  corresponding  bulk  melting  temperature 
(Fig.  17(b))  [208],  It  is  worth  noting  that  the  suppression  of 
the  melting  point  of  the  nanoparticles  and  surface  premelt¬ 
ing  (which  is  a  complexion  transition;  see  Fig.  17(a)  for  an 
example)  are  two  related,  but  not  identical,  phenomena. 
Following  the  same  concept,  Tanaka  et  al.  [209-211]  devel¬ 
oped  phase  diagrams  of  binary  nanoparticles  showing  that 
the  solidus  and  liquidus  lines  are  size  dependent 
(Fig.  17(c)).  As  a  first-order  approximation,  we  can  view 
the  liquid-like  nanolayer  grain  boundary  complexions  as 
interfacial  films  that  are  confined  between  two  grains;  it 
is  thus  not  surprising  that  they  can  be  stabilized  at  temper¬ 
atures  well  below  the  bulk  solidus  line  (Fig.  17(d)). 

Pursuing  along  this  line,  grain  boundary  2-diagrams 
have  been  calculated  to  forecast  the  stability  of  liquid-like 
grain  boundary  complexions  via  extending  bulk  CalPhaD 
(Calculation  of  Phase  Diagram)  methods  to  grain  bound¬ 
aries  [22,31,32,106].  The  basic  idea  is  as  follows.  A  nano¬ 
meter-thick  quasi-liquid  complexion  can  be 
thermodynamically  stabilized  at  an  average  general  grain 
boundary  if  the  excess  free  energy  of  this  “dry”  grain 
boundary  (y^;  noting  that  yGB  is  reserved  for  the  true, 
thermodynamically  equilibrated,  grain  boundary  energy 
that  corresponds  to  the  complexion  that  has  the  lowest 
energy,  i.e.  a  complexion  with  the  “equilibrium”  thickness 
and  a  global  minimum  in  Gx(h)  in  this  model)  is  greater 
than  the  sum  of  two  solid-liquid  interfacial  energies  (2 ysj): 

Ay  =  2ys/ -  y£]  <  0  (13) 

In  the  sharp-interface  (or  force-balance)  model  dis¬ 
cussed  in  the  prior  section  (extended  from  premelting  mod¬ 
els  for  unary  systems  [63,70,7 1 ,212]),  the  relative  interfacial 
energy  can  be  expressed  as: 

G\h)  -  y<°>  =  Ay  +  A G™*  ■  h  +  AG„,V(A)  +  W 

ee  Ay /(/*)  + A G(Ziph-h  (14) 

where  h  is  the  film  thickness  and  A G^rph  is  the  volumetric 
free-energy  penalty  for  forming  an  undercooled  liquid 
(which  can  be  quantified  using  bulk  CalPhaD  methods). 
It  is  important  to  recognize  that  this  liquid-like  complexion 
is  not  a  true  liquid  as  a  partial  order  and  a  through-thick¬ 
ness  gradient  in  the  order  parameter  inevitably  exist.  Thus, 
the  term  AGn  y(h)  can  be  introduced  to  represent  the  excess 
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Fig.  17.  (a)  Premelting  on  the  surface  of  an  ice  crystal  (reprinted  from  Ref.  [345]  with  permission)  and  (b)  suppression  of  the  melting  temperature  of 
nanoparticles  (adapted  from  Ref.  [346])  are  two  related,  but  different,  interfacial  phenomena,  (c)  Tanaka  et  al.  [209-211]  developed  phase  diagrams  of 
binary  nanoparticles,  showing  that  the  solidus  and  liquidus  lines  are  size  dependent;  adapted  from  Ref.  [210],  (d)  Comparison  between  a  grain  boundary 
/.-diagram  [22,31,32,106]  for  the  Ni-W  system  and  an  experimental  HRTEM  image  of  a  nanolayer  complexion  at  a  grain  boundary  in  Ni-doped  W; 
reprinted  front  Ref.  [32]  with  permission. 


free  energies  associated  with  partial  order  and  gradient 
energies,  such  as  those  described  by  the  integrals  in 
Eqs.  (9)  and  (10)  in  the  diffuse-interface  approach.  The 
term  a*nt  represents  all  other  interfacial  interactions,  such 
as  the  oscillatory  structural  solvation  interaction  that 
produces  the  Dillon-Harmer  complexions,  as  discussed 
below,  as  well  as  dispersion  and  electrostatic  interactions 
in  ceramic  systems  -  noting  that  A G,,y{h)  also  represents 
an  interfacial  interaction,  even  though  it  is  listed  separately 
here  for  clarity.  The  boundary  conditions  for  the  above 
equation  dictate  that  (i)  both  AGriy(h)  and  a*int  vanish  as 
h  — ►  +oo  (noting  that  here  we  select  2ysh  as  our  reference 
state  and  the  excess  free  energies  associated  with  partial 
order  and  gradient  energy  terms  for  the  case  of  h  — *  +oo 
are  already  included  in  2ys!;  thus,  by  definition, 
AG,;,v(+  oo )  =  0  [63])  and  (ii)  Gv(0)  =  Thus,  we  can 
define  a  dimensionless  “interfacial  coefficient”, 
f(h)  =  1  +  [A Gny(h)  +  <j*nt  (h)]/Ay,  to  describe  the  details 
of  all  interfacial  interactions;  by  definition,  it  satisfies  the 
boundary  conditions  /(0)  =  0  and  /(+oc)  =  1.  A  detailed 
discussion  about  the  above  formalism,  partial  order/gradi- 
ent  energies  and  other  interfacial  energy  terms,  as  well  as 
the  appropriate  reference  states  for  writing  the  above  equa¬ 
tion  consistently,  can  be  found  in  Ref.  [63]  (noting  that 
A Gny(h)  is  one  of  interfacial  interactions  in  Ref.  [63], 
though  this  term  is  listed  separately  here  to  emphasize 
the  importance  of  partial  order  and  gradient  terms). 


A  thermodynamic  variable,  X,  is  then  defined  to  repre¬ 
sent  the  maximum  thickness  of  a  quasi-liquid  IGF  that 
can  be  stabilized  at  an  average  grain  boundary  (without 
consideration  of  the  details  of  partial  order,  gradient  terms, 
and  other  interfacial  forces): 

X  =  Mux  [  Ay/AGamorph}  (15) 

We  argue  that  1  can  be  rigorously  defined  thermody¬ 
namically  using  a  “reference”  state  of  the  true  liquid;  it 
explicitly  represents  the  maximum  thickness  of  an  ideal 
liquid  film  that  can  be  stabilized  at  the  grain  boundary 
under  the  conditions  that  the  terms  A Gny(h)  and  o*int  {h) 
can  be  ignored,  which  is  rigorously  held  at  the  limit  of 
h  — >  +oo.  The  actual  interfacial  width  should  correspond 
to  the  “equilibrium”  thickness  (hE e)  that  produces  a  mini¬ 
mum  in  the  excess  grain  boundary  energy.  For  simple 
(structureless)  metallic  alloys  where  an  exponentially 
decaying  interfacial  force  is  dominant,  liquid-like  grain 
boundary  complexions  start  to  develop  at  2  >  ^  with  an 
effective  interfacial  width  of  hEQ  «  £  •  ln(z/c),  where  £  is  a 
coherence  length.  Moreover,  an  oscillatory  structural  inter¬ 
action  can  produce  discrete  interfacial  phases,  as  we  have 
discussed  above  (Fig.  11).  For  ceramics,  the  existence  of 
electrostatic  space  charges  and  London  dispersion  forces 
will  further  complicate  the  situation. 

Nonetheless,  the  X  value  represents  the  thermodynamic 
tendency  to  stabilize  a  liquid-like  complexion  at  an  average 
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general  grain  boundary,  and  it  scales  the  actual  interfacial 
width  (although  the  predictability  of  grain  boundary  1-dia¬ 
grams  will  be  less  for  thinner  complexions  with  significant 
level  of  partial  order  and  discrete  thickness).  Thus,  we  can 
quantify  1  and  plot  lines  of  constant  1  in  bulk  phase  dia¬ 
grams  to  construct  grain  boundary  1-diagrams  -  a  (less  rig¬ 
orous  but  robustly  useful)  type  of  grain  boundary 
complexion  diagram.  Examples  are  shown  in  Fig.  17(d) 
and  Fig.  18(a).  The  correctness  and  usefulness  of  these  grain 
boundary  1-diagrams  in  predicting  general  trends  has  been 
systematically  validated  by  the  following  experiments: 

•  The  model  predictions  were  corroborated  with  direct 
HRTEM/Auger  studies  for  several  selected  tempera¬ 
ture/composition  points  in  W-Ni  and  Mo-Ni,  where 
the  majority  component,  i.e.  the  solvent,  is  underlined 

[31,106,203,205,207], 

•  In  a  systematic  model-experiment  comparison  study 
conducted  for  Mo-Ni,  the  estimated  grain  boundary 
diffusivities  as  a  function  of  temperature  and  overall 
composition  correlated  well  with  the  computed  binary 
l-diagram  (Fig.  18)  [31]. 

•  The  computed  grain  boundary  1-diagrams  (with  no  free 
parameters)  correctly  predicted  the  onset  of  activated 
sintering  temperatures  for  W-Ni,  W-Fe,  W  Co,  W-Pd, 


W-Cu,  and  Mo-Ni  22,31,106],  These  results,  along  with 
direct  HRTEM  observations,  demonstrated  that 
enhanced  diffusion  in  nanolayer  complexions  is  the  root 
cause  for  solid-state  activated  sintering  in  refractory  met¬ 
als,  thereby  solving  a  greater  than  50-year-old  scientific 
mystery  [22,203], 

•  The  predicted  grain  boundary  solidus  temperature  was 
consistent  with  a  prior  direct  grain  boundary  diffusivity 
measurement  for  W-Co  using  radioactive  tracers 

[106,213], 

•  Computed  grain  boundary  1-diagrams  were  used  to  rank 
the  effectiveness  of  seven  sintering  aids  (Ni,  Co,  Fe,  Cr, 
Zr,  Nb  and  Pd)  on  the  densification  of  Mo-Si-B  alloys 

[214], 

•  Most  interestingly,  a  counterintuitive  phenomenon  of 
decreasing  grain  boundary  diffusivity  with  increasing 
temperature  predicted  by  the  computed  1-diagram  of 
the  Mo-Ni  system  was  subsequently  corroborated 
experimentally  for  a  Mo  +  0.5  at.%  Ni  alloy  [204], 

It  should  be  emphasized  again  that  these  grain  boundary 
1-diagrams  are  not  yet  rigorous  grain  boundary  complexion 
diagrams,  which  should  have  well-defined  grain  boundary 
transition  lines  and  critical  points  (similar  to  those  shown 
in  Fig.  12(b)).  However,  they  have  proven  to  be  robustly 
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Fig.  18.  (a|  A  computed  grain  boundary  A-diagram  for  Ni-doped  Mo  can  predict  the  general  trends  of  (b)  the  temperature-  and  composition-dependent 
grain  boundary  diffusivities  estimated  from  sintering  experiments.  Reprinted  from  Ref.  [31]  with  permission  from  APS. 
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useful  in  predicting  some  important  trends.  Recent  studies 
on  constructing  binary  1-diagrams,  along  with  efforts  to 
develop  other  types  of  complexion  diagrams,  are  reviewed 
by  Luo  [32],  Current  efforts  are  in  progress  to  compute  such 
grain  boundary  2-diagrams  for  multicomponent  materials 
and  to  further  consider  boundary-to-boundary  variations; 
these  efforts  will  make  grain  boundary  2-diagrams  a  more 
practically  useful  materials  design  tool. 

Future  studies  should  also  be  conducted  to  further 
develop  more  rigorous  grain  boundary  complexion  dia¬ 
grams  with  well-defined  transition  lines  and  critical  points; 
in  fact,  a  true  grain  boundary  complexion  diagram  with  a 
premelting  transition  line  and  critical  point  has  already  been 
calculated  for  Cu-Ag  based  on  the  CalPhaD  data  using  a 
phase-field  model  [61];  however,  the  model  has  unknown  fit¬ 
ting  parameters  and  experimental  validation  has  not  yet 
been  conducted.  Nonetheless,  this  study  signaled  a  good 
start  to  developing  theoretical  and  modeling  approaches  to 
compute  rigorous  grain  boundary  complexion  diagrams. 

Bulk  phase  diagrams  are  arguably  one  of  the  most 
important  materials  design  tools  for  guiding  the  devising 
of  fabrication  protocols  and  predicting  materials  proper¬ 
ties.  We  believe  that  both  rigorous  grain  boundary  com¬ 
plexion  diagrams  and  some  simplified  grain  boundary 
diagrams  (such  as  the  2-diagrams  discussed  above)  can  be 
equally  useful  materials  design  tools  and  make  it  possible 
to  speed  the  development  of  new  materials.  Thus,  develop¬ 
ing  such  diagrams  is  consistent  with  the  goals  of  the  Mate¬ 
rials  Genome  initiative  [215], 

3.  Characterization  of  complexions 

3.1.  Overview 

This  section  will  describe  techniques  that  can  be  used  to 
characterize  the  structure  and  composition  of  grain  bound¬ 
ary  complexions.  Broadly  speaking,  the  techniques  can  be 
divided  into  direct  characterization  techniques,  which  per¬ 
mit  visualization  of  the  grain  boundary  structure  and 
chemistry  as  well  as  quantification  of  compositional  pro¬ 
files  across  the  grain  boundary,  and  indirect  characteriza¬ 
tion  techniques,  which  measure  changes  in  properties  that 
depend  upon  grain  boundary  structure  and  chemistry, 
e.g.  mobility,  diffusivity,  cohesive  strength,  and  electrical 
and  thermal  properties.  No  single  method,  whether  direct 
or  indirect,  can  fully  characterize  a  grain  boundary  com¬ 
plexion.  The  best  approach  to  characterizing  grain  bound¬ 
ary  complexion  transitions  is  to  use  a  combination  of 
techniques. 

3.2.  Direct  characterization  of  complexions 

Direct  characterization  methods,  such  as  Auger  electron 
spectroscopy  (AES),  transmission  electron  microscopy 
(TEM)  and  atom  probe  tomography,  can  provide  informa¬ 
tion  about  the  structure  and  chemistry  of  grain  boundaries. 
AES  is  a  well-established  technique  capable  of  measuring 


the  chemical  composition  of  grain  boundaries  [216,217]. 
TEM  is  a  versatile  and  powerful  method  for  characterizing 
the  structure  and  chemistry  of  grain  boundaries,  and  can 
provide  atomic-level  image  resolution  as  well  as  elemental 
and  chemical  analysis  with  the  complementary  techniques 
of  energy-dispersive  X-ray  spectroscopy  (EDS)  and  elec¬ 
tron  energy  loss  spectroscopy  (EELS).  The  advent  of  aber¬ 
ration-corrected  STEM  in  particular  has  enabled 
groundbreaking  discoveries  in  grain  boundary  complex¬ 
ions,  and  atom  probe  tomography  [218-220]  is  yet  another 
powerful  analytical  technique  which  has  been  applied  to 
study  grain  boundary  chemistry  and  morphology. 

If  a  grain  boundary  can  be  exposed  by  intergranular 
fracture,  the  composition  of  the  exposed  grain  boundary 
region  can  be  studied  by  AES  [207,216],  AES  can  also  be 
used  to  study  the  change  in  composition  as  a  function  of 
depth,  and  thereby  determine  the  relative  surface  enrich¬ 
ment  of  an  element  that  segregates  to  the  grain  boundary. 
Using  sputtering,  AES  has  been  employed  to  demonstrate 
that  multilayer  adsorption  exists  at  grain  boundaries 
[216].  AES  has  also  been  used  in  combination  with  orienta¬ 
tion  imaging  microscopy  (OIM)  in  scanning  electron 
microscopy  (SEM)  to  study  the  anisotropy  of  grain  bound¬ 
ary  segregation  [164,217],  However,  because  AES  can  only 
be  applied  to  study  grain  boundaries  that  fail  by  brittle 
intergranular  fracture,  it  is  best  suited  to  the  study  of  mate¬ 
rials  systems  that  display  some  degree  of  grain  boundary 
embrittlement  of  the  type  that  occurs  via  grain  boundary 
segregation  of  specific  impurity  elements.  Thus  AES  is 
most  useful  for  studying  grain  boundaries  in  multi-compo¬ 
nent  systems  and  additionally  has  a  selection  bias  for  the 
weakest  grain  boundaries. 

TEM  enabled  the  first  direct  observations  of  grain 
boundary  complexions  and  is  still  one  of  the  most  useful 
techniques  for  studying  them.  Referred  to  at  the  time  as 
IGFs,  the  first  grain  boundary  complexions  identified  by 
TEM  were  non-crystalline  films  of  about  1  nm  thickness 
found  at  grain  boundaries  in  Si3N4  [221,222],  The  discov¬ 
ery  of  IGFs  was  a  seminal  achievement  in  grain  boundary 
research.  The  presence  of  the  intergranular  film  was 
inferred  from  the  lack  of  lattice  fringes  in  a  narrow  region 
surrounding  the  boundary  core.  Determining  the  exact 
structure  of  this  non-crystalline  region  remains  an  experi¬ 
mental  challenge.  The  thickness  of  this  grain  boundary 
complexion  was  constant  along  the  boundary,  leading 
Clarke  and  Thomas  [221]  to  conclude  that  it  was  in  equilib¬ 
rium  with  its  surroundings.  Furthermore,  the  IGFs  were 
only  observed  at  high-angle  grain  boundaries  [221,222], 
No  IGFs  were  discovered  at  low-angle  grain  boundaries, 
and  this  absence  is  consistent  with  the  idea  that  complexion 
transitions  occur  preferentially  at  high-energy  grain  bound¬ 
aries.  Although  a  correlation  between  the  amount  of  segre¬ 
gated  impurity  and  the  grain  boundary  energy  is  not 
expected,  the  barrier  for  first-order  complexion  transitions 
could  be  lower  at  higher-energy  boundaries.  IGFs  of  about 
1  nm  thickness  comprising  mostly  Co  were  also  discovered 
at  tungsten  carbide  grain  boundaries  in  WC/Co  composites 
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[223].  In  this  case,  the  carbide  boundaries  could  be  “dry”, 
could  have  a  1  nm  thick  IGF  or  could  have  a  much  thicker 
wetting  layer.  IGFs  have  also  been  imaged  in  metals  such 
as  Ni-doped  W  [207],  In  these  early  studies,  little  structural 
information  could  be  extracted  because  the  IGFs  were  non¬ 
crystalline.  What  could  be  extracted  was  the  thickness  of 
the  layer  and  some  composition  information.  For  example, 
EDS  and  EELS  can  be  used  to  determine  the  elemental 
constituents  in  the  IGF  [224,225,226], 

Certain  requirements  must  be  met  so  that  a  grain 
boundary  complexion  can  be  imaged  in  the  transmission 
electron  microscope,  and  these  requirements  severely 
restrict  the  fraction  of  grain  boundaries  in  a  specimen  that 
can  be  studied.  The  requirements  vary  depending  upon 
whether  image  contrast  arises  predominantly  from  phase 
contrast,  diffraction  contrast  or  Z-contrast,  for  example. 
One  requirement  that  applies  to  virtually  all  TEM  imaging 
and  contrast  modes  is  that  the  grain  boundary  be  in  the 
“edge-on”  condition,  such  that  the  grain  boundary  plane 
is  parallel  to  the  optical  axis  of  the  transmission  electron 
microscope.  As  pointed  out  by  Lou  et  al.  [222],  a  grain 
boundary  complexion  of  2.5  nm  width  in  a  TEM  sample 
of  50  nm  thickness  will  be  obscured  by  the  adjacent  grains 
if  it  is  tilted  more  than  1.5°  away  from  the  exact  edge-on 
condition.  The  edge-on  condition  is  therefore  rather  strin¬ 
gent.  It  can  be  achieved  by  tilting  the  TEM  sample  in  the 
microscope  or  by  searching  for  a  grain  boundary  that  is 
already  in  the  edge-on  condition.  Highly  curved  grain 
boundaries,  i.e.  boundaries  that  have  atomic  steps  sepa¬ 
rated  by  tens  of  nanometers  or  less,  may  be  impossible  to 
image  in  the  edge-on  condition  even  if  the  TEM  sample 
is  extremely  thin.  Thus,  the  edge-on  requirement  introduces 
a  selection  bias  that  favors  atomically  flat  grain  bound¬ 
aries.  While  the  edge-on  requirement  generally  applies  to 
all  S/TEM  (where  S/TEM  refers  to  scanning  and  conven¬ 
tional  TEM)  imaging  modes,  it  is  possible  to  relax  this 
requirement  somewhat  with  the  aid  of  computer  simula¬ 
tions.  For  example,  if  at  least  one  of  the  two  adjacent 
grains  can  be  aligned  to  a  low-index  zone  axis,  HAADF- 
STEM  images  of  curved  grain  boundaries  with  through¬ 
thickness  atomic  steps  may  be  interpreted  with  the  aid  of 
image  simulations  if  adsorption  of  a  high-Z  element  occurs 
on  the  grain  boundary  [227], 

A  second  requirement  that  applies  especially  to  phase 
contrast  imaging  of  grain  boundary  complexions  is  that 
the  adjacent  grains  must  be  aligned  so  that  low-index 
atomic  planes  are  parallel  to  the  optical  axis  of  the  trans¬ 
mission  electron  microscope  such  that  lattice  fringes  can 
be  observed  [222],  This  orientation  permits  the  identifica¬ 
tion  of  a  thick,  disordered  grain  boundary  complexion 
(e.g.  an  IGF)  due  to  the  absence  of  lattice  fringe  contrast 
in  the  grain  boundary  core.  Ideally,  both  grains  should 
be  aligned  to  low-index  zone  axis  conditions,  but  alignment 
such  that  at  least  one  set  of  lattice  planes  is  visible  in  each 
grain  is  usually  sufficient  to  recognize  a  thicker  disordered 
complexion,  such  as  an  IGF.  This  combination  of  require¬ 
ments  is  highly  restrictive  when  one  considers  that  the  grain 


boundary  itself  must  also  be  in  the  edge-on  condition. 
These  restrictions  prohibit  a  large  number  of  general, 
high-angle  grain  boundaries  in  a  polycrystalline  material 
from  being  analyzed,  and  thus  introduce  a  selection  bias 
towards  special  grain  boundaries  with  a  small  population 
in  the  polycrystal.  Since  general  high-angle  grain  bound¬ 
aries  tend  to  have  the  highest  grain  boundary  energy  and 
thus  are  most  likely  to  exhibit  complexion  transitions,  this 
selection  bias  against  them  may  prohibit  important  grain 
boundary  complexions  from  being  identified  even  though 
they  exist  in  the  materials  system. 

The  stringent  requirements  for  high-resolution  TEM 
imaging  of  grain  boundaries  has  led  to  grain  boundary 
studies  being  frequently  done  on  carefully  oriented  bicrys¬ 
tals  with  relatively  high  symmetry.  These  bicrystals  are 
often  fabricated  with  symmetrical  tilt  boundaries  so  that 
the  crystals  have  a  common  zone  axis,  making  it  possible 
to  visualize  atomic  columns  in  both  crystals.  While  these 
bicrystal  studies  are  useful  as  model  systems,  such  ideal 
grain  boundaries  represent  only  a  small  minority  of  grain 
boundaries  in  real  polycrystalline  materials.  To  character¬ 
ize  the  structures  of  complexions,  it  is  necessary  to  target 
the  most  general  boundaries  that  are  difficult  to  image  in 
the  S/TEM  due  to  the  inability  to  achieve  the  edge-on  con¬ 
dition  while  simultaneously  aligning  the  adjacent  grains  to 
low-index  zone  axes.  Crystallographically  “general” 
boundaries  are  most  likely  to  have  thicker,  more  disordered 
complexions,  and  these  boundaries  are  rarely  the  subject  of 
high-resolution  TEM  studies. 

The  recent  availability  of  aberration  correction  for  both 
TEM  and  STEM  has  made  a  dramatic  impact  on  the  study 
of  grain  boundary  complexions.  For  example,  with  aberra¬ 
tion-corrected  STEM,  it  has  been  shown  that  some  degree 
of  order  exists  in  IGFs  in  rare-earth-doped  Si3N4  [144,228]. 
HAADF-STEM  imaging  (i.e.  Z-contrast  imaging)  demon¬ 
strated  that  the  dopants  segregate  preferentially  and  peri¬ 
odically  to  the  interface  between  the  IGF  and  the 
adjacent  crystallite.  Similarly,  aberration-corrected  TEM 
has  demonstrated  that  a  degree  of  order  exists  in  IGFs  at 
gold-sapphire  interfaces  [43],  Both  of  these  electron 
microscopy  studies  were  correlated  with  computer  simula¬ 
tions  that  predicted  the  observed  ordering  within  the  IGF 
[43,144],  Atomistic  computer  simulations  are  becoming 
increasingly  important  in  the  context  of  high-resolution 
S/TEM  because  they  can  dramatically  aid  in  image  inter¬ 
pretation  and  can  provide  meaningful  insight  into  the 
mechanisms  governing  grain  boundary  complexion 
formation. 

It  was  very  challenging,  and  in  some  cases  impossible,  to 
directly  visualize  the  thinner  multilayer  grain  boundary 
complexions  (e.g.  the  bilayer  complexion)  until  aberra¬ 
tion-corrected  STEM  became  available.  For  example,  sin¬ 
gle-  and  multilayer  complexions  were  imaged  in  doped 
alumina  by  Dillon  and  Harmer  [229]  (Fig.  10),  in  Bi-doped 
Ni  by  Luo  et  al.  [48]  (Fig.  10  and  Fig.  19)  and  in  Au-doped 
Si  by  Ma  et  al.  [34]  (Fig.  13).  Dillon  and  Harmer  could  not 
readily  distinguish  between  different  amounts  of  multilayer 
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adsorption  in  the  A1203  system  with  standard  HRTEM 
(corresponding  to  a  clean  complexion,  a  monolayer  com¬ 
plexion  and  a  bilayer  complexion)  [17].  These  complexions 
could  only  be  distinguished  -  and  their  importance  recog¬ 
nized  -  by  the  use  of  aberration-corrected  HAADF-STEM 
imaging.  Similarly,  bilayer  segregation  of  Bi  at  Ni  grain 
boundaries  was  identified  with  aberration-corrected 
HAADF-STEM  [48].  The  segregation  of  Bi  into  a  bilayer 
complexion  and  its  role  in  liquid  metal  embrittlement  in 
the  Ni  Bi  system  had  not  been  previously  recognized 
because  HRTEM  imaging  cannot  show  the  Bi  bilayer  com¬ 
plexion.  Furthermore,  the  Ni-Bi  results  highlight  another 
major  advantage  of  the  aberration-corrected  STEM 
method,  which  is  the  ability  to  study  complexions  at  gen¬ 
eral,  high-angle  grain  boundaries  without  the  necessity  of 
having  the  adjacent  crystallites  aligned  to  low-index  zone 
axis  conditions.  This  is  possible  because  Z-contrast  rather 
than  phase  contrast  is  the  main  contrast  mechanism  in 
HAADF-STEM  imaging.  Thus,  grain  boundary  adsorp¬ 
tion  of  higher-Z  elements  can  be  easily  identified  at  grain 
boundaries  even  if  the  adjacent  crystallites  are  not  properly 
aligned  for  HRTEM  imaging.  The  grain  boundary,  how¬ 
ever,  must  still  be  aligned  to  the  edge-on  condition. 

Despite  the  advances  in  aberration  correction  and  other 
TEM  technology,  grain  boundary  complexions  are  almost 
always  studied  at  room  temperature,  far  away  from  the 
equilibrium  conditions  under  which  they  formed.  Many 
important  grain  boundary  complexions  are  only  stable  at 
high  temperatures,  yet  characterization  is  often  carried 
out  at  room  temperature.  Thus,  the  specimen  must  be 
quenched  to  room  temperature  quickly  enough  to  preserve 
the  high-temperature  complexion. 


It  is  possible  to  study  complexions  under  the  conditions 
at  which  they  are  in  equilibrium,  e.g.  at  high  temperatures 
in  the  TEM,  if  an  in  situ  heating  stage  is  employed.  One 
such  study  demonstrated  that  the  thickness  of  IGFs  in 
Si3N4  changes  as  a  function  of  temperature  [230].  The 
width  of  the  IGF  was  larger  at  950  °C  than  at  room  tem¬ 
perature,  and  hysteresis  was  observed  in  the  IGF  width 
upon  cooling.  These  interesting  results  demonstrate  that 
in  situ  heating  experiments  are  promising,  and  more  efforts 
should  be  aimed  at  studying  grain  boundary  complexions 
in  the  TEM  at  high  temperatures. 

Atom  probe  tomography  combines  a  field  ion  micro¬ 
scope  with  a  mass  spectrometer  to  analyze  composition 
and  structure  on  the  near-atomic  scale  [218-220].  As  atoms 
are  ionized  and  ejected  from  the  sample,  their  trajectories 
are  used  to  determine  their  original  positions  in  the  sample 
and  their  atomic  masses  are  determined  by  mass  spectros¬ 
copy.  Technical  advances  have  led  to  tremendous  increases 
in  the  rate  of  data  collection.  With  these  latest  advances,  it 
is  possible  to  image  ~108  atoms  in  a  period  of  hours, 
approximately  50%  of  all  of  the  ionized  atoms  are  ana¬ 
lyzed,  the  composition  sensitivity  is  in  the  parts  per  million 
range  and  spatial  resolution  is  on  the  order  of  1  nm.  Atom 
probe  tomography  was  used  to  measure  the  grain  bound¬ 
ary  segregation  in  nanocrystalline  Ni-W  alloys  [231]  and 
in  Cu-Ni  alloys  [232],  In  the  Cu-Ni  alloys,  the  grain 
boundary  width  (i.e.  the  width  of  the  segregation  zone  in 
the  grain  boundary)  increases  monotonically  as  a  function 
of  annealing  temperature,  increasing  from  0.7  nm  at  563  K 
up  to  2.5  nm  at  643  K  [232],  with  what  appears  to  be  an 
inflection  point  or  perhaps  a  discontinuity  in  the  plot  of 
thickness  vs.  temperature  at  about  600  K.  The  increase  in 


Fig.  19.  HAADF-STEM  images  of  a  Bi-rich  bilayer  complexion  on  general  grain  boundaries  in  polycrystalline  Ni.  In  (A)  the  grain  boundary  is  still  intact, 
but  (B)  and  (C)  show  the  initiation  of  brittle  fracture.  Reprinted  from  Ref.  [48]  with  permission. 
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the  thickness  of  the  segregation  zone  with  temperature  is 
characteristic  of  a  transition  to  a  thicker  complexion  and 
counter  to  what  is  observed  for  conventional  segregation. 
While  it  is  not  clear  if  these  changes  in  grain  boundary 
width  represent  a  complexion  transition,  this  study  clearly 
highlights  the  unique  capabilities  of  atom  probe  tomogra¬ 
phy  and  the  promise  it  shows  for  studying  grain  boundary 
complexions. 

It  is  often  necessary  to  combine  two  or  more  direct  char¬ 
acterization  techniques  to  gain  a  better  understanding  of  the 
structure  and  chemistry  of  grain  boundary  complexions. 
Each  technique  has  benefits  and  drawbacks,  and  none  can 
provide  complete  information  about  grain  boundary  struc¬ 
ture  and  chemistry.  For  example,  atom  probe  tomography 
has  been  combined  with  TEM  imaging  to  study  grain  bound¬ 
ary  segregation  in  Fe-Si  alloys  [233],  Furthermore,  the  com¬ 
bination  of  atomistic  computer  simulations  with  direct 
characterization  techniques  has  become  increasingly  impor¬ 
tant  in  recent  years  for  interpreting  experimental  results  and 
understanding  the  mechanisms  behind  grain  boundary  com¬ 
plexion  formation  and  stability  [43,144],  The  most  fruitful 
studies  of  grain  boundary  complexions  in  the  future  will 
likely  rely  on  a  suite  of  characterization  techniques  combined 
with  first-principles  computer  simulations. 

3.3.  Indirect  characterization  of  complexions 

When  Hart  predicted  that  grain  boundaries  may 
undergo  complexion  transitions,  he  suggested  that  the 
non-equilibrium  grain  boundary  properties,  such  as  grain 
boundary  mobility,  diffusivity,  sliding  resistance,  and  cohe¬ 
sive  strength,  would  be  most  dramatically  affected  [3-5],  He 
pointed  out  that  these  properties  are  also  much  more  easily 
accessible  to  experiment  than  the  equilibrium  properties, 
such  as  grain  boundary  energy,  width  and  adsorbate  con¬ 
tent.  Early  studies  by  Aust  [53]  on  grain  boundary  mobility 
in  zone-refined  Pb  demonstrated  a  discontinuity  in  the  acti¬ 
vation  energy  of  grain  boundary  mobility  in  carefully  con¬ 
structed  bicrystals. 

While  bicrystal  studies  are  useful  as  model  systems,  most 
engineering  materials  are  polycrystals  with  a  wide  variety 
of  general  grain  boundaries  present.  With  the  use  of 
dual-beam  focused  ion  beam  (F1B)/SEM  to  prepare  site- 
specific  TEM  samples,  it  is  possible  to  correlate  grain 
boundary  mobility  directly  with  atomic  structure  and 
chemistry  of  the  grain  boundary.  The  most  well-known 
example  of  such  a  study  was  done  on  polycrystalline 
alumina,  where  boundaries  with  increasing  mobilities  have 
been  associated  with  complexions  that  have  increasing 
levels  of  disorder  [17],  The  appearance  of  large  grains  in 
annealed  alumina  has  been  used  as  an  indicator  of  a 
complexion  transition.  For  example,  the  fast-moving 
boundaries  around  the  largest  grains  in  Ca-doped  Y203 
have  been  associated  with  an  IGF  complexion  [33],  Using 
this  behavior  to  identify  complexion  transitions,  Dillon 
and  Harmer  [19]  found  that  the  probability  of  complexion 
transitions  increased  with  both  temperature  and  grain 


boundary  adsorbate  concentration  (which  increases  with 
grain  size  in  a  sample  of  constant  bulk  composition).  In 
this  way,  abnormal  grain  growth  in  alumina  was  shown 
by  Dillon  et  al.  [17]  to  result  from  the  coexistence  of  two 
or  more  grain  boundary  complexions  with  differing  mobil¬ 
ity.  Fig.  20  shows  schematically  how  a  thicker,  more  disor¬ 
dered,  higher  mobility  complexion  exists  at  the  grain 
boundaries  of  abnormal  grains  in  alumina,  whereas  thin¬ 
ner,  more  ordered  complexions  exist  at  the  grain  bound¬ 
aries  of  normal  grains. 

The  changes  in  grain  boundary  mobility  can  be  mea¬ 
sured  to  characterize  the  grain  boundary  complexion  tran¬ 
sition.  The  mobility  can  be  measured  in  a  bimodal 
structure  by  characterizing  the  grain  growth  of  both  the 
normal  and  the  unimpinged  abnormal  grains  as  a  function 
of  time  [20,21],  The  size  of  the  normal  grains  is  used  to 
determine  the  driving  force  for  the  growth  of  the  abnormal 
grains.  If  the  grain  size  as  a  function  of  time  obeys  para¬ 
bolic  kinetics,  then  it  is  possible  to  determine  the  product 
of  the  grain  boundary  mobility  and  the  grain  boundary 
energy.  The  grain  boundary  energy  term  is  assumed  to  be 
relatively  constant  compared  to  the  variations  in  the  mobil¬ 
ity.  Mobility  measurements  made  using  this  technique, 
reported  by  Dillon  and  Harmer,  were  used  to  identify  six 
different  mobility  classes  for  grain  boundaries  in  alumina, 
as  shown  in  Fig.  21  [17].  However,  grain  growth  rates 
can  be  deceptive  as  they  are  sensitive  to  impurity  drag,  par¬ 
ticle  drag,  pore  drag  and  variations  in  driving  force  due  to 
stored  elastic  energy. 

Certain  complexion  transitions  can  also  be  detected  by 
differential  scanning  calorimetry  (DSC).  For  example,  in 
Al-Zn  [234]  and  Al-Mg  [235]  alloys,  it  has  been  shown  that 
the  grain  boundaries  and  triple  lines  melt  before  the  bulk  sol¬ 
idus  temperature  is  reached.  The  grain  boundary  premelting 
is  detected  by  the  DSC  as  an  absorption  of  heat  below  the 
solidus  temperature.  The  transition  is  only  detected  in  fine- 


monolayer 


trilayer 


Fig.  20.  A  schematic  diagram  showing  that  abnormal  grains  are 
surrounded  by  thicker,  more  disordered  complexions  (e.g.  the  trilayer 
complexion)  with  higher  mobility  (red),  whereas  the  normal  grains  have 
thinner,  more  ordered  complexions  (e.g.  the  monolayer  complexion)  with 
lower  mobility  grain  boundaries  (blue).  In  general,  it  is  believed  that  AGG 
results  from  such  coexistence  of  two  or  more  grain  boundary  complexions 
with  vastly  different  grain  boundary  mobilities.  Adapted  from  Ref.  [17], 
(For  interpretation  of  the  references  to  color  in  this  figure  legend,  the 
reader  is  referred  to  the  web  version  of  this  article.) 
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Fig.  21.  Grain  boundary  mobility  in  doped  and  undoped  alumina,  showing  six  distinct  regimes  of  mobility  corresponding  to  the  six  discrete  Dillon- 
Harmer  complexions  which  are  shown  schematically.  Adapted  from  Ref.  [17]. 


grained  samples  (10-20  micron  grain  size,  in  this  case)  when 
there  is  enough  grain  boundary  area  to  create  a  detectable 
heat  signature  [235].  Samples  with  large  grains,  produced 
in  situ  by  melting  and  solidifying,  do  not  have  a  detectable 
transition.  The  result  on  a  phase  diagram  is  that  a  separate 
grain  boundary  solidus  appears  at  a  temperature  below  the 
bulk  solidus,  e.g.  as  shown  by  the  solid  blue  premelting/ 
prewetting  ( Tpw)  line  in  the  solid  two-phase  region  in  Fig.  9 
[236].  In  principle,  complexion  transitions  other  than  melt¬ 
ing  should  be  detectable  by  DSC.  However,  small  enthalpy 
changes  associated  with  the  transition  and  the  simultaneous 
occurrence  of  competitive  enthalpic  changes,  such  as  grain 
growth,  may  mask  the  signal. 

In  several  cases,  it  has  been  shown  that  complexion  tran¬ 
sitions  are  associated  with  reductions  in  the  grain  boundary 
energy  [28,33,42],  Therefore,  grain  boundary  energy  mea¬ 
surements  can  be  used  to  characterize  complexion  transi¬ 
tions.  In  practice,  it  is  relatively  simple  to  measure 
changes  in  the  relative  grain  boundary  energy,  but  much 
more  difficult  to  measure  absolute  energies  [237],  As  a 
result,  it  is  usually  the  relative  energy  that  is  measured. 
Thermal  grooves  form  to  balance  the  interfacial  energies 
between  the  surfaces  and  the  grain  boundary.  The  relative 
energy  of  a  grain  boundary  ( ygh )  to  the  adjacent  grain  sur¬ 
face  ('/,)  can  be  expressed  as  a  function  of  the  dihedral 
angle  (¥',)  at  the  groove  according  to  Eq.  (16). 

—  =  2cos  —  (16) 

A  2 

Using  Mullins’s  [238]  analysis,  it  is  possible  to  measure 
the  height  and  width  of  a  thermal  groove  and  solve  for 
the  relative  interfacial  energy  of  the  grain  boundary.  How¬ 
ever,  this  method  includes  a  number  of  approximations 
and  assumptions.  For  example,  it  is  assumed  that  the  two 
surface  energies  are  the  same,  that  the  grain  boundary  is 
normal  to  the  surface  and  that  the  interface  energy  anisot¬ 
ropy  is  small.  While  these  assumptions  will  not  hold  for 
any  single-grain  boundary,  it  has  been  found  that,  for 
many  measurements  of  grain  boundary  dihedral  angles, 


the  mean  value  and  width  of  the  distribution  are  reproduc¬ 
ible  and  meaningful  quantities  [239-241].  In  the  last  dec¬ 
ade,  atomic  force  microscopy  (AFM)  has  greatly 
simplified  these  measurements  and  specific  procedures  have 
been  established  for  making  reliable  grain  boundary  energy 
measurements  using  AFM  [241]. 

As  an  example,  the  data  in  Fig.  22  shows  the  relative 
grain  boundary  energy  distributions  for  two  subsets  of 
grain  boundaries  in  the  same  100  ppm  Nd-doped  alumina 
[30],  These  measurements  were  derived  from  thermal 
grooves  according  to  Eq.  (16).  In  this  case,  the  crystallo¬ 
graphic  characteristics  of  the  junctions  were  not  considered 
and  we  are  concerned  only  with  the  distribution  derived  by 
averaging  over  more  than  200  boundaries.  The  first  subset 
of  boundaries  surrounds  small  grains  growing  normally 
(NGG)  and  the  second  subset  surrounds  grains  growing 
abnormally  (AGG).  The  grain  boundaries  growing  nor¬ 
mally  are  known,  on  average,  to  have  a  monolayer  com¬ 
plexion  with  adsorbed  Nd,  whereas  the  grains  growing 
abnormally  have,  on  average,  a  bilayer  complexion  with 
adsorbed  Nd  and  greater  mobility.  The  results  clearly  show 
that  the  boundaries  enriched  with  Nd  have  a  lower  average 
energy  than  those  with  less  Nd.  Similar  changes  have  been 
observed  in  other  systems  [28,33,30],  Similarly,  a  study  of 
IGFs  that  formed  at  gold-sapphire  interfaces  in  contact 
with  anorthite  glass  has  shown  that  the  IGF  reduces  the 
interfacial  energy  by  about  1 90  m.l  nT  2  compared  to  the 
intrinsic  gold-sapphire  interfacial  energy  [42], 

The  changes  in  the  grain  boundary  energy  have  an  indi¬ 
rect  influence  on  the  distribution  of  grain  boundaries  in  the 
sample.  From  comparative  studies  of  the  anisotropy  of  the 
grain  boundary  character  distribution  (GBCD,  the  relative 
areas  of  different  types  of  grain  boundaries)  and  the  anisot¬ 
ropy  of  the  grain  boundary  energy,  it  has  been  shown  that 
there  is  an  inverse  relationship  between  the  two  quantities 
in  polycrystals  evolving  by  normal  grain  growth,  which 
implies  that  there  are  relatively  more  low  energy  grain 
boundaries  than  there  are  high  energy  grain  boundaries 
[237,242-245],  Therefore,  if  there  is  a  transition  in  the  kinds 
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Fig.  22.  Cumulative  distribution  of  dihedral  angles  in  neodymia-doped 
alumina  annealed  at  1400  °C  with  normal  grain  boundaries  (monolayer 
complexion)  and  abnormal  grain  boundaries  (bilayer  complexion).  Insets 
schematically  illustrate  the  boundary  structure  of  the  two  complexions. 
Reprinted  from  Ref.  [30]  with  permission. 

of  complexions  in  a  sample  and  this  alters  the  grain  bound¬ 
ary  energies,  it  will  also  change  the  GBCD.  This  has  been 
observed  and  measured  in  several  cases  [29,33],  For  exam¬ 
ple,  Fig.  23(a  and  b)  compares  the  microstructure  of  a  Ca- 
doped  yttria  sample  before  and  after  some  of  the  bound¬ 
aries  have  undergone  a  complexion  transition.  The  abnor¬ 
mally  large  grains  have  an  IGF  complexion  and  the  smaller 
ones  have  a  bilayer  complexion. 

When  the  distributions  of  grain  boundary  planes  are 
compared  (Fig.  23(c  and  d)),  the  distribution  in  the  sample 
with  a  normal  grain  size  distribution  is  nearly  isotropic, 
with  a  mild  preference  for  grain  boundaries  with  a  ( 1 1 1 ) 
orientation,  while  in  the  sample  after  the  complexion  tran¬ 
sition  there  is  a  strong  preference  for  grain  boundaries  with 
(100)  planes.  This  suggests  that  the  energies  of  the  (10  0) 
planes  have  been  selectively  reduced  in  energy  compared 
to  the  other  types  of  grain  boundaries.  Similar  effects  have 
been  observed  in  doped  aluminas  [29], 

4.  Terminology  and  categorization  of  complexions 

A  wide  variety  of  terminology  has  been  used  to  describe 
grain  boundary  complexions  and  transitions  between  dif¬ 
ferent  complexions.  The  terminology  can  be  confusing, 
and  sometimes  several  terms  are  used  in  the  literature  to 
describe  the  same  phenomenon.  Furthermore,  there  is  no 
generally  accepted  method  of  categorizing  the  many  differ¬ 
ent  types  of  complexion  transitions  that  may  occur.  Now 
that  it  is  possible  to  see  more  details  of  the  structure  and 
chemistry  of  grain  boundary  complexions  using  advanced 
electron  microscopy,  it  is  becoming  increasingly  important 
to  develop  a  standardized  nomenclature  for  describing  and 
categorizing  complexions.  A  categorization  system  will  be 


especially  important  as  grain  boundary  complexion 
diagrams  are  developed,  so  that  lines  of  transition  can  be 
marked  appropriately. 

In  this  section,  we  will  first  give  some  examples  of  the 
various  terms  that  have  been  used  to  describe  grain  bound¬ 
ary  complexions  and  transitions  between  different  complex¬ 
ions.  Then  we  will  discuss  a  categorization  scheme  for 
complexion  transitions  first  proposed  by  Cahn  [6],  Finally, 
we  will  suggest  a  method  to  categorize  complexions  based 
upon  a  combination  of  their  physical  characteristics  and 
properties. 

4.1.  Common  terminology  used  to  describe  grain  boundary 
complexion  transitions 

Many  descriptive  terms  have  been  used  when  discussing 
grain  boundary  complexion  transitions,  e.g.  structural 
transitions,  faceting  transitions,  defaceting  transitions, 
roughening  transitions,  wetting  transitions,  prewetting 
transitions  and  premelting  transitions.  Countless  variations 
of  these  terms  exist  in  the  literature,  yet  rigorous  definitions 
are  rarely  given  and  therefore  the  terms  can  be  confusing  to 
non-experts.  Complexion  terminology  that  invokes  the 
word  “wetting”  can  be  particularly  misleading,  because  a 
wetting  phase  is  a  bulk  phase  (with  arbitrary  thickness) 
rather  than  an  interfacial  state  (with  a  fixed  equilibrium 
thickness).  We  will  review  some  of  these  terms  here. 

As  discussed  in  Section  2,  a  first-order  complexion  tran¬ 
sition  occurs  when  there  is  a  discontinuity  in  the  first  deriv¬ 
ative  of  y  as  a  function  of  one  of  the  thermodynamic 
parameters  ( T.  P.  pu  R,  or  it),  or  in  a  higher-order  deriva¬ 
tive  for  a  second  or  higher-order  transition  [8,9].  Thus, 
cusps  in  grain  boundary  energy  as  a  function  of  grain 
boundary  misorientation  R  or  inclination  h  represent  com¬ 
plexion  transitions,  as  do  the  disappearance  of  the  cusps  as 
a  bulk  thermodynamic  parameter  is  varied  (e.g.  as  temper¬ 
ature  is  raised)  [8], 

A  structural  transition  is  defined  as  a  complexion  tran¬ 
sition  that  occurs  when  a  bulk  thermodynamic  parameter 
is  varied  ( T,  P,  etc.)  while  all  five  interfacial  thermody¬ 
namic  parameters  (it,  R)  are  held  constant  [8],  Despite  the 
name,  structural  transitions  can  involve  changes  in  the 
chemical  composition  of  the  grain  boundary  in  addition 
to  structural  rearrangement  of  atoms  within  the  grain 
boundary  core.  A  structural  transition  can  be  first-order 
if  it  results  from  the  crossing  of  two  or  more  “competing” 
grain  boundary  free  energy  curves  (as  shown  schematically 
in  Fig.  2),  or  it  can  be  second-order,  for  example,  if  the 
transition  involves  the  gradual  replacement  of  one  type  of 
grain  boundary  repeating  structural  unit  for  another  [8], 
By  definition,  all  structural  transitions  are  congruent  tran¬ 
sitions  and  vice  versa  (congruent  transitions  are  discussed 
in  Section  4.2). 

A  faceting  transition  is  a  complexion  transition  that 
involves  the  restructuring  of  an  initially  singular  or  curved 
grain  boundary  into  a  hill-and-valley  structure  [246]  with  at 
least  two  different  boundary  inclinations,  /q  and  «2,  as  a 
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Fig.  23.  EBSD  data  for  a  Ca-doped  yttria  sample  with  superimposed  reconstructed  grain  boundaries  for  (a)  an  NGG  sample,  which  exhibits  a 
homogeneous  grain  size  distribution  before  some  boundaries  undergo  a  complexion  transition,  and  (b)  an  AGG  sample,  displaying  a  bimodal  grain  size 
distribution  (as  a  result  of  some  boundaries  undergoing  a  complexion  transition  to  a  higher  mobility  complexion)  and  an  apparent  change  in  texture.  The 
grain  boundary  plane  distributions,  independent  of  misorientation,  for  the  (c)  NGG  sample  and  (d)  AGG  sample  plotted  in  the  standard  stereographic 
triangle.  Reprinted  from  Ref.  [33]  with  permission. 


bulk  or  interfacial  thermodynamic  variable  is  varied  ( T,  P. 
(it,  or  R,  for  example)  [8].  The  reverse  transition,  from  fac¬ 
eted  to  planar  grain  boundary,  is  referred  to  as  a  defaceting 
transition.  Reversible  faceting/defaceting  transitions 
induced  by  changes  in  composition  have  been  demon¬ 
strated  in  the  Cu-Bi  system  [247].  Each  facet  represents  a 
complexion  and  is  associated  with  a  specific  free  energy. 
The  coexistence  of  two  such  facets  is  shown  in  Fig.  3(b). 
These  facets,  initially  present  on  a  fine  scale  as  microfacets, 
may  coarsen  into  macrofacets  to  eliminate  the  energy  asso¬ 
ciated  with  the  one-dimensional  complexion  boundary 
between  each  pair  of  facets.  Microfacets  may  only  be  visi¬ 
ble  by  transmission  electron  microscopy  (e.g.  as  shown  in 
Fig.  24(a))  and  therefore  this  type  of  complexion  transition 
can  be  easily  overlooked.  By  definition,  all  faceting  transi¬ 
tions  are  non-congruent  transitions  and  vice  versa  (non- 
congruent  transitions  are  discussed  in  Section  4.2). 

Recent  experimental  evidence  in  Bi-doped  Cu  has  dem¬ 
onstrated  that  faceting  transitions  may  be  associated  with 
changes  in  the  structure  and  chemistry  of  the  grain  bound¬ 
ary  core  [35],  In  this  study,  a  Bi-rich  bilayer  complexion 
was  observed  on  one  of  the  two  facet  inclinations  (see 
Fig.  24).  However,  it  is  not  known  if  the  grain  boundary 
faceting  transition  occurred  simultaneously  with  formation 
of  the  Bi-rich  bilayer  complexion.  It  is  possible  that  the 


Bi-rich  bilayer  complexion  formed  before  grain  boundary 
faceting,  in  which  case  the  bilayer  complexion  formation 
would  be  a  structural  transition,  which  was  then  followed 
by  a  faceting  transition.  We  can  simply  say  that  a  complex¬ 
ion  transition  has  occurred  and  then  describe  its  character¬ 
istic  features,  rather  than  attempt  to  categorize  the 
transition  as  a  faceting  transition  or  structural  transition, 
which  could  be  misleading  because  faceting  and  structural 
transitions  are  mutually  exclusive  by  definition. 

A  grain  boundary  roughening  transition  is  a  complexion 
transition  that  occurs  when  an  initially  smooth  grain 
boundary  becomes  disordered  and  rough  on  an  atomic 
scale.  Roughening  transitions  are  well  known  and  have 
been  extensively  studied  [150,248-251],  They  occur  at  a 
critical  temperature  (which  may  vary  as  a  function  of  it 
and  R).  The  high-temperature  rough  grain  boundary  is  pre¬ 
dicted  to  have  a  significantly  larger  mobility  than  the  low- 
temperature  smooth  grain  boundary  [150],  Different  grain 
boundaries  roughen  at  different  temperatures  and  not  all 
grain  boundaries  may  undergo  a  roughening  transition, 
even  at  temperatures  as  high  as  the  melting  temperature. 
Computer  simulations  predict  that  grain  boundaries  which 
remain  smooth  at  high  homologous  temperatures,  and 
hence  are  low  mobility  boundaries,  are  responsible  for 
grain  growth  stagnation  in  pure  materials  [248];  such 
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Fig.  24.  HAADF-STEM  images  of  faceted  grain  boundary  in  Bi-doped  Cu,  exhibiting  the  coexistence  of  two  complexions  on  different  inclinations  of  a 
faceted  grain  boundary.  There  is  a  Bi-rich  bilayer  complexion  on  one  of  the  two  facet  inclinations  and  a  clean  complexion  on  the  other  facet.  This  type  of 
complexion  coexistence  is  represented  schematically  in  Fig.  3(b)  and  Fig.  26.  Reprinted  from  Ref.  [35]  with  permission. 


stagnation  at  high  temperatures  cannot  otherwise  be 
explained  by  solute  drag  or  other  mechanisms. 

Grain  boundary  roughening  has  long  been  presumed  to 
be  a  second-order  transition,  but  a  recent  simulation  sug¬ 
gests  that  roughening  transitions  can  also  be  first-order 
[150],  Although  roughening  and  faceting  are  distinct  phe¬ 
nomena  [248],  it  is  possible  to  have  combined  roughen- 
ing/defaceting  transitions  in  which  an  initially  faceted 
boundary  (with  two  smooth  facets)  decomposes  upon  heat¬ 
ing  into  an  undulating,  rough  boundary  [1,94].  It  is  also 
possible  that  one  or  both  of  the  two  orientations  in  a 
hill-and-valley  structure  can  be  rough  rather  than  smoothly 
faceted  [251], 

Wetting,  prewetting  and  premelting  transitions  are  clo¬ 
sely  related  phenomena.  The  prototypical  wetting  transi¬ 
tion  involves  three  bulk  phases  -  a  liquid  drop,  /,  a  solid 
substrate,  s,  and  a  vapor  phase,  v,  as  shown  schematically 
in  Fig.  25(a).  Equilibrium  between  the  three  phases  can 
be  described  by  a  balance  of  three  interfacial  energies,  ysh 
ysv  and  y&.  The  contact  angle,  cj),  defines  the  degree  of  wet¬ 
ting.  Poor  wetting  occurs  if  (p  >  90°,  good  wetting  if  <p 
<  90°  and  non-wetting  if  cp  =  1 80°  [252],  Complete  wetting 
occurs  if  cp  =  0°,  as  shown  in  Fig.  25(b).  A  wetting  transi¬ 
tion  occurs  when  the  interface  energies  change  such  that 
the  wetting  phase  with  an  initial  wetting  angle  cp  >  0°  flat¬ 
tens  out  and  covers  the  entire  surface  of  the  substrate  with 
a  wetting  angle  of  <p  =  0°.  Wetting  transitions  need  not 
involve  a  liquid  phase,  and  can  occur  even  if  all  three  of 
the  phases  are  solid.  This  situation  is  shown  schematically 
in  Fig.  25(c  and  d),  which  depicts  a  solid  8  phase  existing  at 
a  8—8  grain  boundary  with  a  finite  contact  angle  (Fig.  25(c)) 
and  then  undergoing  a  wetting  transition  such  that  the  8 
phase  covers  the  entire  grain  boundary,  as  shown  in 
Fig.  25(d).  Although  it  can  be  defined  in  terms  of  bulk 
phases,  this  solid-state  wetting  transition  can  also  be 
viewed  as  a  dissociation  complexion  transition  in  which 


the  initial  s-s  grain  boundary  complexion  dissociates  into 
two  different  8-8  grain  boundary  complexions.  A  dissocia¬ 
tion  transition  is  shown  schematically  in  Fig.  3(c),  and 
occurs  when  the  grain  boundary  inclination  it  remains  con¬ 
stant  while  the  misorientation  R  changes  as  a  single  grain 
boundary  complexion  decomposes  into  two  complexions. 

Some  authors  use  the  convention  that  the  wetting  phase 
described  above  is  non-wetting  if  cp  >  90°  [252],  Confus¬ 
ingly,  the  phrase  “non-wetting  film”  has  been  used  to  refer 
to  nanolayer  complexions  (historically  referred  to  as  IGFs) 
[253],  The  rationale  behind  this  terminology  is  that  a  wet¬ 
ting  film  (which  may  have  an  arbitrary  thickness)  is  a  bulk 
phase  because  wetting  is  a  bulk  phenomenon,  and  therefore 
a  non-wetting  film  (which  has  a  fixed,  finite  thickness)  is  a 
non-bulk  phase  and  is  thus  an  interfacial  phase  -  what  we 
call  a  complexion.  This  terminology  is  especially  confusing 
in  light  of  the  convention  that  a  non-wetting  phase  is 
defined  based  on  a  non-zero  contact  angle  as  discussed 
above  (</>  >  90°  or  cp  =  180°,  depending  on  the  convention). 
The  phrase  “non-wetting  film”  therefore  implies  a  non-zero 
contact  angle  and  suggests  that  the  complexion  does  not 
cover  the  entire  grain  boundary,  which  is  misleading. 
Rather,  it  is  understood  that  a  complexion  completely  cov¬ 
ers  the  entire  interface  at  the  grain  boundary,  because  it  is 
the  interface,  if  two  complexions  coexist  at  equilibrium  (as 
shown  in  Fig.  3(a),  for  example),  then  each  one  does  not 
cover  the  entire  interface  but  there  is  no  notion  of  a  contact 
angle  between  them.  The  complexion  boundary  that  sepa¬ 
rates  these  two  complexions  is  essentially  one-dimensional, 
so  the  notion  of  a  contact  angle  implied  by  the  term  “non¬ 
wetting”  does  not  apply  to  complexions. 

Prewetting  and  premelting  transitions  are  both  complex¬ 
ion  transitions  whose  names  allude  to  the  related  bulk  phe¬ 
nomena  (wetting  and  melting).  Prewetting  describes  a  first- 
order  complexion  transition  in  which  the  interfacial  thick¬ 
ness  and  the  adsorbed  solute  increase  discontinuously. 
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Fig.  25.  Schematic  wetting  diagrams:  (a)  the  prototypical  three-phase  wetting  geometry  when  a  liquid  drop  is  on  a  solid  substrate  in  contact  with  a  vapor; 
(b)  the  geometry  after  a  wetting  transition,  when  the  contact  angle  goes  to  zero;  (c,  d)  analogous  solid-solid  wetting  phenomena  at  a  grain  boundary  in  a 
two-phase  solid,  with  the  5  phase  (c)  partially  wetting  the  e-e  grain  boundary  and  (d)  completely  wetting  the  s-e  grain  boundary. 


Prewetting  is  often  described  as  a  “thin  to  thick  film  tran¬ 
sition”  [23].  Both  the  “thin”  and  the  “thick”  film  are  com¬ 
plexions  rather  than  bulk  films.  Prewetting  transitions  are 
adsorption  transitions,  and  thus  typically  occur  in  multi- 
component  systems,  but  they  may  occur  in  either  one-  or 
two-phase  regions  on  the  bulk  phase  diagram,  as  shown 
in  Fig.  9.  The  prewetting  terminology  derives  from  Cahn’s 
critical  point  wetting  theory  for  fluids  [69],  Prewetting  does 
not  necessarily  imply  any  type  of  atomic  rearrangement  or 
restructuring  in  the  grain  boundary  core,  but  rather  simply 
refers  to  an  abrupt  jump  in  adsorbate  content  of  the  grain 
boundary  and  an  associated  discrete  increase  in  the  thick¬ 
ness.  It  typically  occurs  under  thermodynamic  conditions 
near  a  bulk  wetting  transition,  when  the  phase  that  would 
be  involved  in  the  wetting  transition  is  not  yet  a  thermody¬ 
namically  stable  bulk  phase. 

On  the  other  hand,  premelting  is  a  complexion  transi¬ 
tion  in  which  the  grain  boundary  inclination  remains 
invariant  and  the  core  becomes  disordered  and  liquid-like 
at  high  temperatures  near  the  melting  temperature  (in  con¬ 
trast,  roughening  transitions  result  in  local  changes  in  grain 
boundary  inclination  and  may  occur  at  much  lower  homol¬ 
ogous  temperatures).  Premelting  can  be  considered  as  a 
special  case  of  prewetting.  For  grain  boundaries  in  binary 
or  multicomponent  systems,  the  term  “premelting”  is  used 
to  emphasize  the  structural  disorder,  while  the  term  “prew¬ 
etting”  is  used  to  emphasize  the  compositional  variation 
(particularly  adsorption).  Complexion  transitions  tend  to 
have  elements  of  both  prewetting  (adsorption  transitions) 
and  premelting  (disordering  transitions),  which  has  led 
some  researchers  to  use  the  phrase  “coupled  prewetting/ 
premelting  transition”  to  refer  to  a  complexion  transition 
which  has  discontinuities  in  both  adsorption  and  structure 
[13]. 

4.2.  Congruent,  non-congruent  and  dissociation  transitions 

Cahn  suggested  classifying  grain  boundary  complexion 
transitions  into  two  types,  congruent  and  non-congruent 
[1,6],  During  a  congruent  transition,  the  structure  and/or 


chemistry  of  the  grain  boundary  core  changes  while  the 
morphology  of  the  grain  boundary  remains  invariant,  so 
that  flat  boundaries  remain  flat  and  curved  boundaries 
remain  curved.  Non-congruent  transitions  involve  different 
grain  boundary  inclinations,  it,  and  occur  when  a  grain 
boundary  facets  to  a  new  orientation  or  otherwise  under¬ 
goes  a  change  in  h.  A  third  type  of  transition  described 
by  Cahn  is  the  dissociation  transition,  in  which  one  grain 
boundary  separates  into  two  boundaries  by  the  insertion 
of  a  new  bulk  phase  [6],  The  two  new  boundaries  have  a 
total  free  energy  less  than  the  original  boundary,  and  the 
combined  misorientation  R  is  conserved.  This  transition 
is  commonly  referred  to  as  a  wetting  transition.  Cahn  con¬ 
cluded  that  congruent  complexion  transitions  are  exceed¬ 
ingly  rare,  saying  “The  search  for  congruent 
transformations  has  all  the  aspects  of  a  search  for  a  needle 
in  a  6  +  C  dimensional  haystack”  [6],  This  was  a  rather 
remarkable  conclusion,  because  most  reports  of  grain 
boundary  transitions  implicitly  assume  that  a  congruent 
transition  has  occurred. 

The  coexistence  of  two  complexions  at  equilibrium 
within  the  three  categories  of  transitions  discussed  by  Cahn 
can  be  visualized  according  to  Fig.  3:  (a)  congruent  transi¬ 
tion,  (b)  non-congruent  transition  and  (c)  dissociation  tran¬ 
sition.  As  discussed  by  Rottman  [9],  coarsening  of  the 
complexions  in  Fig.  3(a)  and  (b)  may  occur  to  minimize 
the  energy  associated  with  the  one-dimensional  y.-jl  com¬ 
plexion  boundaries.  If  the  two  facets  in  Fig.  3(b)  have  dif¬ 
ferent  specific  free  energies,  there  would  be  an  additional 
driving  force  for  the  elimination  of  the  higher-energy  facet 
and  the  growth  of  the  lower-energy  facet.  If  this  coarsening 
and  growth  occurred  quickly  enough,  evidence  of  a  non- 
congruent  (faceting)  transition  could  essentially  disappear 
for  a  grain  boundary  in  a  polycrystal  and  the  complexion 
transition  might  appear  to  have  been  congruent  because 
the  boundary  is  not  faceted.  This  possibility  complicates 
experimental  analysis:  if  a  grain  boundary  is  examined  ex 
situ  in  a  scanning/transmission  electron  microscope  and 
there  is  evidence  of  a  grain  boundary  complexion  transition 
(e.g.  multilayer  adsorption)  but  no  observable  faceting,  can 
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it  be  concluded  that  a  congruent  transition  has  occurred 
rather  than  a  non-congruent  transition?  In  light  of  Cahn’s 
comments  on  the  rarity  of  congruent  transitions,  it  seems 
presumptuous  to  classify  complexion  transitions  as  congru¬ 
ent  without  strong  evidence  that  the  transition  did  indeed 
occur  congruently,  without  a  change  in  grain  boundary 
inclination.  Therefore,  while  categorization  based  on  the 
scheme  of  congruent  vs.  non-congruent  transitions  is 
attractive  in  principle,  in  practice  it  is  rather  impractical 
because  it  is  not  always  apparent  which  type  of  transition 
has  occurred  by  studying  the  grain  boundary  ex  situ.  We 
therefore  recommend  a  more  empirical  categorization 
scheme  that  is  based  on  linking  the  structure  and  chemistry 
of  complexions  to  their  properties. 

4.3.  Other  methods  for  categorizing  complexions 

An  ideal  grain  boundary  complexion  categorization 
scheme  would  be  simple  enough  to  be  used  as  an  engineer¬ 
ing  tool  but  robust  enough  to  be  applicable  to  all  types  of 
complexions  in  all  materials  systems.  Unfortunately,  any 
categorization  scheme  will  necessarily  be  a  simplification 
that  does  not  capture  every  detail  and  nuance  of  every  pos¬ 
sible  complexion  or  transition.  A  method  of  categorizing 
complexion  transitions  can  be  very  useful,  especially  with 
regard  to  developing  practical  complexion  transition  dia¬ 
grams.  Despite  the  inherent  shortcomings  of  categoriza¬ 
tion,  we  can  make  some  generalizations  regarding  various 
types  of  grain  boundary  complexions,  and  it  is  possible  in 
some  instances  to  develop  a  categorization  scheme  based 
on  structure-properties  relationships  within  a  given  materi¬ 
als  system. 

Grain  boundary  complexion  transitions  can  be  catego¬ 
rized  broadly  into  two  groups,  intrinsic  and  extrinsic. 
Intrinsic  transitions  occur  in  pure  materials,  including 
unary  systems  and  stoichiometric  or  non-stoichiometric 
compounds,  and  may  include  structural  and  roughening 
transitions.  Extrinsic  transitions  occur  in  non-pure  materi¬ 
als  and  usually  involve  the  adsorption  of  impurity  or  dop¬ 
ant  elements.  Due  to  the  ubiquitous  presence  of  impurities 
in  real  materials,  intrinsic  complexion  transitions  are 
extraordinarily  difficult  to  study  experimentally.  The  mini¬ 
mum  concentration  of  impurities  that  can  induce  an  extrin¬ 
sic  complexion  transition  varies  depending  on  the  materials 
system,  but  as  little  as  30  ppm  CaO  has  been  shown  to 
induce  extrinsic  complexion  transitions  in  A  ECU,  for  exam¬ 
ple  [17]. 

One  of  the  earliest  studies  on  grain  boundary  complex¬ 
ion  transitions  identified  a  discontinuity  in  the  activation 
energy  of  the  grain  boundary  migration  rate  in  high-purity, 
zone-refined  Pb  [53].  The  study  was  done  on  Pb  bicrystals, 
and  the  discontinuity  was  assumed  to  be  evidence  of  an 
intrinsic  transition  between  a  low-temperature  and  high- 
temperature  grain  boundary  structure.  However,  impuri¬ 
ties  may  have  contributed  to  this  transition,  i.e.  it  may 
not  have  been  an  intrinsic  transition.  Modern  computer 
simulations  can  offer  more  insight  into  intrinsic  transitions. 


For  example,  as  discussed  earlier,  a  recent  MD  simulation 
of  X 5  grain  boundaries  in  pure  copper  demonstrated  a  first- 
order  grain  boundary  structural  transition  at  elevated  tem¬ 
peratures,  as  well  as  the  coexistence  of  two  complexions  on 
a  single  grain  boundary  separated  by  a  one-dimensional 
complexion  boundary  [39], 

Intrinsic  grain  boundary  roughening  transitions  have 
also  been  predicted  in  pure  metals  by  computer  simula¬ 
tions,  and  these  transitions  are  associated  with  abrupt  dis¬ 
continuities  in  grain  boundary  mobility  [150],  Some  grain 
boundaries  do  not  roughen  even  at  high  temperatures  near 
the  melting  point,  and  this  is  believed  to  be  responsible  for 
grain  growth  stagnation  [248],  Such  roughening  transitions 
would  be  difficult  or  impossible  to  study  by  standard  room- 
temperature  TEM,  which  is  the  predominant  mode  of  grain 
boundary  complexion  analysis  today.  In-situ  heating 
experiments  in  the  transmission  electron  microscope  could 
conceivably  demonstrate  defaceting/roughening  transi¬ 
tions,  and  indeed  this  has  been  done  for  multicomponent 
systems  previously  [247],  However,  such  experiments 
invariably  contain  impurities,  which  is  not  the  case  for 
computer  simulations,  and  therefore  would  not  be  very 
useful  for  studying  intrinsic  complexion  transitions.  Impu¬ 
rities  will  typically  be  present  in  TEM  samples  prepared 
using  a  FIB  microscope,  due  to  Ga  implantation  from 
the  ion  beam. 

There  has  not  been  enough  work  done  on  intrinsic  com¬ 
plexion  transitions  to  suggest  a  useful  categorization 
scheme.  Work  is  ongoing  in  this  area,  and  useful  advances 
are  likely  to  be  dominated  by  computer  simulations,  which 
can  be  impurity-free. 

Significantly  more  research  has  been  done  on  extrinsic 
complexion  transitions.  Since  most  engineering  materials 
are  multi-component  systems,  extrinsic  complexion  transi¬ 
tions  are  generally  of  greater  engineering  interest  than 
intrinsic  complexion  transitions.  A  useful  categorization 
scheme  for  extrinsic  grain  boundary  complexion  transitions 
will  likely  be  empirical,  system-specific  and  simple  enough 
so  that  complexion  transitions  can  be  overlaid  on  a  bulk 
phase  diagram,  but  not  so  simple  that  critical  information 
relevant  to  changes  in  materials  properties  is  omitted. 

To  develop  a  categorization  scheme,  the  structure  and 
chemistry  of  grain  boundaries  in  a  given  materials  system 
must  be  studied  experimentally,  typically  by  electron 
microscopy,  and  correlated  to  grain  boundary  properties 
of  interest.  Such  a  direct  correlation  between  grain  bound¬ 
ary  structure  and  properties  was  carried  out  in  doped 
A1203  by  Dillon  and  Harmer  [17],  TEM  samples  of  high- 
mobility  grain  boundaries  from  fast-growing  abnormal 
grains  and  normal  grain  boundaries  with  lower  mobility 
were  prepared  via  the  in  situ  FIB  lift-out  method  so  that 
grain  boundary  structure  and  chemistry  were  directly  cor¬ 
related  to  grain  boundary  mobility.  Six  distinct  mobility 
regimes  were  discovered,  as  shown  in  Fig.  21,  and  these 
regimes  correlate  to  six  categories  of  grain  boundary  com¬ 
plexions  shown  in  Fig.  10.  Of  course,  by  simplifying  a  vast, 
multi-dimensional  phase  space  into  only  six  categories, 
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some  information  is  necessarily  lost.  Nevertheless,  there  is 
strong  experimental  evidence  that  this  Dillon-Harmer  cat¬ 
egorization  scheme  works  very  well  for  understanding  how 
complexion  transitions  influence  grain  boundary  mobility 
in  doped  A1203  [17,20,21], 

All  six  complexion  categories  will  not  necessarily  exist 
in  other  materials  systems,  but  many  of  them  have  been 
shown  to  exist  in  several  other  materials.  Likewise,  six 
distinct  structures  are  not  sufficient  to  describe  the  types 
of  complexions  in  any  arbitrary  system.  In  a  system  such 
as  Bi-doped  Ni,  for  example,  only  two  major  types  of 
grain  boundary  complexions  were  identified  (a  clean  com¬ 
plexion  with  no  observable  Bi  adsorption,  and  a  bilayer 
complexion  rich  in  adsorbed  Bi).  These  two  types  of  com¬ 
plexions  correlate  strongly  to  grain  boundary  cohesive 
strength  [48].  Although  other  complexion  categorization 
schemes  might  make  sense  for  Ni-Bi  or  doped  A1203  if 
grain  boundary  properties  other  than  cohesive  strength 
and  mobility  are  studied,  these  two  examples  illustrate 
how  grain  boundary  structure  and  chemistry  can  be 
correlated  with  grain  boundary  properties  to  develop 
an  empirical,  system- specific  complexion  categorization 
scheme. 

Another  categorization  scheme  for  grain  boundary  com¬ 
plexions  has  been  previously  proposed  in  which  three  cate¬ 
gories  exist:  dry,  moist  and  wet,  as  defined  by  Cannon  et  al. 
[254—256],  Despite  the  terminology  that  implies  liquid 
phases  are  involved,  these  three  categories  usually  refer  to 
grain  boundaries  and  other  solid-solid  interfaces.  The  dry 
category  corresponds  to  the  clean  complexion  in  the 
Dillon-Harmer  categorization  scheme,  in  which  the  grain 
boundary  is  free  of  solute  adsorption  or  contains  a  very 
small  amount  of  solute.  The  moist  category  encompasses 
complexions  that  display  multilayer  adsorption,  which 
corresponds  to  the  Dillon-Harmer  bilayer,  trilayer  and 
nanolayer  complexions.  The  wet  category  describes  the 
case  where  a  bulk  wetting  film  exists  at  the  boundary, 
and  this  corresponds  directly  to  the  Dillon-Harmer  wetting 
category. 

The  terms  and  definitions  discussed  in  this  section  are 
summarized  in  Table  2.  This  table  is  designed  as  a  quick 
reference  guide  for  terms  that  have  been  used  to  describe 
complexions  and  complexion  transitions.  In  an  effort  to 
unify  terminology  in  grain  boundary  complexions  research, 
we  recommend  using  a  subset  of  these  terms  to  describe 
grain  boundary  complexions  and  transitions  between  dif¬ 
ferent  complexions.  The  recommended  terminology  is  sum¬ 
marized  in  Fig.  26:  complexion  transitions  are  shown 
schematically  in  the  top  section  of  the  figure  along  with 
the  preferred  terms,  while  methods  of  categorizing  com¬ 
plexions  are  shown  schematically  in  the  lower  section. 
While  the  terminology  in  Fig.  26  is  not  necessarily  sufficient 
to  describe  all  possible  complexion  phenomena,  it  provides 
a  framework  within  which  specific  complexion  behavior 
and  phenomena  may  be  discussed.  For  example,  the 
roughening  complexion  transition  is  not  explicitly  shown, 
but  there  could  be  a  “congruent,  intrinsic  roughening 


complexion  transition”  in  a  pure  material  in  which  the 
macroscopic  grain  boundary  geometry  (character)  does 
not  change  but  the  boundary  becomes  disordered  and 
rough  on  the  atomic  (microscopic)  scale  upon  heating  to 
a  critical  temperature.  Similarly,  one  can  imagine  other 
types  of  complexion  phenomena  that  can  be  described  in 
greater  detail  within  this  framework. 

Some  of  the  terms  are  specific  to  grain  boundary  com¬ 
plexions,  while  other  terms  can  be  used  for  any  type  of 
complexion,  whether  it  is  a  surface  complexion,  a  grain 
boundary  complexion  or  a  phase  boundary  complexion. 

5.  The  effect  of  complexions  on  properties,  processing  and 
microstructure 

Efforts  to  better  understand  anomalous  microstructural 
development  and  materials  properties  have  motivated 
much  of  the  development  of  the  science  of  grain  boundary 
complexions  [17,47,203,222,257],  Grain  boundaries  may 
undergo  complexion  transitions  independently  of  bulk 
phase  transformations,  and  thus  grain  boundary  related 
properties  may  change  unexpectedly  and  unpredictably  as 
temperature,  pressure  or  chemical  potential  is  varied.  A 
complexion  transition  can  affect  equilibrium  properties 
such  as  grain  boundary  energy,  entropy,  enthalpy,  defect 
density  and  adsorbate  concentration.  Complexion  transi¬ 
tions  may  also  induce  discontinuous  changes  in  the  non¬ 
equilibrium  properties  of  grain  boundaries,  such  as  mobil¬ 
ity,  cohesive  strength  and  sliding  resistance.  It  is  the  discon¬ 
tinuity  in  properties  that  is  the  hallmark  of  a  complexion 
transition,  and  is  the  reason  why  such  transitions  play  a 
large  role  in  the  properties,  processing  and  microstructure 
development  of  materials  [5],  The  effect  of  a  grain  bound¬ 
ary  transition  on  these  properties  is  not  easy  to  predict 
and  may  be  difficult  to  rationalize  on  the  basis  of  classical 
models  for  grain  boundary  structure-property  relation¬ 
ships.  For  example,  a  grain  boundary  complexion  transi¬ 
tion  in  Ga-doped  Al  increased  grain  boundary  mobility 
by  approximately  an  order  of  magnitude  [258],  in  stark 
contrast  to  the  classical  impurity  drag  theory,  which  pre¬ 
dicts  that  grain  boundary  segregation  will  always  reduce 
grain  boundary  mobility  [259], 

Despite  their  relatively  small  volume  fraction,  grain 
boundaries  often  play  a  dominant  role  in  materials  proper¬ 
ties.  Grain  boundary  engineering  was  first  proposed  by 
Watanabe  in  1984  [260],  and  has  since  been  considered 
by  others.  While  grain  boundary  engineering  recognizes 
the  effects  of  grain  boundary  solute  segregation,  the 
emphasis  is  usually  on  the  geometric  structure  of  grain 
boundaries  and  a  categorization  scheme  based  on  the  coin¬ 
cident  site  lattice  model  that  compares  “low-angle”  and 
“high-angle”  boundaries,  or  “special”  to  “general”  bound¬ 
aries.  Clearly,  grain  boundary  complexion  transitions 
deserve  attention  in  the  context  of  grain  boundary  engi¬ 
neering,  especially  given  the  non-classical  behavior  induced 
by  some  grain  boundary  complexion  transitions.  In  partic¬ 
ular,  grain  boundary  complexion  engineering  may  provide 
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new  opportunities  to  engineer  nano-  and  microstructures, 
control  mass  transport  dependent  properties  and  tailor 
mechanical  properties.  In  this  section  we  will  discuss  the 
effect  of  complexion  transitions  on  grain  boundary  proper¬ 
ties,  bulk  properties,  and  materials  processing. 

5.1.  Grain  boundary  properties 

During  a  complexion  transition,  the  slope  of  the  grain 
boundary  energy  (as  a  function  of  temperature,  pressure, 
chemical  potential,  etc.)  will  exhibit  special  behavior  such 
that  a  discontinuity  may  appear  in  the  first  derivative  (or 
a  higher  order  derivative)  of  the  grain  boundary  energy 
curve.  This  discontinuity  in  the  derivative  of  the  grain 
boundary  energy  is  the  defining  characteristic  of  a  com¬ 
plexion  transition.  However,  it  is  difficult  and  time  consum¬ 
ing  to  measure  the  grain  boundary  energy  curves. 
Therefore,  discontinuities  in  grain  boundary  properties, 
structure  and  chemistry  can  serve  as  excellent  proxies  for 
identifying  grain  boundary  complexion  transitions. 

One  early  such  study  measured  the  grain  boundary 
energy  in  Pb  by  observing  the  dihedral  angles  of  a  tricrystal 
as  a  function  of  temperature  [261],  Based  on  these  measure¬ 
ments,  which  showed  a  discontinuous  change  in  dihedral 
angle  at  a  specific  temperature,  it  was  claimed  that  a  grain 
boundary  transition  had  occurred  and  was  accompanied  by 
a  discontinuous  change  in  grain  boundary  energy.  Because 
grain  boundary  energy  must  be  continuous  even  at  transi¬ 
tion  points,  this  work  has  been  criticized  by  others 
[5,6,262].  Cahn  discussed  the  controversy  of  discontinuities 
in  dihedral  angles  and  pointed  out  that  Gleiter’s  work 
showed  evidence  of  a  grain  boundary  faceting  transition, 
which  is  itself  evidence  of  a  grain  boundary  complexion 
transition,  and  furthermore  that  a  discontinuity  in  the  dihe¬ 
dral  angle  (which  was  incorrectly  assumed  to  prove  a  dis¬ 
continuity  in  energy)  accompanying  such  a  faceting 
transition  is  not  inconsistent  with  thermodynamic  princi¬ 
ples  [6].  This  early  study  in  Pb  tricrystals  and  the  subse¬ 
quent  criticism  and  analysis  demonstrates  the  difficulty 
involved  in  properly  identifying  grain  boundary  complex¬ 
ion  transitions  by  experimentally  measuring  the  grain 
boundary  energy.  A  number  of  studies  have  observed  hys¬ 
teresis  in  complexion  transitions  upon  heating  and  cooling 
[28,108,112],  Such  hysteresis  results  from  a  complexion 
transition  occurring  under  non-equilibrium  conditions 
(e.g.  at  a  higher  or  lower  temperature)  and  can  lead  to 
the  measurement  of  a  discontinuous  change  in  energy. 

A  more  recent  study  of  grain  boundary  energy  in  Bi- 
doped  Cu  demonstrated  that  a  grain  boundary  complexion 
transition  can  be  identified  by  measuring  the  dihedral  angle 
of  a  grain  boundary  groove  using  AFM  [263],  The  dihedral 
angle  measurements  in  the  Cu  Bi  bicrystals  showed  two 
discontinuities  in  the  slope  of  the  grain  boundary  energy 
as  a  function  of  Bi  content.  Based  on  the  transition  from 
approximately  ~1  ML  of  Bi  adsorption  to  ~2  ML  of  Bi 
adsorption,  which  was  measured  with  AES,  it  was  con¬ 
cluded  that  the  second  slope  discontinuity  was  associated 


with  a  grain  boundary  complexion  transition,  whereas 
the  first  slope  discontinuity  was  an  artifact  arising  from  a 
surface  complexion  transition  that  changed  the  dihedral 
angle  of  the  grain  boundary  groove.  Grain  boundary 
energy  measurements  using  dihedral  angles  are  relative 
measurements  that  compare  surface  energy  to  grain  bound¬ 
ary  energy.  Therefore,  surface  complexion  transitions  can 
create  discontinuities  in  the  slope  of  the  (relative)  grain 
boundary  energy  graphs  derived  from  these  measurements. 
This  study  on  Cu-Bi  bicrystals  [263]  highlights  the  problem 
of  relying  exclusively  on  dihedral  angle  measurements  of 
grain  boundary  grooves  to  identify  grain  boundary  com¬ 
plexion  transitions.  Concurrent  measurements  of  at  least 
one  additional  grain  boundary  property  are  typically 
needed  to  confirm  the  results  of  such  studies. 

The  difficulty  in  experimentally  measuring  grain  bound¬ 
ary  energy  has  led  researchers  to  measure  other  properties 
that  are  strongly  correlated  with  grain  boundary  complexion 
transitions,  such  as  grain  boundary  mobility  and  diffusivity. 
Discontinuities  in  these  properties  are  indirect  indicators  of 
complexion  transitions  in  which  an  increase  in  mass 
transport  kinetics  accompanies  increases  in  the  structural 
or  chemical  width  of  the  grain  boundary  [17,22,47, 
106,203,214].  Grain  boundary  diffusivity  measurements  in 
polycrystalline  Cu-Bi  alloys  demonstrated  an  abrupt 
increase  in  diffusivity  by  two  orders  of  magnitude  at  a  Bi  con¬ 
centration  less  than  the  bulk  solidus  concentration  [162], 
This  Bi  concentration  was  in  agreement  with  the  concentra¬ 
tion  at  which  a  discontinuity  in  the  slope  of  grain  boundary 
energy  was  measured  in  a  previous  experiment  [263].  These 
results  provide  confirmation  of  Hart’s  [5]  prediction  that 
grain  boundary  complexion  transitions  can  induce  large  dis¬ 
continuities  in  non-equilibrium  grain  boundary  properties. 
Although  discontinuities  in  the  grain  boundary  structure 
and  chemical  width  in  the  Cu-Bi  alloys  were  not  confirmed 
directly  by  TEM  in  these  studies,  the  observed  discontinuity 
in  diffusivity  as  a  function  of  composition  is  likely  associated 
with  a  complexion  transition  that  has  been  identified  by 
STEM  in  another  study  [35], 

Grain  boundary  mobility  measurements  can  also  be 
used  to  identify  complexion  transitions.  To  carry  out  a 
mobility  measurement  on  a  single  grain  boundary,  a  special 
bicrystal  geometry  is  typically  used  to  control  the  driving 
force  throughout  the  experiment  [264],  This  type  of  bicrys¬ 
tal  measurement  was  used  to  show  that  small  (10  ppm) 
additions  of  Ga  to  Al  increased  the  grain  boundary  mobil¬ 
ity  by  about  an  order  of  magnitude  [258],  The  activation 
enthalpy  and  the  preexponential  factor  both  changed  with 
Ga  doping  as  well,  suggesting  a  change  in  the  mechanism 
of  grain  boundary  motion  due  to  a  grain  boundary  com¬ 
plexion  transition.  Grain  boundary  mobility  can  also  be 
deduced  using  grain  growth  studies  on  polycrystalline  sam¬ 
ples,  as  was  done  with  doped  A1203  [17],  In  this  study, 
grain  boundary  mobility  spanning  more  than  three  orders 
of  magnitude  was  directly  correlated  with  six  types  of  grain 
boundary  complexions  in  polycrystalline  A1203  by  S/TEM 
studies  of  individual  grain  boundaries  (Fig.  21). 
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Table  2 

Terminology  related  to  complexion  transitions  and  methods  of  categorizing 

Complexion  Defined  by  Congruent  transition 

transitions  geometry 

Non-congruent 

transition 

Structural  transition 


Faceting  transition 


Dissociation  transition 


Defined  by  structure  Premelting  transition 

and/or  composition 

Prewetting  transition 


Adsorption  transition 

Complexion  Defined  by  Intrinsic 

categories  composition 

Extrinsic 


Defined  by  thickness  Dry 

and  composition 

Moist 


Wet 

Dillon-Harmer  Clean  complexion 

complexions, 

defined  by  thickness 

and  composition 

Monolayer  complexion 
Bilayer  complexion 
Trilayer  complexion 
Nanolayer  complexion 

Wetting 

Defined  by  degree  Ordered 

of  periodicity 


complexions. 

A  complexion  transition  that  occurs  without  a  change  in  grain 
boundary  character  (ff  and  h  remain  invariant);  typically  involves 
changes  in  atomic  structure  and  composition  of  the  grain  boundary 
core 

A  complexion  transition  that  results  in  a  change  in  grain  boundary 
character  ( R  and/or  h  change) 

A  complexion  transition  that  occurs  when  a  bulk  thermodynamic 
parameter  is  varied  ( T ,  P ,  ph  etc.)  while  all  five  interfacial 
thermodynamic  parameters  (»,  R)  are  held  constant  [8] 

A  complexion  transition  in  which  a  single  complexion  decomposes  into 
two  complexions:  during  a  faceting  transition,  the  grain  boundary 
plane  normal  h  decomposes  into  h\  and  ii2.  the  area-weighted  average 
of  which  is  equivalent  to  n 

A  complexion  transition  in  which  a  single  complexion  decomposes  into 
two  complexions:  during  a  dissociation  transition,  a  single  grain 
boundary  dissociates  into  two  new  interfaces,  separated  by  a  new  bulk 
phase,  with  misorientation  R,  and  R2,  whose  combined  misorientation 
is  equal  to  the  original  misorientation  R.  Also  known  as  a  “wetting 
transition” 

The  formation  of  a  disordered,  liquid-like  film  on  a  crystalline  surface 
( or  at  a  grain  boundary  or  phase  boundary)  at  a  temperature  below  the 
melting  temperature  (or  solidus)  of  the  bulk  crystalline  phase 
Occurs  when  a  nanolayer  complexion  of  fixed  equilibrium  thickness 
forms  at  the  interface  in  the  thermodynamic  vicinity  of  a  wetting 
transition,  i.e.  near  the  temperature  or  composition  at  which  a  wetting 
transition  would  occur 

A  dramatic  change  (usually  a  first-order  transition)  in  the  composition 
of  an  interface  in  which  the  relative  amount  of  solute  increases  or 
decreases  significantly 

Any  complexion  that  exists  in  pure  systems  such  as  an  elemental  metal; 
the  composition  of  the  complexion  is  identical  to  the  bulk  composition 
(although  the  density  may  be  different) 

Any  complexion  that  exists  in  a  non-pure  system,  e.g.  in  a  system  that 
is  intentionally  doped  with  additional  elements  or  a  system  that 
contains  unintentional  impurities,  such  that  the  complexion 
composition  is  in  general  not  equal  to  the  bulk  composition 
A  complexion  with  no  adsorbed  solute  or  very  little  adsorption 
(corresponds  to  the  monolayer  complexion  in  the  Dillon-Harmer 
scheme) 

A  complexion  with  multilayer  solute  adsorption  (corresponds  to  the 
bilayer,  trilayer  and  nanolayer  complexions  in  the  Dillon-Harmer 
scheme) 

Refers  to  the  existence  of  a  bulk  wetting  film  (solid  or  liquid)  at  a 
boundary  (corresponds  to  the  wetting  category  in  the  Dillon-Harmer 
scheme) 

A  complexion  that  is  structurally  abrupt  with  a  core  thickness  that  is 
not  detectably  widened  by  solute  segregation.  Solute  segregation  is  not 
necessarily  entirely  absent,  but  is  minimal  or  is  not  observed  at  all,  and 
does  not  lead  to  an  increase  in  thickness  of  the  grain  boundary  core, 
phase  boundary  core  or  surface 

A  complexion  in  which  adsorbed  solute  is  observed  and  the  majority  of 
the  solute  is  confined  to  a  thickness  of  a  single  atomic  layer 
A  complexion  in  which  most  of  the  adsorbed  solute  occupies  a 
thickness  equal  to  two  atomic  layers 

A  complexion  in  which  most  of  the  adsorbed  solute  occupies  a 
thickness  equal  to  three  atomic  layers 

A  complexion  in  which  the  adsorbed  solute  occupies  a  thickness 
greater  than  three  atomic  layers,  but  which  is  still  finite,  fixed  and 
governed  by  equilibrium  thermodynamics.  Equivalent  to  IGF 
Refers  to  the  existence  of  a  bulk  wetting  film  (solid  or  liquid)  at  a  grain 
boundary  or  interface.  There  are  two  complexions  (one  on  each  side  of 
the  wetting  film) 

A  complexion  which  has  a  recognizable  degree  of  long-range  structural 
or  chemical  periodicity 


(continued  on  next  page) 
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Table  2  ( continued ) 


Disordered 

A  complexion  with  no  recognizable  long-range  periodicity  in  either 
structure  or  composition 

Related  terminology 

Monolayer 

The  term  “monolayer”  is  used  in  surface  science  to  quantify  the 
coverage  of  a  surface  by  atoms  or  molecules  arranged  in  a  two- 
dimensional  layer.  There  are  several  different  conventions  for  defining 
monolayer.  A  common  convention  is  that  a  monolayer  is  “a  number 
density  equal  to  that  of  the  atoms  in  a  single  atomic  layer  of  the 
substrate  material  parallel  to  the  surface”  [104],  This  particular 
definition  breaks  down  for  general  grain  boundaries,  which  have  two 
different  atomic  planes  or  layers  parallel  to  the  boundary  plane,  and 
thus  another  convention  must  be  used.  The  lone  word  “monolayer” 
should  not  be  confused  with  the  phrase  “monolayer  complexion”, 
which  is  one  of  the  Dillon-Harmer  complexions 

Segregation 

An  equilibrium  phenomenon  that  occurs  in  multicomponent  materials, 
causing  the  composition  of  grain  boundaries  to  differ  from  the  overall 
composition  at  equilibrium 

Adsorption 

The  term  “adsorption”  is  often  used  to  refer  to  specific  interfacial 
excess,  T,  which  appears  in  Gibbs  adsorption  isotherm: 
dy  =  —'2AiTidiii  [252],  “Adsorption”  is  also  used  descriptively  to  refer 
to  the  phenomenon  of  the  interfacial  composition  being  different  than 
the  bulk  composition.  Although  it  is  sometimes  used  interchangeably 
with  “segregation”  when  speaking  of  grain  boundaries,  “adsorption” 
was  originally  used  to  discuss  this  phenomenon  when  it  occurs  at 
surfaces 

Intergranular 

A  film  of  material  at  a  grain  or  phase  boundary  that  is  approximately 

film  (IGF) 

1-2  nm  thick;  IGFs  have  been  widely  observed  in  various  ceramics  that 
contain  impurities  such  as  SbN4  with  Si02  impurities,  ZnO  with  Bi2C>3 
impurities,  and  SrTiC>3  with  Ti02  impurities.  An  IGF  is  referred  to  as  a 
nanolayer  complexion  under  the  Dillon-Harmer  categorization 
scheme 

The  work  of  adhesion,  Wad,  of  a  grain  boundary  is  the 
reversible  work  necessary  to  convert  a  grain  boundary  into 
two  free  surfaces.  For  brittle  fracture,  it  may  be  described 
by  the  relationship 

Wad  =  2ys-ygb  (17) 

Flere,  ys  and  ygb  are  the  surface  and  grain  boundary  ener¬ 
gies,  respectively.  The  reduction  in  grain  boundary  energy 
associated  with  disordering  complexion  transitions  sug¬ 
gests  that  the  work  of  adhesion,  and  therefore  the  grain 
boundary  strength,  should  increase.  However,  in  many  sys¬ 
tems  studied  experimentally,  such  complexion  transitions 
have  been  associated  with  grain  boundary  embrittlement 
[48,265],  Local  configurational  considerations  could  out¬ 
weigh  classical  thermodynamic  considerations  in  such  a 
case.  For  example,  high  enthalpy  bonds  forming  with  the 
adjacent  lattice  could  primarily  contribute  to  lowering  the 
boundary  energy,  while  percolating  low  enthalpy  bonds 
in  the  interface  could  induce  embrittlement.  This  mecha¬ 
nism  has  been  proposed  by  Luo  et  al.  [48]  to  explain 
embrittlement  in  Ni-Bi.  Alternatively,  the  grain  boundary 
energy  may  be  reduced  by  an  increase  in  entropy.  An  entro¬ 
py  increase  may  result  from  an  increase  in  the  free  volume 
of  the  boundary,  point  defects  or  site  disorder,  for  example, 
which  could  also  decrease  the  strength  of  the  boundary. 
The  effect  of  complexions  on  strength  and  toughness  is  dif¬ 
ficult  to  generally  rationalize  on  the  basis  of  Eq.  (17)  be¬ 
cause  the  energies  of  the  two  resultant  surfaces  are 
difficult  to  predict. 


The  grain  boundary  self-diffusion  coefficient,  DGB , 
may  be  described  in  a  manner  analogous  to  the  bulk 
self-diffusion  coefficient  in  terms  of  the  atomic  jump  fre¬ 
quency,  v,  the  lattice  parameter,  aQ,  the  free  energy  of 
vacancy  formation,  A  6'/;/,,  the  activation  energy  barrier 
to  atomic  migration,  A  Gnub,  and  a  geometric  term,  g, 
as  follows: 

Dgb  =  gva20  exp(-A Gffi/RT)  exp(-A Gm,b/RT)  (18) 

Borisov  et  al.  [266]  first  postulated  that  the  free 
energy  of  activation  for  grain  boundary  diffusion, 
AGb(=AGfb  +  A G„ub),  is  proportional  to  the  corresponding 
free  energy  of  activation  for  diffusion  in  the  lattice,  AG/, 
minus  the  Gibb’s  free  energy  of  the  grain  boundary,  y. 

Although  this  rule  is  empirical,  it  provides  reasonable 
predictability  in  alloys  [267].  A  reduction  in  grain  bound¬ 
ary  energy  dominated  by  enthalpy  should  increase  the  acti¬ 
vation  energy  and  reduce  the  grain  boundary  diffusivity. 
However,  a  decrease  in  energy  dominated  by  increasing 
entropy  should  result  in  enhanced  diffusivity.  This  was 
recently  demonstrated  in  the  Ni-Bi  and  Cu-Bi  systems, 
where  an  increase  in  the  entropy  of  vacancy  formation, 
of  32  and  71  J  mol-1  K-1,  respectively,  was  associated  with 
thermally  induced  complexion  transitions  [268], 

A  chemically  induced  transition  in  segregation  level  and 
an  associated  transition  in  diffusivity  have  been  also 
reported  in  Cu-Bi  alloys  [162],  The  observed  discontinuous 
enhancement  in  diffusivity  as  a  function  of  composition  is 
likely  associated  with  a  complexion  transition  [35].  Indirect 
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Complexion  Transitions 


Congruent 


Faceting 


Non-Congruent 


Dissociation 


Grain  boundary  character  (R,  h )  remains 
constant  while  atomic  structure  and/or 
chemistry  of  the  grain  boundary  core 
changes.  During  a  congruent  transition,  one 
complexion  transforms  into  another. 


Grain  boundary  facets  into  two  new 
inclinations,  «,  and  n2,  while  the  misorientation 
R  remains  constant.  In  general,  the  atomic 
structure  and  chemistry  of  the  core  will  also 
change.  During  a  faceting  transition,  one 
complexion  decomposes  into  two. 

A  bulk  phase  completely  wets  the  original 
grain  boundary,  creating  two  new  grain 
boundaries  with  the  same  inclination  h  but 
different  misorientations  RI  and  R2.  During  a 
dissociation  transition,  one  complexion 
decomposes  into  two. 
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Dillon-Harmer 

Discrete 

Complexions 


Complexion  Categories 


1  interface  ^  2  interfaces 


Clean 


Monolayer 


Bilayer 


Trilayer 


Nanolayer 


Wetting 


~\ 


Dillon-Harmer  discrete  complexions  are  categorized  based  on  the  thickness  of  the  grain  boundary  core  converted  to  an 
integer  number  of  atomic  layers.  The  clean  complexion  has  very  little  or  no  adsorbed  solute  (and  hence  may  be  extrinsic 
or  intrinsic)  while  the  thicker  discrete  complexions  are  extrinsic  with  one  or  more  layers  enriched  in  solute.  Discrete 
complexions  may  have  varying  degrees  of  structural  and  chemical  order.  The  transition  from  a  single  complexion 
(interface)  to  two  complexions  (interfaces)  occurs  between  the  nanolayer  complexion  and  wetting. 


Intrinsic 

vs. 

Extrinsic 


Intrinsic  complexions  occur  in 
pure  single-component  systems 
and  are  free  of  adsorbed  solute, 
!  * !  * ! ' !  ’ !  ‘ *  but  may  have  different  atomistic 
configurations  and  different 
■  'I  grain  boundary  character. 

Intrinsic 


Extrinsic  complexions  contain 
adsorbed  solute  atoms  that 
influence  the  atomistic  structure, 
degree  of  crystalline  order,  and 
preferred  grain  boundary  character. 


Extrinsic 


Order 

vs. 

Disorder 


V 


Structurally 

Ordered 


Chemically 

Ordered 


A  structurally  ordered 
complexion  has  recognizable 
periodic  atomic  order. 


A  chemically  ordered 
complexion  has  a  periodic 
array  of  adsorbate  atoms. 


Structurally 

Disordered 


A  structurally  disordered 
complexion  has  very 
little  periodic  order. 


Chemically 

Disordered 


A  chemically  disordered 
complexion  has  adsorbate 
atoms  that  occupy  non¬ 
periodic  sites. 


J 


Fig.  26.  Schematic  diagrams  showing  different  types  of  complexion  transitions  and  different  ways  of  categorizing  complexions.  The  terminology  used  in 
these  diagrams  is  the  recommended  and  preferred  terminology  for  describing  complexions  and  transitions.  Terms  not  shown  in  these  diagrams  to  describe 
complexion  phenomena  in  more  detail  (e.g.  “roughening”)  can  be  used  within  the  framework  of  terminology  presented  here. 


indicators,  such  as  sintering  rate  or  grain  boundary  mobil¬ 
ity,  have  also  indicated  an  enhancement  in  mass  transport 
kinetics  associated  with  complexion  transitions  that 
increase  the  structural  or  chemical  width  of  the  grain 


boundary  [17,22,47,106,203,214],  The  simultaneous 
decrease  in  grain  boundary  energy  and  increase  in  diffu¬ 
sion-dependent  mass  transport  kinetics  suggests  that  these 
complexion  transitions  primarily  impact  the  entropic  con- 
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tribution  to  diffusivity.  For  ionic  crystals,  the  response  of 
the  anion  and  cation  sublattices  have  not  been  isolated  in 
any  given  system.  It  remains  unclear  how  anion  and  cation 
doping  might  affect  each  sublattice,  either  together  or 
separately. 

In  principle,  grain  boundary  complexion  transitions 
may  impact  many  other  grain  boundary  properties  as  well, 
including  electrical,  thermal,  optical,  chemical  and  mag¬ 
netic  properties.  The  role  of  complexions  in  affecting  these 
various  properties  has  not  been  thoroughly  investigated, 
although  there  is  some  evidence  that  complexion  transi¬ 
tions  are  correlated  with  marked  changes  in  these  proper¬ 
ties.  For  example,  intergranular  films  in  ZnO  introduce 
barriers  to  electrical  conductivity  that  produces  a  non-lin¬ 
ear  I-V  response  [269,270].  The  electrical  properties  of 
individual  grain  boundaries  in  MnZn  ferrites  containing 
no  solute  segregation,  CaO  segregation  and  intergranular 
films  have  been  characterized  [271,272].  The  grain  bound¬ 
aries  with  no  solute  segregation  displayed  the  highest  con¬ 
ductivity,  while  the  segregated  boundaries  exhibited  the 
lowest  conductivity.  The  resistance  of  X2Ru207  (X  =  Pb 
or  Bi)  samples  containing  intergranular  films  is  highly  sen¬ 
sitive  (three  orders  of  magnitude)  to  the  chemistry  of  those 
intergranular  films  [273],  Low-temperature  electrical  resis¬ 
tivity  measurements  in  Cu-Bi,  annealed  at  high  tempera¬ 
tures  and  then  quenched,  demonstrated  special  behavior 
under  conditions  in  which  a  dramatic  change  in  Bi  adsorp¬ 
tion  at  the  grain  boundary  occurred  [274],  suggesting  that 
grain  boundary  complexion  transitions  in  metallic  systems 
may  also  produce  measurable  changes  in  electrical  proper¬ 
ties.  One  might  envision  that  higher  conductivity  complex¬ 
ions  may  exist  in  certain  materials  that  could  enhance  bulk 
conductivity.  Grain  boundary  complexion  transitions 
occur  often  in  doped  titanates,  such  as  BaTi03  and  SrTi03, 
and  have  been  utilized  to  tailor  microstructural  develop¬ 
ment  [275,276].  SrTi03  grain  boundaries  with  intergranular 
films  may  display  varistor  behavior  [277].  However,  their 
effects  on  the  ferroelectric  and  dielectric  response  of  the 
grain  boundary  have  not  been  studied  in  detail. 

The  thermal  conductivity  of  grain  boundaries  is  also 
affected  by  complexion  transitions.  A  decrease  in  the  thick¬ 
ness  of  a  Y203-based  intergranular  film  in  AIN  was  corre¬ 
lated  with  an  increase  in  bulk  thermal  conductivity  (24%) 
[278].  However,  these  results  were  not  completely  isolated 
from  the  influence  of  second  phase  and  grain  size.  Inter¬ 
granular  films  have  also  been  associated  with  a  similar 
decrease  in  the  thermal  conductivity  of  SiC  [279,280],  Com¬ 
plexion  transitions  should  affect  color  centers  at  grain 
boundaries  as  they  influence  the  average  coordination 
and  bond  distances  of  atoms  at  the  interface.  The  effect 
on  local  bonding  has  been  demonstrated  through  X-ray 
absorption  fine  spectra  characterization  [167,281].  Optical 
mapping  of  complexions  could  be  an  ideal  approach  to 
characterizing  their  distributions  in  bulk  microstructures. 
Related  techniques  have  been  utilized  to  map  grain  bound¬ 
ary  networks  in  three  dimensions  [282,283],  With  the  pro¬ 
liferation  of  nanostructured  materials  containing  larger 


interfacial  volume  fractions,  the  need  to  understand  the 
effect  of  complexion  transitions  on  interfacial  properties 
will  continue  to  grow.  Investigating  interfacial  structure- 
property  relationships  in  the  context  of  complexions  pro¬ 
vides  ample  future  research  opportunities. 

5.2.  Microstructural  development 

Grain  boundary  complexion  transitions  have  received 
the  most  attention  as  a  result  of  their  influence  on  micro- 
structural  development.  In  fact,  the  first  direct  observation 
of  intergranular  films  was  motivated  by  their  role  in  dra¬ 
matically  influencing  the  sintering  behavior  of  Si3N4 
[221,222],  Subsequently,  complexion  transitions  were 
implicated  in  activated  sintering,  activated  grain  growth 
or  abnormal  grain  growth  in  A1203  [17,257,284-286], 
BaTi03  [276,287],  Mo  [204,205],  SiC  [288,289],  SrTi03 
[275,290],  W[106,207],  Y203  [33,291]  and  ZnO 

[47,292,293],  Higher  entropy  complexions  (i.e.  complexions 
that  are  thermodynamically  stable  at  higher  temperatures 
at  fixed  T,  P  and  composition)  are  known  to  exist  in  sys¬ 
tems  such  as  Al  [265],  AIN  [278,294],  Cu  [35,295-297],  Ni 
[48],  Si  [34],  Ti02  [36],  X2Ru207  (X  =  Pb  or  Bi)  [273], 
ZnMn  ferrite  [271,272]  and  Zr02  [298],  In  general,  higher 
entropy  complexions  are  associated  with  increased  rates 
of  mass  transport  and  high  diffusivity,  and  therefore  are 
likely  to  affect  abnormal  grain  growth,  activated  grain 
growth  or  activated  sintering  exhibited  by  those  systems. 

While  high  diffusivity  complexions  will  generally  reduce 
the  characteristic  temperature  at  which  sintering  becomes 
practical,  in  general  they  will  not  necessarily  enhance  den- 
sification  because  similar  surface  complexions  -  which  may 
also  have  increased  diffusivity  -  have  been  observed  in  most 
systems  exhibiting  higher  entropy  grain  boundary  com¬ 
plexions  (e.g.  A1203,  Ti02  and  ZnO)  [24,42,109,198,299], 
The  classic  competition  between  surface  diffusion  and  grain 
boundary  diffusion  in  promoting  coarsening  and  densifica- 
tion,  respectively,  complicates  a  simple  prediction  of  how 
interfacial  complexions  will  influence  the  early  stages  of 
sintering  in  a  particular  system  [300-302],  In  the  final  stages 
of  sintering,  the  grain  boundary  complexion  will  also  sig¬ 
nificantly  impact  the  competition  between  grain  growth 
and  densification.  The  influence  of  the  complexion  transi¬ 
tions  on  the  thermodynamic  driving  forces  [ygb  and  ys) 
must  be  considered  along  with  their  distribution  amongst 
grain  boundaries  and  surfaces  [29], 

Abnormal  grain  growth  results  from  a  difference  in  the 
average  grain  boundary  velocities  of  the  boundaries  sur¬ 
rounding  a  particular  grain  and  the  rest  of  the  neighboring 
population  [303].  The  difference  in  velocity  between  normal 
and  abnormal  grain  boundaries  is  typically  at  least  a  factor 
of  two.  This  velocity  difference  may  result  from  a  variety  of 
factors  related  to  the  local  driving  force  (e.g.  strain  energy 
or  a  non-linear  driving  force  associated  with  interface-con¬ 
trolled  processes)  or  the  grain  boundary  mobility  (e.g. 
localized  complexion  transitions  or  a  persistent  mobility 
advantage  for  a  specific  misorientated  grain  in  a  crystallo- 
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graphically  textured  material)  [304-307].  Here  we  distin¬ 
guish  between  abnormal  grain  growth,  in  which  a  small 
subset  of  the  population  experiences  rapid  growth,  and 
activated  grain  growth,  in  which  a  large  portion  of  the  pop¬ 
ulation  exhibits  a  discontinuous  change  in  the  average 
grain  growth  rate  as  a  function  of  annealing  time  or  tem¬ 
perature.  Systems  exhibiting  normal  grain  growth  have 
self-similar  grain  size  distributions,  systems  growing  abnor¬ 
mally  are  never  self-similar  and  systems  undergoing  acti¬ 
vated  grain  growth  are  not  self-similar  for  a  transient 
period. 

Complexions  promote  abnormal  grain  growth  when 
the  grain  boundaries  adjacent  to  a  small  fraction  of  the 
grains  (typically  <  10_3-10~5)  undergo  a  transition  that 
dramatically  increases  their  mobility  relative  to  the 
average  surrounding  population  [19,303,304,307],  This 
phenomenon  results  in  a  bimodal  (or  multimodal)  grain 
size  distribution  that  characterizes  abnormal  grain 
growth.  The  specific  boundaries  at  which  these  transitions 
occur  will  be  influenced  by  chemical  heterogeneities  and, 
more  importantly,  anisotropy  in  grain  boundary  energy 
and  structure.  While  this  transition  decreases  the  interfa¬ 
cial  energy,  and  therefore  the  thermodynamic  driving 
force  for  growth,  the  increase  in  mobility  typically  domi¬ 
nates  the  grain  growth  rate  [17,19,29],  These  rapidly 
growing  grains  become  abnormally  large  with  increasing 
time.  For  a  chemically  induced  complexion  transition, 
the  growing  grain  will  continuously  accumulate  solute  at 
the  grain  boundary.  This  excess  solute  might  precipitate 
in  the  lattice  behind  the  moving  boundary,  partition  to 
neighboring  boundaries,  precipitate  on  the  boundary  or 
thicken  the  grain  boundary  complexion.  In  BaTi03,  con¬ 
tinued  abnormal  grain  growth  has  been  associated  with 
thickening  of  intergranular  films,  but  not  to  an  extent  that 
can  account  for  all  of  the  excess  solute  [276],  In  A1203, 
precipitation  of  excess  solute  has  been  observed  behind 
rapidly  migrating  abnormal  grain  boundaries.  However, 
when  large  abnormal  grains  impinge  in  this  system,  grain 
boundary-precipitate  depinning  becomes  more  difficult 
and  precipitates  can  persist  at  the  boundaries  [21],  Inves¬ 
tigations  of  solute  redistribution  to  neighboring  grain 
boundaries  have  been  limited.  Such  redistribution  is  nec¬ 
essary  for  the  observation  that  a  particular  complexion 
type  can  seemingly  surround  an  entire  grain,  in  spite  of 
the  fact  that  it  likely  nucleated  at  an  individual  boundary, 
triple  line  or  quad  junction  [28],  Continued  accumulation 
of  solute  may  partially  explain  why  abnormal  grain 
growth  persists  even  as  the  misorientation  of  the  bound¬ 
ary  varies  during  grain  growth.  However,  hysteresis  asso¬ 
ciated  with  complexion  transitions  may  play  a  role  in 
stabilizing  either  or  both  the  normal  grain  population 
and  the  abnormal  grain  population  [108,308], 

When  grain  boundary  complexion  transitions  occur  at  a 
larger  fraction  of  grain  boundaries  (typically  >  10— 3— 10  3), 
almost  the  entire  population  rapidly  re-establishes  a  uni- 
rnodal  grain  size  distribution  [20,33].  During  this  transi¬ 
tion,  the  size  distribution  may  deviate  significantly  from 


log-normal  behavior,  but  may  re-establish  such  a  distribu¬ 
tion  with  continued  annealing.  Such  behavior  has  been 
observed  in  Y203,  where  a  10  °C  change  in  annealing  tem¬ 
perature  changes  the  average  grain  size  of  a  log-normal  dis¬ 
tribution  by  an  order  of  magnitude  [33],  This  transition 
may  also  be  captured  as  a  function  of  annealing  time. 
These  complexion  transitions  may  not  be  obvious  absent 
a  detailed  investigation  over  a  broad  range  of  annealing 
times  and  temperatures,  with  data  collected  at  reasonably 
small  time  and  temperature  intervals.  Similar  activated 
grain  growth  behavior  has  been  observed  in  Si02-doped 
A1203  [20],  It  has  also  been  observed  in  A1203  that,  when 
abnormal  grains  impinge  in  the  presence  of  high  dopant 
concentration,  bulk  second-phase  pinning  will  reduce  the 
growth  rate  of  these  boundaries  [21].  This  effect  previously 
suggested  that  inconsistent  results  existed  with  regard  to 
whether  intergranular  films  enhance  or  suppress  grain 
boundary  mobility. 

The  above  discussion  highlights  the  importance  of 
understanding  how  complexion  transitions  occur  spatially, 
temporally  and  as  a  function  of  temperature.  In  A1203,  the 
number  density  of  abnormal  grains  increases  linearly  with 
grain  size  at  a  fixed  temperature  and  exponentially  with 
increasing  temperature  [19].  A  qualitative  correlation  exists 
between  the  relative  grain  boundary  anisotropy  in  the  sys¬ 
tem  and  the  temperature  dependence  of  the  complexion 
transition  rate  [29],  The  grain  size  dependence  may  result 
from  the  fact  that  new  grain  boundaries  are  formed  ran¬ 
domly  during  grain  growth  [19], 

Grain  growth  suppression  has  attracted  significant 
interest  in  recent  years  for  applications  in  nanostructure 
stabilization  [309-313].  Complexion  transitions  should 
be  considered  as  important  microstructural  phenomena 
in  this  context.  Preventing  complexion  transitions  that 
enhance  diffusivity  might  be  an  effective  method  of  sup¬ 
pressing  diffusion-controlled  grain  growth  processes. 
However,  the  reduction  in  grain  boundary  energy  associ¬ 
ated  with  complexion  transitions  may  also  reduce  the 
driving  force  for  grain  growth.  In  fact,  a  reduced  driving 
force  for  grain  growth  achieved  by  lowering  the  grain 
boundary  energy  via  judicious  selection  of  dopants  is  cur¬ 
rently  believed  to  be  the  primary  contributor  to  nanocrys¬ 
talline  thermal  stability.  The  concept  that  grain  growth  is 
suppressed  by  the  stabilization  of  complexions  with  dra¬ 
matically  low  mobility  (rather  than  low  grain  boundary 
energy)  is  a  competing  explanation  that  remains  to  be 
explored. 

Strategies  to  drive  the  grain  boundary  energy  toward 
zero  have  produced  some  success  in  suppressing  grain 
growth.  Recently,  a  significant  discontinuous  decrease  in 
grain  boundary  mobility  was  observed  with  increasing  tem¬ 
perature  in  SrTi03  [314],  The  exact  mechanism  remains 
controversial,  although  the  experimental  results  indicate 
that  the  phenomenon  relates  to  a  grain  boundary  complex¬ 
ion  transition  in  which  the  atomic  structure  and  chemistry 
of  the  grain  boundary  changes.  Abnormal  grain  growth  in 
SrTi03  and  BaTi03  without  intergranular  films  is  difficult 


P.R.  Cantwell  et  al.  I  Acta  Materialia  62  (2014)  1-48 


43 


to  explain,  but  may  still  relate  to  transformations  in  the 
local  structure  and  chemistry  [27,314], 

Activated  sintering,  a  process  associated  with  significant 
sub-eutectic  enhancement  in  sintering  rate,  was  poorly 
understood  for  decades  [315-317],  Three  decades  ago,  the 
process  became  associated  with  grain  boundary  complex¬ 
ion  transitions  that  enhance  diffusivity  in  ceramics 
[221,292],  The  existence  of  intergranular  films  enables  the 
practical  pressureless  sintering  of  refractory  non-oxides 
such  as  AIN,  SiC  and  Si3N4  [278,292,318],  Activated  sinter¬ 
ing  promoted  by  complexions,  specifically  intergranular 
films,  has  been  noted  in  various  oxides  [47,273,319],  The 
concept  has  recently  been  extended  to  explain  activated  sin¬ 
tering  in  alloys  [22,32,203,320],  Complexions  might  be 
related  to  other  anomalous  sintering  behavior  that  has 
been  reported  in  the  literature,  where  intergranular  films 
are  not  observed  (e.g.  MgAl204  [321,322]  or  Y3A15012 
[323-325]),  but  more  investigation  is  required.  Recently, 
it  was  observed  in  systems  with  retrograde  solubility  that 
disordered  complexions  become  less  stable  with  increasing 
temperature  [204],  This  effect  reduces  the  diffusivity  and 
sintering  rate  with  increasing  temperature.  The  result  high¬ 
lights  the  role  of  the  thermodynamic  stability  of  interfacial 
complexions  in  introducing  complex  microstructural  evolu¬ 
tion  and  response. 

Grain  boundary  complexions  will  also  likely  affect  other 
grain-boundary-related  processing  phenomena,  such  as 
grain  boundary  sliding,  superplastic  deformation  and  wet¬ 
ting  [298,326,327],  The  role  of  complexions  in  affecting 
these  processes  has  not  been  thoroughly  explored.  The  sta¬ 
bility  of  complexions  may  be  sensitive  to  process  variables 
such  as  pressure  [296,328],  electric  field  or  magnetic  field 
that  have  not  been  investigated  in  detail. 

5.3.  Bulk  properties  and  behavior 

Thermodynamically  induced  grain  boundary  complex¬ 
ions  may  alter  fundamental  grain  boundary  properties, 
which  ultimately  affects  the  performance  of  a  bulk  poly¬ 
crystalline  aggregate.  While  the  role  of  complexions  in 
affecting  diffusional  transport  during  materials  processing 
is  reasonably  well  studied,  their  impact  on  diffusion-depen- 
dent  properties  such  as  oxidation,  creep  and  ionic  conduc¬ 
tivity  remains  somewhat  unexplored.  Creep  in  systems 
containing  intergranular  films  such  as  oxide-doped  silicon 
nitride  has  been  linked  to  enhanced  transport  in  these  films 
[298,329,330],  The  creep  and  oxidation  rates  in  alumina  are 
known  to  vary  by  several  orders  of  magnitude  at  a  single 
temperature  as  a  function  of  dopant  type,  dopant  concen¬ 
tration  and  processing  history  [37],  The  system  also  dis¬ 
plays  various  complexions  as  a  function  of  these  same 
parameters,  but  the  phenomena  have  not  been  directly  cor¬ 
related.  The  published  literature  often  focuses  on  condi¬ 
tions  that  produce  desirable  properties  (e.g.  low  oxidation 
rate  or  low  creep  rate),  which  may  explain  the  limited 
investigation  of  complexion-enhanced  oxidation  rate  or 
creep  rate.  In  order  to  identify  the  complexion  transitions 


responsible  for  improved  properties  such  as  low  oxidation 
and  creep  rates,  regimes  of  both  good  and  bad  properties 
must  be  explored,  with  special  attention  paid  to  the  divid¬ 
ing  line  between  these  regimes. 

Silica-rich  intergranular  films  have  been  observed  in  zir- 
conia  utilized  as  ion  conducting  electrolytes  [298,331], 
These  films  have  been  demonstrated  to  reduce  anion  con¬ 
duction  [331],  However,  similar  films  may  also  enhance  cat¬ 
ion  conduction.  For  example,  intergranular  films  enhance 
proton  conduction  in  LaP04  [332].  Similarly,  lithium  phos¬ 
phate-based  surficial  films  on  lithium  ion  battery  cathode 
particles  have  been  shown  to  enhance  battery  discharge 
rate  [299,333,334],  Enhanced  surface  diffusion  has  been 
invoked  as  a  possible  explanation,  although  other  theories 
have  also  been  proposed  [334], 

Complexion  transitions  have  been  linked  to  embrittle¬ 
ment  in  Al-Ga  [265],  Ni-Bi  [48],  SiC  [335,336]  and  Si3N4 
[337-339],  and  complexion  transitions  likely  explain  the 
anomalous  embrittlement  observed  in  A1203  [340],  Chemi¬ 
cally  and  thermally  induced  transitions,  analogous  to  com¬ 
plexion  transitions,  in  the  structure  of  dislocations,  triple 
lines  and  four-grain  junctions  may  also  occur.  Such  transi¬ 
tions  have  not  been  well  explored,  but  could  play  an  impor¬ 
tant  role  in  affecting  the  nucleation  and  propagation  of 
cracks.  In  many  systems,  avoiding  transitions  to  disordered 
complexions  may  be  a  critical  feature  of  effectively  process¬ 
ing  materials  with  desired  properties.  However,  embrittling 
complexions  have  been  exploited  to  control  crack  propaga¬ 
tion  around  large  abnormal  grains  in  a  scheme  that  utilizes 
those  grains  as  toughening  fibers  [336,338],  This  can  intro¬ 
duce  r-curve  toughening  behavior  into  intrinsically  brittle 
materials.  Cermets  such  as  WC-Co  contain  intergranular 
films  at  grain  boundaries  [223,341].  However,  their  impact 
on  the  relative  strength  and  toughness  of  these  materials  is 
not  fully  understood. 

6.  Conclusion 

The  term  “complexion”  has  been  introduced  to  differen¬ 
tiate  thermodynamically  stable  interfacial  states  of  matter 
from  bulk  phases  [13],  Although  the  concept  of  grain 
boundary  complexion  transitions  was  proposed  in  the  late 
1960s  [3],  sufficient  experimental  evidence  that  conclusively 
demonstrated  the  existence  of  such  transitions  did  not  exist 
until  the  1980s  [8].  In  recent  years,  complexion  transitions 
have  been  shown  to  dramatically  affect  microstructure 
development  and  properties  of  polycrystalline  materials. 
Grain  boundary  complexion  transitions  offer  mechanistic 
explanations  for  many  materials  phenomena,  such  as 
abnormal  grain  growth  [17],  activated  sintering  [47]  and 
liquid  metal  embrittlement  [48],  Many  other  properties, 
such  as  thermal,  electrical,  optical  and  magnetic  properties, 
may  also  be  affected  by  grain  boundary  complexion  transi¬ 
tions,  and  these  properties  should  be  explored  in  future 
studies. 

Much  of  the  progress  in  the  study  of  grain  boundary 
complexions  can  be  attributed  to  the  dramatic  enhance- 
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ment  in  atomic  resolution  capability  of  modem  S/TEMs 
that  occurred  with  the  commercial  introduction  and  wide 
availability  of  spherical  aberration  correction  during  the 
past  decade.  Nevertheless,  despite  the  enormous  progress 
that  has  been  made  in  the  field,  several  outstanding  chal¬ 
lenges  in  the  study  of  grain  boundary  complexions  remain. 
In  closing,  we  will  highlight  a  few  of  these  challenges  and 
suggest  avenues  for  future  research. 

From  a  materials  engineering  perspective,  one  of  the 
most  important  tasks  is  the  development  of  grain  bound¬ 
ary  complexion  diagrams  for  technologically  relevant 
materials  systems.  Efforts  in  this  area  are  already  under¬ 
way.  Lines  of  complexion  transition  have  been  determined 
by  theory  and  experiment,  and  overlaid  onto  phase  dia¬ 
grams  [22,31,32,106,161,342],  The  majority  of  this  work 
has  been  carried  out  on  binary  systems,  although  some 
results  for  ternary  systems  exist  [342],  Grain  boundary 
transitions  in  multi-component  systems  remain  largely 
unexplored,  even  though  they  represent  the  majority  of  engi¬ 
neering  materials.  Future  research  should  focus  on  deter¬ 
mining  complexion  transition  lines  in  model  binary  and 
ternary  systems  as  well  as  in  multi-component  engineering 
materials,  and  should  be  done  in  parallel  with  experiments 
that  determine  the  direct  relationship  between  each  type  of 
complexion  and  its  properties.  Such  experiments  are  time- 
consuming  and  challenging,  but  grain  boundary  complexion 
diagrams  overlaid  onto  bulk  phase  diagrams  will  be  most 
useful  if  it  is  known  how  each  complexion  transition  affects 
grain  boundary  properties  of  interest,  such  as  diffusivity, 
conductivity  and  mechanical  strength. 

Many  important  grain  boundary  complexions  are  only 
stable  at  elevated  temperatures,  but  most  experimental 
studies  to  date  have  relied  on  room-temperature  observa¬ 
tions  of  grain  boundary  complexions.  The  structure  and 
chemistry  of  the  high-temperature  complexions  is  presum¬ 
ably  “frozen”  in  place  by  rapid  quenching  of  the  specimen 
to  room  temperature.  However,  in  almost  every  case 
explored  to  date  it  is  unknown  whether  important  changes 
have  occurred  in  the  structure  and  chemistry  of  the  grain 
boundary  complexion  during  cooling.  In  situ  heating 
experiments  in  the  S/TEM  are  the  next  frontier  in  complex¬ 
ions  research,  and  may  prove  to  be  one  of  the  best  experi¬ 
mental  methods  for  constructing  complexion  diagrams.  In 
situ  experiments  allow  the  observation  of  grain  boundary 
complexions  under  nominally  equilibrium  conditions. 
Some  in  situ  observations  of  complexions  at  elevated  tem¬ 
peratures  have  already  been  made  [230],  but  this  area  of 
study  remains  largely  unexplored.  Issues  that  could  compli¬ 
cate  the  analysis  of  in  situ  heating  TEM  experiments 
include  contamination  from  TEM  specimen  preparation 
(especially  Ga  contamination  from  TEM  specimen  prepa¬ 
ration  with  the  FIB),  the  close  proximity  of  free  surfaces 
on  each  side  of  the  thin  TEM  specimen,  and  electron 
beam-specimen  interaction  effects.  The  relative  importance 
of  each  of  these  effects  in  disturbing  the  complexion  under 
observation  is  not  clear  and  must  be  determined  by  careful 
experimentation.  Furthermore,  the  magnitude  and  influ¬ 


ence  of  each  effect  will  depend  upon  the  properties  of  the 
material  and  the  complexion  under  study. 

It  is  believed  that  complexion  transitions  occur  via  a 
nucleation  and  growth  process.  However,  the  mechanism 
and  kinetics  of  the  process  are  not  yet  known.  Some  studies 
have  explored  the  kinetics  of  wetting  transitions  [343],  but 
analogous  experimental  studies  on  complexion  transitions 
have  not  yet  been  done.  A  recent  computer  simulation  of 
a  Cu  225(3 10)  grain  boundary  showed  that  a  structural 
complexion  transition  nucleates  at  the  intersection  of  the 
grain  boundary  and  the  free  surface  then  grows  inward 
along  the  boundary  away  from  the  surface  [39],  Research 
into  the  nucleation  and  growth  mechanisms  of  complex¬ 
ions  could  potentially  explain  the  spatial  inhomogeneity 
of  complexion  transitions  and  may  also  provide  novel 
insight  into  methods  of  inhibiting  complexion  transitions, 
which  could  be  useful,  for  example,  in  the  thermal  stabiliza¬ 
tion  of  nanocrystalline  materials. 

A  great  deal  of  progress  has  been  made  during  the  past 
few  decades  in  understanding  grain  boundary  complexion 
transitions.  These  transitions  can  strongly  influence  the 
properties  of  polycrystalline  materials  and  should  be  taken 
into  account  when  anomalous  or  unexpected  materials 
behavior  occurs.  The  advent  of  aberration-corrected 
S/TEM  and  other  advanced  characterization  techniques, 
combined  with  ever  more  powerful  and  realistic  computer 
simulations,  offers  great  promise  for  the  advancement  of 
grain  boundary  complexions  research.  By  better  understand¬ 
ing  and  controlling  grain  boundary  complexion  transitions,  it 
will  be  possible  to  advance  the  science  of  grain  boundary 
complexion  engineering  and  develop  new  materials  with 
novel  behavior  and  dramatically  improved  properties. 
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